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1.0  INTRODUCTION 


1 . 1  Summary 

Continuous  fiber  composite  laminates  are  known  to  undergo  a 
substantial  amount  of  complex  load-induced  damage  which  can  adversely 
afpect  component  performance  [1].  Therefore,  it  is  desirable  to  develop 
new  models  capable  of  accounting  for  the  effect  of  damage  on  material 
properties. 

This  report  documents  research  completed  during  the  second  year  of  a 
three  year  effort  under  AFOSR  grant  no.  AFOSR-84-0067  and  originally 
detailed  under  Texas  A&M  Research  Foundation  proposal  no.  RF- 8 3 14  and 
dated  October  1983.  The  objective  of  this  research  is  to  develop  an 
accurate  damage  model  for  predicting  strength  and  stiffness  of  continuous 
fiber  composite  media  subjected  to  fatigue  or  monotonic  loading  and  to 
verify  this  model  with  experimental  results  obtained  from  composite 
specimens  of  selected  geometry  and  makeup  to  be  described  herein. 

1 . 2  Statement  of  Work 

The  following  is  a  brief  summary  of  work  to  be  performed  under  the 
present  grant: 

1  )  develop  constitutive  equations  relating  stresses  to  strains  and 
damage  internal  state  variables  (ISV)  which  may  be  used  in  a  stress 
gradient  field; 

2)  develop  ISV  growth  laws  as  a  function  of  load  history  for  matrix 
cracking,  interlaminar  fracture,  etc.; 

3)  develop  finite  element  algorithms  capable  of  evaluating  ply 
properties  in  damaged  components. 

4)  perform  experiments  on  components  with  selected  stacking 
sequences  in  order  to  verify  the  model. 

2.0  RESEARCH  COMPLETED  TO  DATE 
2 . 1  Summary  of  Completed  Research 

The  following  research  has  been  completed  during  the  second  year: 

1 )  an  experimental  investigation  or  damage  mechanisms  in 

graphite/epoxy  laminates  (Appendix  6.1); 

2)  the  theoretical  relationship  between  ^racture  mechanics  and 
matrix  crack  damage  (Appendix  6.3);  and 

3)  a  theoretical  and  experimental  investigation  of  curved  matrix 
cracks  in  crossply  laminates  (Appendix  6.4). 

In  addition,  the  following  research  is  well  underway  at  the  end  of  the 
second  year: 


-1- 


h  .  b  )  ; 


1  '  damage  ’liimoi  1  i  nt;  curved  cracks  >n  urosaply  laminates  (Appendix 

J>  exper imentai  haracler >  z.n  >  on  or  interply  delami  nat  ’ ons ; 

j)  damage  modelling  o'"  interply  delami.nati  on  (Appendix  6.6);  and 

4)  a  study  of  model  accuracy  in  a  stress  gradient  field. 

2 . 2  Experimental  investigation  or  Damage  Mechanisms  in  Grapdite/Epoxy 
Laminates 


A  literature  survey  during  the  first  year  of  research  [1]  indicated 
that  no  experimental  data  were  currently  available  in  the  open  literature 
which  reported  all  of  the  data  necessary  to  characterize  the  damage  model 
developed  under  this  grant.  It  was  therefore  necessary  to  initiate  an 
experimental  effort  as  part  of  the  current  research.  This  research  has 
thus  far  resulted  in  one  M.S.  Thesis  (Appendix  6.1).  The  results  of  this 
experimental  effort  are  briefly  reviewed  in  this  section. 

It  was  decided  that  the  model  should  be  first  applied  '’or  a  single 
damage  mode.  The  damage  mode  of  interest  is  matrix  cracking.  Matrix 
cracks  are  cracks  in  any  ply  that  run  parallel  to  fibers  in  the  matrix 
material.  This  damage  mode  occurs  in  all  practical  continuous  fiber 
composite  structures.  This  phenomenon  was  investigated  experimentally  in 
seven  different  laminates  subjected  to  quasi-static,  stepwise  loading.  All 
the  laminates  were  cross-ply  laminates,  consisting  of  0°  plies  (fibers  in 
the  principle  load  direction)  and  90°  plies.  This  type  of  laminate  was 
chosen  to  promote  the  growth  of  matrix  cracks  and  to  suppress  other  damage 
modes.  A  comprehensive  data  base  is  essential  not  only  to  develop  the 
present  model  but  to  insure  that  the  experimental  results  are  useful  for 
future  research  needs.  Primary  attention  was  placed  'on  identifying  the 
mechanisms  or  initiation  and  growth  of  matrix  cracks.  Secondary  emphasis 
was  placed  on  observing  the  effect  of  cracking  on  material  response. 

A  qualitative  analysis  was  made  to  determine  the  significance  or  the 
damage  events  observed.  These  included  crack  density,  stiffness  reduction, 
crack  shapes,  development  op  crack  surface  area,  and  residual  strain.  An 
erf,ort  was  made  to  interpret  the  effects  or  stress  level,  specimen 
thickness,  geometry,  and  stacking  sequence  on  damage  rormation. 

In  this  experimental  program,  two  distinct  types  of  transverse  cracks 
were  observed  (See  Fig.  1.).  Typical  transverse  cracks  formed  at  regular 
intervals  along  the  length  o*'  the  specimen  and  generally  spanned  the  90° 
layer  perpendi cular  to  the  laminate  plane.  The  accepted  explanation  ror 
this  widely  observed  phenomenon  is  the  concept  of  shear  lag.  This  type  of 
transverse  crack  was  termed  a  "straight  crack."  The  second  type  of  crack 
had  many  features  that  indicated  that  it  possessed  a  totally  different 
growth  mechanism.  The  most  obvious  or  these  features  was  the  incidence 
angle  to  the  plane  or  the  laminate  at  which  the  crack  initially  propagated. 
This  crack  was  termed  an  "angle  crack.”  Several  important  observations  or 
this  phenomenon  are  enumerated  as  follows: 

1.  The  angle  always  formed  in  conjunction  witn  a  straight  crack. 

An  angle  crack  did  nut  form  as  an  isolated,  independent  crack. 


The  crack  was  always  oriented  toward  the  straight  crack. 


3.  The  angle  crack  formed  within  a  characteristic  distance  f'rom  a  ‘ 

straight  crack  or  another  angle  crack.  This  distance  was  pound  to  be  less 

than  or  equal  to  the  thickness  or  the  90°  ply  group.  . 

« 

4.  The  proportion  of  angled  cracks  to  straight  cracks  increased  ■ 

with  an  increase  in  the  90°  ply  thickness.  J 

1 

5.  Two  angle  cracks  at  opposite  interfaces  would  commonly  \ 

intersect  at  the  midplane  to  form  a  "curved  crack."  j 

6.  The  acoustic  emission  of  an  angle  crack  was  audibly  different 
from  that  o^  a  straight  crack. 

7.  Upon  unloading,  a  number  of  angled  cracks  could  be  seen  in 

intermediate  stages  of  growth.  This  was  not  the  case  with  straight  cracks.  j 

i 

Two  damage  phenomena  were  found  to  be  associated  with  the  cracks.  j 

Edge  delaminations  formed  as  bands  across  the  width  of  the  specimens.  In  1 

each  case,  the  delamination  propagated  along  the  angle  crack  interface.  j 

Longitudinal  splits  formed  over  the  delamination  areas  with  one  end  of  the  I 

split  terminating  at  the  angle  crack  interface. 

Several  observations  were  made  concerning  laminate  response.  A 
laminate  with  alternating  0°  and  90°  plies  developed  less  stiffness  loss 
than  a  similar  laminate  with  the  90°  plies  grouped  together  (See  Fig.  2.).  i 

Possible  causes  are  the  difference  in  90°  constraint  and  the  difference  in  I 

crack  opening  displacement.  Another  observation  was  that  the  strain  at 
which  first  ply  failure  occurred  was  inversely  related  to  the  number  of  90° 
plies  grouped  together. 

2.3  Relationship  Between  Fracture  Mechanics  and  Matrix  Crack  Damage 

In  order  to  demonstrate  the  potential  of  the  model  to  predict 
degraded  stiffness  components,  a  specialized  form  o'"  the  constitutive  model 
was  developed  for  the  case  of  transverse  matrix  cracks  in  90°  layers  or 
cross-ply  laminates.  For  this  special  case,  only  one  component  of  the 
Internal  state  variable  (ISV)  tensor  must  be  specified.  Specific  values  of 
the  ISV  corresponding  to  various  amounts  of  matrix  crack  damage  can  be 
determined  by  using  fracture  mechanics  to  relate  the  ISV  to  the  strain 
energy  release  rate  for  crack  extension  during  load-up.  The  relationship 
between  the  ISV  and  the  matrix  crack  strain  energy  release  rate  is 
developed  in  the  following  discussion  (Complete  details  are  given  in 
Appendix  6.3).  c 

The  local  energy  due  to  cracking,  u^,  is  related  to  the  ISV  ror 
matrix  cracks  as  follows  [See  Appendix  6.2]: 


u 


c 

L 


(1  ) 


where  is  a  material  constant,  e„  is  the  strain  In  the  90°  layer 

perpendicular  to  the  fibers  and  a.-,  is  che  ISV  for  straight-through  matrix 


FiK.  2.  Relationship  between  crack  surface  area  and 
stiffness  reduction  for  the  [02/902]s 
and  [0/90/0/90]  laminates 


racks  >n  yo°  layers.  Also,  the  local  energy  due  Lo  matrix  crack  growtn  < 
r’-'laled  Lo  tuo  strain  energy  r'eLease  rate,  G,  by 


a 


L 


(2) 


where  is  the  matrix  crack  surface  area  and  V  is  the  local  volume.  By 
equating  these  two  expressions  we  obtain  the  following  relationship  between 
the  ISV  and  G  for  crack  extension  during  load-up: 


1 


VL  e2 


load-up 


(3) 


For  a  matrix  crack  extending  completely  through  the  thickness  o' 
90°  layer,  the  strain  energy  release  rate  is  defined  as 


G 


(4) 


where  B  is  the  thickness  of  tne  90°  layer  and  G  is  the  strain  energy  in  the 
90°  layer.  It  can  be  shown  that  equation  (4)  reduces  to  the  following 
expression  for  a  crack  growing  in  the  presence  of  existing  cracks: 


G  =  B  [U  f (S/B)  ]  ,  (5) 

where  U  is  the  strain  energy  density  in  the  90°  layer  in  the  absence  of 
cracks  and  f(S/B)  is  a  '’unction  that  depends  on  the  crack  spacing,  S,  and 
the  layer  thickness.  Assuming  linear  elastic  behavior  of  the  90°  layer  in 
the  absence  of  cracks,  the  strain  energy  density  is  given  by 


U 

o 


(6) 


where  E  3  is  the  initial  undamaged  modulus  of  the  90°  layer  transverse  to 
the  '’ibers.  Because  of  the  complexity  of  the  interaction  between  cracks, 
the  '’unction  '’(S/B)  was  determined  from  experimental  tests  of  cross-ply 
laminates.  The  experimentally  based  expression  for  f(S/B)  together  with 
equations  (3),  (5)  and  (6)  are  used  to  specify  the  value  o'1  the  ISV  as  a 
'’unction  o'’  matrix  crack  surface  area  during  load-up. 

Although  it  is  possible  for  matrix  crack  surface  area  to  increase 
luring  unloading,  in  the  current  development  this  ef °ect  ’s  assumed  to  be 
negligible.  There'’ ore,  on  unloading  a0  depends  only  on  the  crack-closure 
displacement  and  would  go  to  zero  on  complete  crack  closure.  Assuming  that 


'he  crack-closure  d' splucoment  is  linear  w«  th  str'i's  i  m  *  1  the  matrix  :t 
.'dt,f'aeo  area  > a  constant ,  on  unloading  the  ISV  *  s  expressed  .o; 


unload  in.: 


-sere  c  is  a  constant  and  depends  only  on  the  matrix  crack  damage  state. 

Considering  a  tensile  test  with  a  load  and  unloading  cycle,  at  the 
instant  of  load  reversal  the  expressions  for  the  ISV  for  load-up  and 
unloading  must  be  equal.  Therefore,  by  equating  equations  (3)  and  (7)  at 
the  same  value  of  strain  yields  a  relationship  between  the  unloading 
constant,  c,  the  strain  energy  release  rate  and  the  associated  damage 
state.  This  relationship  is  given  by 


f  (S/B)dS.) 


.  (8) 


L  J22 


Since  the  local  volume  ror  matrix  crack  damage  is  the  thickness  of  1  ply 
multiplied  by  the  length  and  width  or  the  specimen  over  which  the  damage 
state  is  determined,  the  thickness  to  volume  ratio  in  equation  (8)  is  given 
by 


where  n  is  the  number  of  consecutive  90°  layers  and  the  damage  is  specified 
over  a  1.0"  width  and  1.0"  length.  Equation  (8)  reduces  to 


1  22 

C  =  2  ”  J 

22  0 


’(S/B)dS- 


It  hr.s  been  previously  shown  that  the  effective  stiffness  or  a  laminate  <s 
related  to  the  partial  derivative  of  the  ISV  with  respect  to  strain.  Using 
equation  (7)  above,  gives 


,  (ID 


which  5r»  a  constant  during  unloading.  Therefore,  equations  (10)  and  (11) 
are  used  to  predict  the  effective  axial  stiffness  (E  )  of  cross-ply 
laminates  as  a  '’unction  of  tiie  matrix  crack  damage  state. 


.’.4  Model  1  <  ng  A  ;c  1  1 1  SI  i  f f ness  Loss  Due  to  Matrix  Cracks 

As  described  in  detail  in  Append’ x  6.3,  equation  (11)  may  be  utilized 
to  predict  stirfness  loss  in  crossply  laminates  with  matrix  cracks.  by 
using  standard  laminate  analysis  techniques,  the  following  reduced 
stiffness  equations  may  be  derived: 


where  are  the  components  of  the  effective  laminate  stiffness,  n  is  the 
number  of  plies,  and  overbars  denote  quantities  measured  in  laminate 
coordinates. 

Equations  (12)  were  utilized  to  predict  the  damage  dependent  reduced 
axial  stiffness  op  severa^  crossply  laminates.  The  reduced  stiffnesses 
were  compared  to  experimental  results  obtained  for  graphite/epoxy  coupons 
composed  of  Hercules  AS4/3502.  A  typical  result  is  shown  in  Fig.  3.  The 
model  appears  to  be  accurate  in  predicting  stiffness  loss  ''or  the  crossply 
laminates  considered  herein.  Further  results  are  contained  in  Appendix 
6.3. 


2.5.  Modelling  Out-of-Plane  Stiffness  Loss  Due  to  Curved  Cracks 


In  Section  2.2  experimental  evidence  was  reviewed  which  indicated 
that  curved  matrix  cracks  can  occur  in  significant  quantities  in  crossply 
laminates.  As  a  measure  of  the  capability  of  the  model  to  predict 
stiffness  loss  components  other  than  the  axial  stiffness,  it  was  decided  to 
use  the  model  to  predict  out-of-plane  stiffness  loss  due  to  the  curved 
cracks.  Although  this  research  is  still  underway  at  this  time,  the 
theoretical  procedure  has  been  completed.  The  procedure  will  be  briefly 
reviewed  here.  Further  details  are  given  in  Appendix  6.5. 

Consider  now  a  local  volume  element  with  n  cracks  as  shown  in  Fig. 
4.  Equations  (1)  may  be  written  in  the  following 'form: 
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where  matrix  crac'  ing  is  designated  by  the  superscript  1.  Now  define 
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so  that  equation  (5)  may  be  wr’tten 
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Now,  since  a . ^  Is  a  second  order  tensor, 
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where  ,  are  the  direction  cosines  relating  the  coordinates  of  the  kth 
crack  top  the  laminate  coordinates.  Differentiating  (16)  with  respect  to 
the  midplane  strain  components  gives 
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Transforming  to  the  coordinates  of  the  crack  gives 
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However,  since  it  is  assumed  that  3  a?2^z22  is  the  on-*-y  non-zero  component, 
the  above  reduces  to 
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The  above  equations  may  be  utilized  to  obtain  the  last  term  in  reduced 
stiffness  equations  ( ^4 ) ,  where  it  is  assumed  that  .  a0  ,  /  c  , ,  ,,  is 
independent  of  crack  orientation  and  is  obtained  from  [0,90^0]  experimental 
data  [2]. 


fy> 


An  oxporUn-jnt.il  study  of  matrix  crack  damage  reveal od  that  the  curved 
cracks  formed  ,ipter  the  straight  cracks  and  '’oi lowed  a  repeatable  pattern 
op  location  and  orientation  relative  to  the  straight  cracks.  For  a 
[0,90,,]  laminate  the  average  location  of  the  curved  crack  relative  to  the 
straight  crack  was  approximately  two  ply  thicknesses  away.  Furthermore,  it 
was  found  that  the  curved  cracks  were  never  an  isolated  phenomenon. 
Therefore,  it  was  postulated  that  the  growth  mechanism  for  curved  cracks  is 
driven  by  the  stress  state  resulting  from  the  formation  of  the  straight 
crack.  The  cracks  were  found  to  exhibit  slow  stable  crack  growth  behavior 
through  the  thickness  unlike  the  brittle  fracture  behavior  of  the  straight 
cracks.  An  analytical  investigation  of  the  curved  crack  phenomenon  was 
performed  using  a  finite  element  model  of  a  crossply  laminate  with  straight 
cracks.  The  analytical  results  for  a  [0,90  ]  indicated  that  the  principal 
stress  and  strain  undergo  significant  rotation  as  well  as  attaining  a 
maximum  value  at  approximately  one  ply  thickness  away  from  the  cracks  (see 
Fig.  5.),  thus  substantiating  the  postulated  growth  mechanisms.  Growth 
mechanisms  for  both  types  of  cracks  were  experimentally  related  to  the 
number  of  consecutive  90°  plies,  the  ratio  of  0°  to  90°  plies,  and  the 
applied  stress  level.  Further  details  of  this  study  are  given  in  Appendix 
6.4. 

2.6  Modelling  Axial  Stiffness  Loss  Due  to  Interlaminar  Delamination 

As  described  in  detail  in  Appendix  6.6  the  axial  stiffness  loss  in 
crossply  laminates  with  both  matrix  cracks  and  interlaminar  delaminations 
is  given  by 

K  -  Ai i  -£<s,>k  u 

k^l 

where  A1  ^  is  the  undamaged  axial  stiffness,  the  second  term  denotes  the 
axial  stiffness  loss  due  to  matrix  cracking,  where  a  represents  the  ISV 
for  matrix  cracking,  and  the  last  term  corresponds  to  the  stiffness  loss 
due  to  interlaminar  delamlnation,  where  a-}  represents  the  ISV  for 
delamination.  In  deriving  equation  (20)  a  different  local  volume  must  be 
used  for  obtaining  the  locally  averaged  ISV  for  interlaminar  delamination 
than  for  matrix  cracks.  Here  the  entire  laminate  is  treated  as  the  local 
volume.  For  interlaminar  delaminaton  there  are  only  three  nonzero 
components  of  the  tensor  valued  ISV.  However,  for  predicting  axial 
stiffness  only  is  of  interest. 

The  ISV  for  interlaminar  delamination  is  related  to  the  surface  area 
following  a  similar  procedure  for  matrix  cracking.  Incorporating  fracture 
mechanics  and  thermodynamics  gives 


Gd  dS  ,  (21  ) 
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D 

where  h  is  the  total  ply  thickness,  V^  is  the  local  volume,  S?  is  the 
damage  surface  area,  and  is  the  energy  release  rate  associated  with 
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ABSTRACT 


An  Investigation  of  Damage  Accumulation  in 
Graphite/Epoxy  Laminates.  (August  1985) 

Robert  Gerald  Norvell,  B  .  S  .  ,  Texas  A&M  University 

Co-Chairs  of  Advisory  Committee:  Dr.  David  H.  Allen 

Dr.  Richard  A.  Schapery 

The  objective  of  this  investigation  has  been  to 
identify  the  mechanisms  of  initiation  and  growth  of  matrix 
cracks  in  graphite/epoxy  laminates  and  to  identify  the 
effect  of  matrix  cracking  on  material  response.  An 
extensive  experimental  data  base  was  produced  for  use  in 
the  development  of  a  damage  model  and  for  model  verifica¬ 
tion. 

An  as  yet  unreported  form  of  transverse  cracking  has 
been  observed.  Two  distinct  forms  of  transverse  cracks 
were  found,  each  clearly  having  its  own  mechanisms  of 
initiation  and  growth.  Subsequent  damage  modes  associated 
with  transverse  cracks,  such  as  longitudinal  splitting  and 
delamination,  also  developed  seperate  forms  corresponding 
to  the  transverse  crack  variations. 
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INTRODUCTION 


The  advantages  of  composite  materials  for  use  in 
aerospace  vehicles  are  numerous.  Properties  such  as  high 
strength- to-weight  ratio  and  reduced  part  count  make 
composite  laminates  desireable  replacements  for  metallic 
structures.  However,  most  applications  have  been  re¬ 
stricted  to  non-critical  structures  because  of  a  lack  of 
understanding  of  the  behavior  of  composite  laminates  as 
their  ultimate  strength  is  approached.  This  problem 
arises  due  to  the  complex  microstructure  of  the  composite 
which  lends  itself  to  complicated  forms  of  damage.  Damage 
is  defined  as  any  fracture  process.  Unlike  metals,  com¬ 
posite  laminates  are  subject  to  an  array  of  possible 
damage  modes.  One  damage  mode  may  interact  with  another 
damage  mode  or  may  contribute  to  the  growth  of  new  modes. 
All  thi3  leads  to  a  very  complex  damage  state  which  is 
dependent  on  stacking  sequence.  Because  of  the  major 
differences  in  damage  growth  between  composites  and 
metals,  linear  elastic  fracture  mechanics  of  non¬ 
interacting  cracks  and  other  standard  analytical  methods 
have  failed  when  applied  to  composites.  To  maximize  the 
use  of  composites,  it  is  of  utmost  importance  to  under¬ 
stand  the  mechanisms  of  damage  and  to  develop  models  to 
predict  strength  reduction  and  changes  in  material 

This  thes i s  follows  the  style  and  format  of  the  Jour¬ 
nal  of  Composite  Materials. 


properties  due  to  these  mechanisms. 

The  research  described  here  represents  the  experimen¬ 
tal  portion  of  a  damage  model  for  fiberous  composites 
currently  being  developed  at  Texas  A&M  University  [1,2]. 
The  objective  of  this  research  has  been  to  identify  the 
parameters  controlling  damage  accumulation  in  laminated 
composites.  The  damage  mode  of  interest  is  matrix  crack¬ 
ing.  Matrix  cracks  are  cracks  in  any  ply  that  run  paral¬ 
lel  to  fibers  in  the  matrix  material.  This  damage  mode 
occurs  in  all  practical  continuous  fiber  composite 
structures.  This  phenomenon  has  been  investigated 
experimentally  in  seven  different  laminates  subjected  to 
quasi-static,  stepwise  loading.  All  the  laminates  are 
cross-ply  laminates,  consisting  of  0°  plies  (fibers  in  the 
principle  load  direction)  and  90°  plies.  This  type  of 
laminate  has  been  chosen  to  promote  the  growth  of  matrix 
cracks  and  to  suppress  other  damage  modes.  A  comprehen¬ 
sive  data  base  is  essential  not  only  to  develop  the 
present  model  but  to  insure  that  the  experimental  results 
are  useful  for  future  research  needs.  Primary  attention 
has  been  placed  on  identifing  the  mechanisms  of  initiation 
and  growth  of  matrix  cracks.  Secondary  emphasis  has  been 
placed  on  observing  the  effect  of  cracking  on  material 
response. 

A  qualitative  analysis  has  been  made  to  determine  the 
significance  of  the  damage  events  observed.  These  include 


crack  density,  stiffness  reduction,  crack  shapes,  develop¬ 
ment  of  crack  surface  area,  and  residual  strain.  An 
effort  has  been  made  to  interpret  the  effects  of  stress 
level,  specimen  thickness,  geometry,  and  stacking  sequence 


on  damage  formation. 


LITERATURE  SURVEY 


Current  literature  concerning  damage  growth  mecha 
nisms  and  damage  accumulation  in  advanced  fiberoua  com¬ 
posites  is  extensive  and  diverse.  Work3  concerning  matrix 
cracking,  both  analytical  and  experimental,  are  discussed 
in  the  following  paragraphs. 

One  of  the  earliest  investigations  of  matrix  cracking 

was  reported  by  Vasilev,  et  al .  [3]  in  which  a  closed  form 

solution  was  derived  for  simplified  governing  equations 

for  a  [0/90]  laminate  with  a  transverse  crack  in  the  90° 

3 

plies  subjected  to  a  load  in  the  0°  direction.  This 
analysis  predicts  distributions  of  the  two  in-plane  normal 
stresses  and  the  interlaminar  3hear  stress  as  a  function 
of  distance  from  the  crack  in  both  the  0°  and  90°  plies. 

Tsai  and  Hahn  [4,5]  developed  a  simplified  elasticity 
model  of  a  cracked  90°  ply  within  a  laminate  which  was 
used  to  predict  matrix  crack  spacing  and  the  overall 
stress-strain  behavior  of  a  C  0 / 9  0  ]  s  glass/epoxy  laminate. 
Their  model  involved  a  gradual  failure  mode,  employing  a 
secant  modulus  which  is  reduced  through  a  parameter  repre¬ 
senting  the  degree  of  degradation  occurring  in  the  90° 
plies  as  a  function  of  the  applied  stress  when  subjected 
to  uniaxial  tension.  This  is  significantly  different  from 
the  method  proposed  by  Jones  [6]  in  which  the  entire 


lamina  stiffness  matrix  is  set  to  zero  with  the  onset  of 
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cracking. 

Kim  &  Hahn  [7]  have  shown  the  importance  of  environ¬ 
mental  conditioning  on  composite  laminates.  They  investi¬ 
gated  the  effect  of  residual  stress  on  the  onset  of  matrix 
cracking  which  they  refer  to  as  first  ply-failure  (FPF). 
Variation  of  residual  stresses  was  accomplished  by  mois¬ 
ture  absorbsion.  They  noted  that  in  fiber  dominated 
laminates  (containing  0°  fibers)  the  absorbed  moisture 
increases  both  the  FPF  stress  and  the  ultimate  tensile 
strength.  Herakovich  and  coworkers  [8]  investigated  the 
characteristics  of  thermally  induced  transverse  cracks  in 
graphite/epoxy  cross-ply  and  quasi- isotropic  laminates, 
both  experimentally  and  analytically.  The  state  of  stress 
in  the  vicinity  of  a  transverse  crack  and  the  influence  of 
transverse  cracking  on  the  laminate  coefficient  of  thermal 
expansion  was  predicted  using  a  generalized  plane  strain 
finite  element  analysis  and  a  modified  shear  lag  analysis. 

Garrett  and  Bailey  [9]  investigated  the  thickness 
effect  of  90°  plies  in  cross-ply  laminates.  Using  speci¬ 
mens  of  glass  fibers  and  polyester  resin,  they  varied  the 
transverse  ply  thickness  between  .75  mm  and  3  mm.  They 
observed  that  for  a  given  strain,  a  higher  crack  density 
occurred  as  the  number  of  transverse  layers  decreased. 
Cracks  were  seen  on  the  edge  of  the  specimen  running  at 
45°  to  the  main  transverse  cracks.  These  cracks  formed  at 
high  strains  and  failure  occurred  soon  thereafter.  No 
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attempt  was  made  to  explain  these  angled  cracks.  A  simple 
theory  was  advanced  to  predict  the  crack  spacing  for  a 
given  applied  stress  and  transverse  ply  thickness  /hen 
longitudinal  and  transverse  plies  are  elastically  bonded. 
This  theory  is  based  on  earlier  work  by  Aveston  and  Kelly 

Cio]. 

Flaggs  and  Kural  [11]  reported  experimental  work  in 
which  they  investigated  the  uniaxial  tensile  loads  at 
which  transverse  cracking  initiated  in  the  90°  laminae  of 
[±0/9On]g  T300/934  laminates.  The  transverse  strengths  of 
the  90°  laminae  corresponding  to  these  loads  were  then 
calculated  based  upon  lamination  theory  predictions  in¬ 
cluding  effects  of  residual  thermal  stresses.  The  stress 
at  which  transverse  cracking  initiated  was  found  to  be  a 
function  of  both  laminae  thickness  and  the  orientation  of 
adjacent  ±0  laminae. 

A.S.D.  Wang  and  coworkers  [12-16]  analyzed  the  growth 
mechanisms  of  transverse  cracking  and  delamination  with 
the  application  of  linear  elastic  fracture  mechanics  and 
finite  element  calculations  of  the  strain  energy  release 
rate.  A  series  of  gr a  phi t e / e po xy  laminates  in  the  form  of 
[±25/90n]3>  n- 1 /2 , 1 , 2 , 3 . ^ . 6 , 8 ,  were  examined  both 
analytically  and  experimentally.  To  predict  transverse 
crack  spacing,  the  strain  energy  release  rate  associated 
with  the  formation  of  a  single  transverse  crack  in  the  90° 


ply  was  calculated  using  a  crack  closure  method.  They 


then  postulated  that  saturation  spacing  is  that  spacing  at 
which  the  strain  energy  release  rate  associated  with  the 
formation  of  a  second  crack  is  equal  to  that  associated 
with  the  formation  of  the  first  crack.  They  assumed  that 
the  values  of  critical  strain  energy  release  rate  associ¬ 
ated  with  transverse  fracture  and  interply  delamination 
were  the  same.  Their  experimental  results  showed  that 
transverse  cracks  did  not  occur  for  n-1/2  and  n-1  lami¬ 
nates  until  the  onset  of  delamination.  For  n-2  and 
greater,  the  transverse  crack  density  increased  with  90° 
ply  thickness. 

Relfsnider  [17]  developed  a  one-dimensional  shear-lag 
model  which  can  approximate  normal  stress  in  the  load 
direction  in  the  various  plies  of  a  laminate.  He  contends 
that  shear  deformations  in  any  given  ply  are  restricted  to 
a  thin  area  in  the  region  of  ply  interfaces.  In  the 
model,  this  thin  region  is  assumed  to  transfer  all  the 
loads  between  plies  through  shear.  The  model  predicts  the 
longitudinal  stress  field  in  the  region  of  a  transverse 
crack,  and  with  a  critical  stress  criterion,  the  model 
predicts  the  crack  spacing  in  the  laminate. 

Talug  [18]  formulated  a  finite  difference  solution  to 
determine  how  transverse  cracking  in  off-axis  plies  af¬ 
fects  the  stiffness  of  the  overall  laminate.  This  tech¬ 
nique  yields  a  full-field  solution  where  all  six  stress 
components  are  determined  throughout  the  laminate  in  the 


presence  of  a  crack.  With  this  information,  the  distance 
from  the  crack  over  which  the  stresses  are  redistributed 
to  reach  laminate  analysis  values  can  be  determined. 

Also,  the  out  of  plane  stresses,  which  may  initiate  delam¬ 
inations,  are  determined. 

Reifsnider,  et  al .  [19,20]  introduced  the  concept  of 
a  "characteristic  damage  state  (CDS),"  a  condition  in 
which  a  saturation  of  matrix  cracking  i3  ultimately 
achieved  in  both  static  and  cyclic  loading  and  for  which  a 
characteristic  spacing  of  cracks  is  observed.  The  CDS  is 
said  to  be  independent  of  load  history,  and  depends  only 
upon  the  material,  geometry,  and  stacking  sequence.  From 
a  mechanics  standpoint,  Reifsnider  and  Masters  [21]  pro¬ 
pose  the  CDS  has  the  same  significance  as  the  single  crack 
for  homogeneous  materials  with  regard  to  the  well-defined 
physical  state  from  which  the  fracture  event  develops. 

They  also  point  out  that  the  CDS  is  different  in  that  the 
final  fracture  does  not  occur  by  the  further  development 
of  that  state  as  with  the  single  crack  situation.  Bader, 
et  al .  [22]  also  reported  on  the  formation  of  a  character¬ 
istic  spacing  of  transverse  cracks  in  angle-ply  laminates. 
Highsmith  and  Reifsnider  [23]  have  examined  changes  in 
stiffness  due  to  matrix  cracks.  They  proposed  a  model  for 
treating  the  problem  analytically  wherein  elements  of  the 
lamina  stiffness  tensor  are  systematically  reduced  in 
relation  to  crack  density  and  the  laminate  stiffness  is 


is 
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then  determined  from  a  laminated  plate  theory. 

Kelly  [24]  investigated  transverse  cracking  and  in¬ 
cluded  the  effects  of  curing  stress.  He  hy;othesized  the 
existence  of  longitudinal  splits  caused  by  a  mismatch  in 
Poisson's  ratio  between  plies.  He  suggested  that  the 
splitting  would  be  complete  when  a  saturation  of  trans¬ 
verse  cracking  occurred.  Laws,  Dvorak,  and  Hejazi  [25] 
developed  constitutive  equations  for  fibrous  composites 
which  contain  groups  of  longitudinal  slits.  They  inves¬ 
tigated  the  case  in  which  the  cracks  and  fiber  diameter 
were  of  similar  size  and  and  the  special  case  of  small 
diameter  fibers  and  large  cracks.  The  overall  elastic 
moduli  and  compliance  components  of  the  cracked  composite 
were  obtained  using  a  variation  of  the  self-consistent 
method . 

Stinchcomb,  et  al .  [26]  studied  the  effects  of  ply 
constraint  on  damage  development  in  composite  materials. 
They  determined  that  through-the-thi ckness  constraints 
control  the  pattern  and  spacing  of  cracks  which  in  turn 
determines  the  state  of  stress  and  state  of  strength  in 
the  damaged  laminate.  They  also  determined  that  con¬ 
straint  situations  that  produce  greatest  static  strength 
do  not  minimize  the  extent  of  damage  that  develops  during 
loading.  Likewise,  minimum  damage  situations  do  not  cor¬ 
respond  to  maximum  strength  cases. 

Nuismer  et  al .  [27-28]  introduced  a  continuum  theory 


’  v, 


for  the  constitutive  behavior  of  a  damaged  composite  ply. 

A  model  of  the  growth  of  a  matrix  crack  in  a  90°  ply 
subjected  to  transverse  tension  was  deveioped  and  incor¬ 
porated  into  the  theory.  This  model  is  based  on  fracture 
mechanics  and  an  approximate  stress  analysis.  The  model 
predictions  indicate  that,  after  damage  begins,  the  ply 
constitutive  behavior  is  laminate  dependent. 

Internal  state  variable  theories  have  been  recently 
incorporated  into  cumulative  damage  models.  These  models 
are  based  on  the  notion  that  one  or  more  internal  state 
varibles  can  be  used  to  describe  the  state  of  damage 
within  a  composite  laminate.  These  variables  may  be 
scalar,  vector,  or  tensor  valued.  Talreja  [29]  proposes 
that  damage  be  modeled  by  a  set  of  independent  vector 
fields  based  on  continuum  mechanics.  The  variables  would 
be  proportional  to  the  crack  density,  crack  surface  area, 
and  a  shape  parameter.  Allen,  et  al.  [1,2]  propose  a 
model  using  tensor  valued  internal  state  variables  re¬ 
flected  through  the  local  constitutive  equations.  A  set 
of  local  functions  describing  the  energy  of  cracking  is 
suggested  for  characterizing  the  damage  variable.  The 
local  variables  are  globally  averaged  over  the  entire 
laminate  so  that  the  effects  of  damage  are  averaged  over  a 
local  volume  element. 


EXPERIMENTAL  TECHNIQUE 


Laminate  Selection 

The  principle  objective  of  this  investigation  was  to 
observe  the  effects  of  matrix  cracking  on  the  stiffness  of 
a  laminate.  The  second  objective  was  to  glean  as  much 
information  as  possible  about  the  damage  mechanisms  of 
matrix  cracking.  To  meet  these  objectives,  several  fac¬ 
tors  were  considered  in  selection  of  laminates  to  be 
studied.  Laminate  sequences  that  included  90°  plies  were 
an  obvious  choice  since  matrix  cracking  is  the  predominant 
form  of  damage  in  this  lay-up.  To  this  end,  cross-ply 
laminates  consisting  of  only  0°  and  90°  plies  were  desir¬ 
able.  Another  consideration  was  to  preclude  other  damage 
modes.  The  last  consideration  was  to  use  a  realistic 
material  system.  In  this  case,  a  graphite/epoxy  system 
was  desired. 

In  graphite/epoxy  systems,  the  ratio  of  longitudinal 
modulus  to  transverse  modulus  (E,/E2)  is  very  large  as 
compared  to  a  glass/epoxy  system.  The  load  carried  by  90° 
plies  is  small,  and  even  with  total  degradation  of  the  90° 
layers,  the  decrease  in  modulus  is  small.  For  example,  a 
[0/90]  laminate  would  have  a  stiffness  loss  of  about  6 
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percent.  This  might  suggest  that  a  ±0°  combination  might 
be  substituted  for  the  0°  plies  to  decrease  the  fiber 
dominance.  Unfortunately  this  greatly  increases  the 


tensile  interlaminar  stresses  and,  thus,  the  tendency  for 
delamina t i on .  With  this  in  mind,  the  decision  was  made  to 
use  the  laminates  with  only  0°  and  90°  combinations. 

Seven  laminates  were  chosen  for  this  investigation. 
These  are  listed  in  Table  1.  The  laminate  types  are  iden¬ 
tified  by  letter  designations.  The  [0/90  ]  series  was 

n  3 

used  to  investigate  the  effect  of  the  thickness  of  the  90° 
layer  on  laminate  response  and  damage  growth.  The  [0/90] 

3 

and  [02/902]  laminates  were  included  for  scaling  compari- 

son.  The  [0/90/0]  and  [0/90/0/90]  were  studied  for  the 

s  s 

effect  of  0°  ply  constraint  on  crack  growth  in  the  90° 
plies.  There  were  numerous  other  points  of  comparison  for 
this  laminate  group. 

Two  other  laminates  were  tested.  A  [0,]t  (Type  Z) 
laminate  and  a  [90#]t  (Type  N)  laminate  were  used  to 
determine  laminate  properties. 

Specimen  Fabrication 

The  material  used  throughout  this  investigation  was 
AS4/3502  graphite/epoxy  manufactured  by  Hercules  Corp. 

The  material  was  supplied  as  a  preimpregnated  tape  in  12 
inch  wide  rolls.  The  material  was  stored  at  0°F  until 
ready  to  be  used.  Table  2  lists  the  various  properties  of 
the  material.  Laminates  of  12  inch  by  12  inch  size  were 
fabricated  locally  using  an  air  cavity,  heated  platen 
press  with  a  microprocessor  to  control  heating  and  cooling 


Table  2.  Material  Properties  for  Hercules  AS4/3502 


Lamina  Properties 

Longitudinal  Modulus,  EtI  21.0  x  10s  ±  2.0%  psi 
Transverse  Modulus,  EJ2  1.39  x  10s  ±  2.1%  psi 
Shear  Modulus,  G12‘  .694  x  10s  psi 


Poisson's  Ratio 


Longitudinal  Strength, 

F 

tu  t 

326000 

± 

3.5% 

psi 

Transverse  Strength, 

F 

tu  j 

1  1085 

± 

9.8% 

psi 

Long.  Failure  Strain, 

€tUj 

.0144 

± 

4.6% 

in/in 

Tran.  Failure  Strain, 

etUj 

.00773 

+ 

6.7% 

in/  i  n 

16  4  g/ra1  2 
.58  ±  .03 
35  ±  3 


1 Theore  t ical 

2Hercules  Inc.  Specifications 


rates,  vacuum  and  cavity  pressure.  The  laminates  were 
cured  following  the  cure  cycle  in  Figure  1.  Each  stacked 
laminate  was  weighed  before  ani  after  the  curing  cycle. 

In  this  way  the  amoun.t  of  resin  bled  from  the  plate  was 
monitored.  This  information,  in  conjunction  with  resin 
density  and  resin  content,  was  used  to  calculate  the  fiber 
volume  fraction  after  cure.  These  data  are  presented  in 
Table  3*  By  placing  a  .025  inch  thick  sheet  of  milled 
teflon  on  each  side  of  the  stacked  laminate  before  curing, 
a  fine  matte  finish  was  left  on  the  surface  of  the  lami 
laminate.  In  this  way  the  cross-hatch  pattern  of  peel-ply 
or  bleed  cloth  was  avoided.  Likewise,  the  wrinkles  or  air 
tracks  sometimes  left  by  teflon  release  film  were  also 
avoided.  This  was  important  so  as  not  to  create  stress 
concentrations  or  patterns  on  the  X-ray  exposures.  Each 
laminated  plate  was  inspected  by  visual  observation  for 
general  condition  and  surface  anomalies,  and  by  ultrasonic 
C-Scan  for  void3,  inclusions,  and  delaminations. 

A  typical  specimen  was  identified  by  its  laminate 
type  and  the  position  it  was  cut  from  the  laminated  plate. 
For  example,  specimen  C-3  was  the  third  specimen  cut  from 
the  [0/90,]  laminate. 

3 

Specimens  were  cut  from  the  plate  using  a  Micro-Matic 
Precision  Slicing  and  Dicing  Machine.  The  blade  was  a  6 
inch  diameter,  .05  inch  thick,  diamond-edge  cut-off  wheel. 
The  wheel  was  operated  at  2850  rpm  with  a  feed  rate  of  1.3 


Table  3 
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Post-cure  Resin  Content 


LAMINATE 
WEIGHT  g, 
UNCURED 


LAMINATE 
WEIGHT  g 
CURED 


RESIN 

LOST 

g 


RESIN 
CONTENT, 
%  WT. 


in/tnin  and  was  cooled  with  ethylene  glycol.  The  nominal 
dimensions  for  all  specimens  were  1  inch  in  width  and  11 
inches  in  length.  Widtl  measurements  were  made  with  a 
dial  caliper  at  three  locations  in  the  one  inch  gage 
length  and  averaged.  Thickness  measurements  were  made 
with  a  .0001  inch  resolution  micrometer  at  nine  locations 
in  the  one  inch  gage  length  and  averaged. 

The  0°  unidirectional  specimens  were  provided  with  2 
inch  end  tabs  made  from  woven  cross-ply  glass/epoxy 
material.  Strips  of  tabbing  material  were  tapered  at  15° 
and  bonded  to  the  laminated  plate  with  film  adhesive 
manufactured  by  3~M  Corporation  before  cutting  the  speci¬ 
mens.  This  adhesive  required  1  hour  of  cure  time  at 
250°F.  Because  these  unidirectional  specimens  are  sub¬ 
jected  to  extremely  high  loads,  tabs  were  used  to  prevent 
the  wedge  action  grips  from  digging  into  the  graphite/ 
epoxy  and  to  reduce  the  grip  effect.  None  of  the  other 
specimens  had  end  tabs.  However  the  ends  of  the  specimens 
(approximately  2  in.)  were  coated  with  a  layer  of  MS  907 
two-part  epoxy  manufactured  by  M i 1 1 er-S t ephenson  Chemical 
Company.  This  epoxy  remains  fairly  ductile  after  cure, 
resists  chipping,  and  bonds  well  with  the  graphite/epoxy. 
The  epoxy  coating  was  lightly  sanded  to  a  uniform  thick¬ 
ness.  320  grit  emery  cloth  was  placed  between  the  grips 
and  the  specimen  during  loading  to  improve  the  gripping 
friction.  The  epoxy  layer  and  emery  cloth  serve  to  pro- 


tect  the  specimen  from  the  wear  due  to  the  coarse  diamond 


serration  of  the  grips.  The  two  specimen  types  are  shown 
schematically  in  Figu;  e  2. 

Displacement  measurements  were  made  for  all  tests  by 
means  of  an  extensometer .  The  extensometer ,  used  to 
measure  longitudinal  displacement,  had  a  1  inch  gage 
length  and  was  attached  to  the  center  portion  of  the 
specimen.  The  knife  edges  of  the  extensometer  sat  in 
narrow,  V-shaped  channels  machined  into  aluminum  tabs 
(Figure  3).  The  tabs  were  bonded  to  the  specimen  using 
RTV  silicone  adhesive  manufactured  by  General  Electric 
Corp.  The  extensometer  mounting  arrangement  is  3hown  in 
Figure  4. 

To  prepare  the  specimens  for  edge  replication  (an  NDE 
technique  to  be  described)  the  edges  of  all  the  specimens 
were  polished.  Each  specimen  was  polished  first  with  5.0 
micron  aluminum  oxide  slurry  on  a  felt  polishing  cloth. 
Approximately  400  strokes  on  each  edge  were  used.  Then 
each  was  polished  with  1  micron  slurry  in  a  similar 
fashion.  Approximately  200  strokes  were  needed.  To  check 
polishing  quality,  the  specimens  were  periodically  exam¬ 
ined  under  a  lOOx  power  microscope. 

During  preparation  and  prior  to  being  tested,  the 
specimens  were  stored  at  room  temperature  in  an  airtight 
chamber  containing  desicant  in  order  to  minimize  environ- 
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Edge  Replication 

Surface  replication  is  a  technique  used  to  duplicate 
the  surface  topography  of  a  material  onto  another  medium. 
The  technique  was  orignially  developed  and  used  in  metal¬ 
lography,  and  was  extended  to  use  with  composites  by 
Stalnaker  and  Stinchcomb  [30].  Replicas  are  produced  by 
pressing  cellulose  acetate  tape  softened  by  exposure  to 
acetone  to  the  edge  of  the  composite  specimen.  After  the 
acetone  evaporates  and  the  acetate  hardens,  a  permanent 
record  of  the  surface  features  is  obtained.  This  method 
provides  details  of  the  damage  state  at  the  edge  of  the 
laminate  and,  when  taken  at  various  times  during  loading, 
records  the  history  of  damage  development.  In  paricular, 
details  such  as  number,  spacing,  shape,  and  location  of 
transverse  cracks,  and  delaminations  that  manifest  them¬ 
selves  at  the  edge  are  readily  observable.  An  example  is 
shown  in  Figure  5. 

The  quality  of  the  details  that  could  be  seen  on  an 
edge  replica  was  greatly  increased  when  the  edge  was  pol¬ 
ished  as  described  previously.  Details  were  further  im¬ 
proved  by  making  the  replica  while  the  specimen  was  still 
under  load.  This  served  to  open  the  cracks  to  allow 
better  penetration  of  the  softened  acetate.  After  a  load 
step,  the  specimen  was  reloaded  to  a  modest  fraction  of 
the  previous  load  and  a  replica  was  made.  To  produce  a 


replica,  strips  of  acetate  were  taped  to  an  appropriately 
sized  long  slender  block  with  a  smooth  surface.  The  block 
was  held  lengthwise  against  the  specimen  in  the  load 
frame,  tilted  away  from  the  specimen  slightly,  and  a  small 
stream  of  acetone  was  allowed  to  run  the  length  of  the 
tape.  This  was  quickly  repeated  as  the  block  tilted 
lengthwise  in  the  other  direction.  This  technique  facili¬ 
tated  an  even  coverage  of  acetate  tape  by  acetone.  The 
block  was  held  firmly  in  place  for  approximately  one 
minute.  When  the  tape  was  removed,  its  quality  was 
checked  in  a  25x  power  microfiche  reader  before  proceeding 
with  the  next  load  step.  With  experience,  it  is  possible 
to  place  two  replicas  side  by  side  on  one  strip  of  tape, 
one  replica  of  each  side  of  the  specimen.  When  viewed  in 
the  microfiche  reader,  one  can  easily  ascertain  the  con¬ 
tinuity  of  a  transverse  crack  from  one  edge  of  the  speci¬ 
men  to  the  other. 

To  prevent  the  excess  acetone  from  running  down  the 
specimen  and  softening  the  epoxy  tab  coatings,  an  "apron" 
of  masking  tape  was  placed  at  the  botton  of  the  specimen 
before  making  an  edge  replica.  The  cellulose  acetate  used 
was  No.  11340  Replicating  Tape  manufactured  by  Ernest  F. 
Ful lam ,  Inc . 

X-ray  Radiography 

X-ray  Radiography  is  becoming  a  common  technique  in 
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the  investigation  of  damage  in  advanced  composites.  This 
method  of  non-destructive  evaluation  provides  an  image  of 
the  internal  structure  of  a  material.  The  image  is  pro¬ 
duced  on  film  as  the  material  under  examination  selective¬ 
ly  absorbs  or  transmits  the  radiation.  A  variation  in 
radiation  absorbtion  in  a  composite  material  can  be  caused 
by  the  presence  of  voids,  cracks,  delaminations,  and  in¬ 
herent  flaws.  It  is  obvious  that  this  can  be  an  important 
tool  in  damage  characterization.- 

A  problem  arises  in  that  the  components  of  a 
graphite/  epoxy  material  system  are  organic,  and  thus 
their  absorptivity  of  radiation  is  only  slightly  different 
from  the  surrounding  air.  A3  a  consequence,  an  X-ray 
opaque  substance,  or  enhancing  agent,  must  be  applied  to 
the  composite  to  distinguish  a  damaged  region  from  an 
undamaged  region.  When  the  enhancing  agent  is  introduced 
into  damaged  regions,  the  resulting  X-ray  image,  or  film 
exposure,  shows  these  regions  distinctly. 

A  candidate  for  an  enhancing  agent  must  contain  some 
type  of  "heavy"  element  such  as  iodine,  zinc,  lead,  or 
barium.  Several  enhancing  agents  are  used  both  in  re¬ 
search  and  in  Industry;  tetrabromethane  (TBE),  diiodo- 
butane  (DIB),  di iodomethane ,  and  zinc  iodide  are  examples. 
All  of  these  are  considered  toxic  to  humans  except  zinc 
iodide.  Rummel  and  coworkers  [31]  found  that  a  zinc 
iodide  solution  was  comparable  to  the  other  enhancing 
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agents  in  both  penetration  into  a  composite  material  and 
its  opacity  to  X-rays.  They  also  found  no  influence  of 
the  zinc  iodide  on  the  mechanical  response  of  a  graphite/ 
epoxy  material  system. 

Several  Investigators  have  suggested  a  working  solu¬ 
tion  consisting  of  60  grams  of  zinc  iodide,  10  ml  of 
water,  10  ml  of  isopropyl  alcohol,  and  10  ml  of  Kodak 
"Photo-Flo  600."  This  solution  was  used  for  all  radio¬ 
graphs  in  this  investigation.  The  zinc  iodide  was  ob¬ 
tained  from  Fisher  Scientific  Company.  The  Photo-Flo  acts 
as  a  wetting  agent,  reducing  surface  tension  of  the  solu¬ 
tion  when  applied  to  a  crack.  It  can  be  obtained  from 
most  photographic  supply  outlets. 

In  X-ray  radiography,  it  is  important  to  achieve  as 
high  a  resolution  as  possible  to  increase  the  detail  of 
the  damage  seen  in  a  radiograph.  One  cause  of  decreased 
resolution  is  the  penumbra  effect.  In  the  case  of  visible 
light,  the  penumbra  is  the  area  of  partial  illumination 
between  a  perfect  shadow  and  full  light.  Figure  6  shows  a 
schematic  of  the  arrangement  used  to  make  a  radiograph; 
where  a  denotes  the  size  of  the  anode,  or  source  of 

c;r\cc|.j 

X-rays,  1  is  the  distance  from  the  ebj  oe-t  to  the  ,  d 

is  the  distance  from  the  object  to  the  film,  c  is  the  arc 
length  between  two  intersecting  X-ray  beams,  and  9  is  the 
angle  between  those  beams.  To  maximize  the  resolution  one 
must  minimize  the  arc  length  c.  This  can  be  accomplished 
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by  simply  moving  the  film  as  close  as  possible  to  the 
object,  i.e.  decrease  d.  This  does  not  eliminate  the 
problem  altogether  because  the  object  still  has  a  finite 
depth.  Another  method  is  to  increase  the  source- to-obj ect 
distance.  This  decreases  the  angle  0.  A  third  method  to 
reduce  the  penumbra  effect  is  to  use  a  small  source, 
thereby  decreasing  both  a  and  9. 

There  are  other  nongeometric  means  to  increase  reso¬ 
lution.  In  an  X-ray  exposure,  one  can  only  discern  de¬ 
tails  as  small  as  the  grain  on  the  film.  By  using  a 
fine-grained  film,  more  detailed  features  can  be  seen. 
Another  means  to  increase  resolution  is  to  use  a  film  with 
high  emulsion  contrast.  Scattered  radiation  which  may  fog 
an  exposure  can  be  reduced  by  using  a  collimator  or  dia¬ 
phragm  to  limit  the  exposed  radiographic  field. 

The  voltage  which  operates  the  X-ray  tube  head,  mea¬ 
sured  in  kilovolts  (kVp),  determines  the  intensity  or 
"speed"  of  the  X-ray  radiation.  High  kVp  creates  "hard" 
or  intense  X-rays  while  low  kVp  produces  "soft"  X-rays. 

If  an  exposure  is  made  at  a  kVp  that  is  too  high,  the  rays 
pass  through  the  object  unhindered.  The  film,  when  devel¬ 
oped,  will  be  totally  exposed  as  if  no  object  was  in 
the  path  of  the  X-rays.  If  the  operating  kVp  is  too  low, 
no  X-rays  penetrate  the  object,  and  therefore  no  internal 
damage  will  be  seen  at  all.  The  proper  kVp  setting  de¬ 
pends  on  the  absorptivity  of  the  object  material.  When 
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a  proper  voltage  is  found,  the  developed  film  should  be 
black  in  the  area  surrounding  the  object,  gray  in  the  area 
covered  by  the  object,  and  white  in  the  damaged  areas  that 

I 

have  enhancing  agent  present.  In  this  way  the  voltage 
controls  the  contrast  of  the  radiograph.  Another  factor 
that  controls  the  quality  of  the  image  is  the  time  of  the 
exposure.  This  determines  the  "number"  of  X-ray3  that 
pass  through  the  object.  The  correct  time  of  exposure  is 
dependent  on  the  operating  voltage  and  the  source  distance 
to  the  object.  The  time  of  exposure  is  related  to  the 
square  of  the  distance.  It  the  distance  to  the  object  is 
doubled,  an  identical  exposure  can  be  made  by  increasing 
the  time  of  exposure  by  a  factor  of  four. 

The  X-ray  unit  used  was  a  Phillips  Model  K  1 40  Be 
with  a  140  kVp  range  and  Beryllium  window.  The  nominal 
f ilm-to-source  distance  was  37  inches.  A  lead  bronze 
diaphragm  with  a  radiographic  field  of  4  inches  by  19 
inches  was  used.  The  tube  current  was  2  mA  with  a  .6  mm 
by  .6  mm  focal  spot. 

The  zinc  iodide  solution  was  introduced  to  the  dam¬ 
aged  areas  through  the  edge  of  the  specimen.  The  infil¬ 
tration  of  the  zinc  iodide  can  be  accelerated  by  putting 
the  specimen  under  a  tensile  load.  Loading  the  specimen 
reopens  the  matrix  cracks  and  provides  a  capilary  path  for 
the  solution  to  travel  to  other  damage  forms.  This  load¬ 
ing  was  accomplished  with  a  small  load  frame  i ncor por a t i ng 


Instron  grips  rated  for  5000  pounds  (Figure  7).  Load  was 
applied  by  turning  a  nut  which  drew  a  threaded  shaft 
conrected  to  a  grip  through  a  hollow  shaft.  The  specimen 
was  placed  in  the  load  frame  with  the  edge  upward  and 
loaded  to  approximately  300  pounds.  Two  strips  of  sili¬ 
cone  rubber  were  clamped  to  each  face  of  the  specimen 
approximately  1/4  inch  beyond  the  edge.  Silicone  grease 
compound  from  General  Electric  was  placed  between  the 
strips  at  each  end  of  the  specimen  to  form  a  trough.  The 
zinc  iodide  solution  was  injected  into  this  trough  and 
allowed  to  remain  for  15  minutes  as  gravity  assisted  the 
infiltration  (Figure  8).  This  procedure  was  repeated  for 
the  other  edge  of  the  specimen.  The  residual  solution  was 
then  removed  and  the  specimen  edges  and  surfaces  were 
cleaned  with  methyl  ethyl  ketone. 

All  X-ray  exposures  were  made  while  the  specimens 
were  still  under  load  in  the  small  load  frame.  The  tube 
head  was  mounted  on  an  optical  stand  with  three  axis 
control  for  exact  positioning  of  the  tube  head  (Figure  9). 
The  frame  and  X-ray  tube  head  were  completely  enclosed  in 
a  lead  shield  during  operation  (Figure  10).  The  optimum 
operating  voltage  for  the  laminates  used  was  found  to  be 
30  kVp.  The  film  used  was  Kodak  Industrex  M-type,  double 
emulsion,  high  resolution  film.  Exposure  times  were  30 
seconds  for  4  ply  laminates  and  40  seconds  for  10  ply 
laminates.  An  example  of  the  details  visible  in  an  X-ray 


radiograph  are  shown  in  Figure  11. 

Mechanical  Testing 

Mechanical  testing  was  conducted  in  the  Materials 
Characterization  Laboratory,  Department  of  Aerospace  Engi¬ 
neering,  Texas  A&M  University.  Figure  12  shows  the  array 
of  equipment  used  in  this  investigation.  All  tests  were 
performed  on  a  mechanical  screw-driven  Instron  1125 
Universal  Testing  System  with  wedge  action  grips.  Tests 
were  run  at  a  constant  cross-head  speed  of  .05  inches  per 
minute  as  recommended  in  ASTM  Standard  3039-76.  Extensom- 
eters  were  used  to  measure  both  longitudinal  and  trans¬ 
verse  displacement.  These  extensometers  were  manufactured 
by  Material  Testing  System,  Inc.  (MTS).  An  MTS  632.11B-20 
extensometer  with  a  one  inch  gage  length  was  used  to 
measure  the  longitudinal  displacement.  This  extensometer 
was  seated  on  the  specimen  in  the  aluminum  grooved  tabs 
described  previously.  It  was  held  in  place  by  elastic 
bands.  An  MTS  632.11C-20  extensometer  (metric)  was  used 
to  measure  the  transverse  displacement.  The  flat  sides  of 
the  knife  blades  were  placed  against  the  edges  of  the 
specimen.  An  elastic  band  placed  around  the  "legs”  of  the 
extensometer  kept  the  blades  in  compression  against  the 
specimen  edges.  Another  elastic  band  was  used  to  position 
the  extensometer  on  the  the  specimen.  Cardboard  spacers 
were  used  to  increase  the  gage  length  of  the  transverse 


extensometer  so  that  the  zero  point  or  center  of  the 
response  range  of  the  device  would  be  used  when  mounted. 
The  extensometer s  were  calibrated  just  prior  to  mounting 
for  each  load  step  (Figure  13).  The  extensometers  were 
mounted  on  opposite  sides  at  the  center  portion  of  the 
specimen  (Figure  14). 

The  monitoring  of  stiffness  degradation  was  a  major 
point  of  interest  in  this  investigation.  For  this  reason, 
it  was  imperative  to  insure  the  accuracy  of  the  stiffness 
measurement.  Load  measurement  is  easily  accomplished  with 
a  load  cell  and  is  subject  to  few  error  possibilities.  The 
most  common  methods  of  measuring  strain  are  through  the 
use  of  resistance  strain  gages  and  extensometers.  There 
are  advantages  in  both  techniques.  Strain  gages  are 
permanently  fixed  to  the  specimen  surface  while  extensome¬ 
ters  are  not.  Depending  on  the  mounting  technique,  exten- 
sorneters  may  slip  during  loading.  Strain  gages  are 
sensitive  only  to  the  material  beneath  the  gage.  In  a 
metal,  the  deformation  field  is  continuous  even  down  to 
the  microscopic  level.  A  statistically  representative 
volume  of  the  material  is  several  orders  of  magnitude 
smaller  than  the  size  of  the  strain  gage.  Of  course,  this 
kind  of  deformation  would  be  best  measured  by  resistance 
strain  gages.  In  a  composite  laminate,  the  size  of  the  a 
single  local  damage  phenomenon,  such  as  a  transverse 
crack,  can  be  on  the  same  order  of  magnitude  as  the  size 
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of  a  strain  gage.  This  could  lead  to  a  misinterpretation 
of  the  total  response  of  the  material  system.  Extensome- 
ters,  on  the  other  hand,  integrate  the  material  response 
over  a  comparatively  large  gage  length.  Extensometer 
measurements  are  more  indicative  of  the  total  specimen 
response.  For  this  reason  extensometers  were  used  exclu¬ 
sively. 

The  response  of  the  extensometers  was  amplified  by  a 
signal  conditioner  manufactured  in-house.  The  load  cell 
response  was  processed  by  the  conditioner  built  into  the 
Instron  controller.  The  longitudinal  displacement  and  the 
load  signals  were  monitored  continuously  and  recorded  on 
an  analog  X-Y  plotter. 

A  digital  data  acquisition  system  was  used  to  record 
four  channels  of  data  for  each  test.  The  signals  that 
were  recorded  were  the  load,  longitudinal  displacement, 
transverse  displacement,  and  time.  The  computer  was  a  DEC 
PDP  11/23  Plus  manufactured  by  Digital  Equipment  Corpora¬ 
tion.  The  A-to-D  conversion  board  was  manufactured  by  the 
AD AC  Corporation.  The  test  parameters  entered  into  the 
acquisition  program  were  the  calibration  factors  for  the 
load  cell  and  extensometers,  the  specimen  dimensions, 
cross-head  rate,  and  the  rate  at  which  data  points  were 
recorded.  After  data  reduction,  any  two  of  the  following 
test  outputs  could  be  plotted  against  one  another:  stress, 
longitudinal  strain,  transverse  strain,  Poisson's  ratio, 
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and  time.  Since  the  tests  were  run  at  constant  cross-head 
speed,  this  capability  was  needed  to  determine  stress  and 
strain  rates  for  each  test.  The  digitized  information  was 
plotted  on  a  Hewlett-Packard  7470A  plotter. 

The  Instron  wedge  action  grips  had  the  capability  to 
test  two  inch  wide  specimens.  These  were  used  instead  of 
one  inch  wide  grips  to  prevent  a  loss  of  gripping  effect 
at  the  edges  of  the  one  inch  wide  specimens.  Aluminum 
templates  were  used  to  keep  the  specimen  centered  in  the 
grips.  Once  a  specimen  was  placed  in  the  lower  fixed  grip 
with  a  template,  it  was  aligned  with  the  grip  insert  face 
with  a  level  and  the  grip  was  tightened.  A  steel  straight 
edge  was  used  to  align  the  upper  grip  insert  face  with  its 
lower  counterpart. 

Three  types  of  tests  were  conducted  for  each  laminate 
group.  In  the  first  type  of  test,  two  specimens  of  each 
laminate  were  tested  by  monotonically  loading  to  failure. 
The  specimen  was  fully  instrumented  with  longitudinal  and 
transverse  extensometer s .  The  first  audible  acoustic 
emission  was  recorded  and  any  significant  patterns  in  the 
emissions  were  noted.  This  test  supplied  information  on 
first  ply  failure,  ultimate  strength,  ultimate  strains, 
and  modulus  for  an  undamaged  laminate.  This  information 
set  the  framework  for  the  following  tests. 

In  the  second  type  of  test,  a  single  specimen  from 
each  laminate  was  repeatedly  loaded  and  unloaded,  each 
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succeeding  load  step  greater  than  the  previous  one.  Only 
load  and  longitudinal  strain  were  measured.  An  edge 
replica  was  taken  over  the  center  4  inches  of  the  specimen 
after  each  load  step  when  the  extensometer  was  removed. 

At  no  time  was  the  specimen  removed  from  the  grips.  The 
number  of  load  steps  performed  were  typically  between  8 
and  12  3teps.  This  test  indicated  the  crack  patterns  and 
modulus  changes  to  expect  in  futher  testing. 

The  last  set  of  tests  involved  a  procedure  similar  to 
the  one  described  above,  except  that  an  X-ray  radiograph 
was  made  upon  unloading  from  each  step.  Of  course  this 
necessitated  the  removal  of  the  specimen  from  the  grips 
also.  In  most  cases,  the  longitudinal  and  transverse 
strains  were  monitored.  The  load  level  to  which  a  speci¬ 
men  was  taken  was  based  on  the  amount  of  audible  damage 
that  could  be  heard  when  loading  a  specimen.  The  object 
of  this  was  to  develop  a  consistent  increase  in  damage 
(i.e.  transverse  cracks)  for  each  load  step. 

Care  was  taken  to  avoid  holding  a  specimen  at  extreme 
loads  where  cracking  occurs  spontaneously  without  the 
addition  of  load.  Attempting  to  make  an  edge  replica  at 
this  peak  load  would  further  complicate  this  situation. 
Once  a  test  had  reached  a  desired  load  level,  it  was 
immediately  returned  to  zero  load  at  the  same  cross-head 
rate  as  loading.  The  extensometer  was  then  removed  and 
the  specimen  was  reloaded,  typically  to  300  pounds.  Edge 


replicas  were  made  of  the  full  seven  inches  between  the 
grips,  on  both  edges  of  the  specimen. 

Determination  of  Stiffness 

As  mentioned  before,  an  important  aspect  to  this 
investigation  was  to  monitor  stiffness  degradation.  For 
the  step-wise  loading  test,  stiffness  was  measured  only  on 
the  unloading  part  of  the  test.  In  this  way,  the  stiff¬ 
ness  reflected  the  effects  of  damage  accumulation  during 
that  particular  load  step. 

A  least-squares  linear  curve  fit  was  used  to  deter¬ 
mine  the  stiffness.  A  computer  program  was  created  to  do 
the  numerical  fit  to  the  digitized  data.  The  fit  was  made 
between  specified  high  and  low  load  values. 

All  laminate  types  tended  to  have  a  stress-strain 
curve  that  was  concave  up.  As  a  consequence,  a  stiffness 
calculated  from  the  1 eas t- s quar es  fit  was  dependent  on  the 
specified  load  levels.  A  systematic  method  of  deriving 
the  stiffness  was  needed. 

From  preliminary  tests,  it  was  found  that  a  typical 
one  inch  wide  0°  lamina  could  carry  an  ultimate  load  of 
appr oxima tel y  1500  pounds.  Because  all  the  laminates  in 
this  investigation  were  fiber  dominated,  this  information 
was  useful  in  determining  the  ultimate  load  of  each  of  the 
laminates.  For  example,  the  failure  load  of  a  [02/902] 

s 

laminate  was  approximately  6000  pounds.  It  was  found  that 


the  first  indication  of  damage  occurred  at  about  half  of 
this  estimated  ultimate  load  for  all  the  laminates.  This 
made  a  convenient  high  load  for  the  least-squares  calcula¬ 
tion  and  was  used  in  all  cases  for  the  determination  of 
stiffness . 

Since  the  tests  were  conducted  at  a  fixed  cross-head 
rate,  the  stress  or  strain  rates  were  not  necessarily  con¬ 
stant.  Figure  15  shows  an  example  of  this.  When  a  test 
was  initiated  or  neared  complete  unloading,  the  stress  and 
strain  rates  were  affected  as  the  pressure  of  the  grip 
wedges  eased.  The  low  load  for  the  least-squares  fit  was 
chosen  to  avoid  this  condition. 

Measurement  of  Hesldual  Strain 

One  of  the  goals  of  this  investigation  was  to  monitor 
residual  strain  after  each  load  step.  This  proved  to  be 
anomalous.  When  a  specimen  was  pulled  to  the  desired 
load,  it  was  immediately  unloaded.  At  this  point  the 
stress-strain  curve  would  sometimes  "loop"  back  on  itself. 
That  is,  the  strain  signal  would  decrease  sharply,  while 
the  load  signal  decreased  only  slightly.  This  is  illus¬ 
trated  in  Figure  16.  This  phenomenon  was  not  consistent 
from  test  to  test.  Since  the  residual  strain  upon  unload¬ 
ing  was  very  small,  it  was  of  the  same  order  of  magnitude 
as  this  loop.  An  inertial  effect,  backlash  in  the  drive 
train,  or  movement  of  the  extensometer  may  be  explanations 
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Loop  formed  in  stress-strain 
during  load  reversal 


curve 


for  this  occurrence. 


Measurement  of  Surface  Area 

Another  objective  of  this  investigation  was  ;o  quan¬ 
titatively  determine  the  surface  area  produced  during 
damage  accumulation.  Previous  research  concentrated  on 
determining  the  number  of  cracks  produced  during  loading 
and  the  resulting  crack  density  (number  of  cracks  per  unit 
length).  This  quantity  would  then  be  compared  to  the 
stiffness  loss  of  the  laminate.  A  simple  modification  was 
made  to  this  procedure  to  obtain  surface  area.  Thi3  was 
done  by  taking  the  product  of  the  number  of  cracks,  the 
width  of  the  specimen,  the  thickness  of  the  90°  group,  and 
a  factor  of  2  (for  each  face  of  a  crack). 

Preliminary  tests  indicated  that  many  transverse 
cracks  did  not  grow  self similarly .  From  the  inspection  of 
edge  replicas,  it  was  found  that  many  cracks  had  irregular 
shapes  and  curvatures,  and  in  some  incidences  the  cracks 
did  not  propogate  through  the  full  thickness  of  the  trans¬ 
verse  ply  groups.  X-ray  radiographs  indicated  that  some 
cracks  did  not  span  the  full  width  of  the  specimen.  It 
would  seem  that  an  incomplete  crack  would  not  have  the 
same  impact  on  the  response  of  the  laminate  as  would  a 
full  crack.  Simply  counting  it  as  just  another  crack  and 
using  the  aforementioned  procedure  could  be  misleading. 

To  avoid  this  mi srepr esentat i on ,  crack  surface  area 


5 


measurements  were  made  on  a  crack  by  crack  basis.  In  this 
way  compensation  could  be  made  for  incomplete  cracks  and 
crack  curvature.  The  basis  for  the  area  measurement  w? s 
to  determine  crack  length  through  the  thickness  of  the 
laminate.  Crack  lengths  were  measured  from  edge  replicas 
with  a  lOOx  power  microscope.  The  eyepiece  of  the  micro¬ 
scope  incorporated  a  calibrated  100  division  scale.  At 
the  magnification  used,  there  were  100  divisions  per 
millime  ter,  or  approximately  25  divisions  per  lamina 
taickness . 

Crack  lengths  were  measured  on  both  sides  of  the 
specimen  and  then  averaged.  In  this  way  changes  in  the 
length  (through  the  thickness)  across  the  width  of  the 
specimen  could  be  accounted  for.  By  doubling  the  length 
and  multiplying  by  the  width  of  the  specimen,  the  crack 
surface  area  was  approximated.  Since  a  nondestructive 
technique  was  not  available  to  measure  the  surface  area 
directly,  it  was  assumed  that  this  would  be  a  suitable 
approximation. 

The  angle  at  which  cracks  grew  across  the  thickness 
of  the  specimen  wa3  also  of  interest.  A  video  camera  was 
attached  to  the  microscope  and  the  picture  of  the  edge 
replica  was  displayed  on  a  large  video  moniter.  Angles 
could  be  measured  accurately  using  a  large  homemade  pro- 


51 


$ 


£ 


a 


« 


'"J* 


1 


9 

i 


£ 

fi 


RESULTS  AND  DISCUSSION 

The  results  are  presented  in  three  parts.  First,  a 
general  overview  of  the  test  results  is  presented.  Sec¬ 
ond,  observations  are  presented  for  each  laminate  type. 
Third,  important  comparisons  are  made  for  similar  lami¬ 
nates. 

Overview  of  Results 

As  mentioned  before,  the  tests  in  this  investigation 
were  designed  to  promote  the  growth  of  matrix  cracks.  In 
addition  to  matrix  cracking,  some  other  forms  of  damage 
did  occur.  These  included  delaminations,  fiber  breakage, 
fiber  splitting,  and  transverse  crack  branching,  depending 
on  the  laminate  type. 

A  coordinate  system  shown  in  Figure  17  has  been  used 
to  help  describe  the  various  damage  phenomenon.  The  x- 
and  y-axes  define  the  lamina  plane  with  the  x-axis  aligned 
with  the  longitudinal  or  loading  direction.  The  z-axis  is 
perpendicular  to  the  lamina  plane  and  goes  from  front 
surface  to  back  surface. 

A  compilation  of  the  ultimate  strengths  and  ultimate 
strains  of  all  the  laminates  tested  is  shown  in  Table  4. 
The  data  for  un i d i r ec t i onal  0°  and  90°  laminates  are  in¬ 
cluded.  All  laminates  had  an  ultimate  strain  less  than 
that  of  the  unidirectional  0°  laminate  due  to  the  complex 
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Table  4.  Laminate  Ultimate  Properties 


TYPE 

LAY-UP 

ULTIMATE 

STRENGHTH 

( k3i  ) 

ULTIMATE 

STRAIN 

* 

Z 

[o.]g 

326.0 

± 

1  1  . 

1  .  44 

±  .07 

N 

C90.]3 

11.1 

± 

1  . 1 

.77 

±  .05 

A 

[0/90„]a 

51  .2 

± 

3.0 

1  .22 

±  .12 

B 

[0/90  $  ]g 

74.3 

± 

1  . 1 

1  .33 

±  .09 

C 

[0/90  2  ]g 

106.9 

± 

1  .  4 

1  .41 

±  .03 

D 

[ J/90]a 

150.3 

± 

6.5 

1  .  32 

±  .07 

E 

[02/902]a 

159.9 

t 

6 . 0 

1  .  38 

±  .05 

F 

[0/90/0] 

3 

181.7 

± 

9.8 

1  .  32 

±  .08 

G 

[0/90/0/90] 

s 

165.9 

± 

8.6 

1  .  40 

±.11 
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localized  stress  states  associated  with  various  damage 
modes  present  in  the  transverse  layers.  The  [0/90„] 

3  , 

(Type  A)  had  the  lowest  ultimate  strain  due  to  extensive 
damage  formation.  The  transverse  crack  development  was 
irregularly  spaced  and  delamination  was  widespread.  No 
consistent  pattern  could  be  found  for  the  remaining  lami¬ 
nates  and  further  comparison  of  the  laminates  with  respect 
to  ultimate  strain  was  inconclusive. 

During  each  test  of  an  undamaged  specimen,  the  load 
at  the  first  audible  acoustic  emission  was  recorded. 

Table  5  lists  this  information  for  each  laminate  type.  j 

For  comparison,  the  ultimate  strength  and  strain  data 

for  a  unidirectional  90°  laminate  is  included  in  the 

table.  The  first  "pop"  was  assumed  to  be  the  formation  of 

the  first  transverse  crack.  This  was  termed  as  the  first  j 

ply  failure  (FPF).  This,  of  course,  is  a  misnomer  since 

the  ply  is  able  to  carry  a  load  along  the  remainder  of  its 

length.  In  general,  the  strain  at  FPF  was  inversely 

related  to  the  number  of  90°  plies  in  a  group.  This  is 

easily  seen  in  laminates  types  A  through  D.  The  [ 0 / 9  0  2  ]  ■ 

s 

and  [02/902]  laminates  both  had  four  90°  layers  and  had 

9 

very  nearly  the  same  FPF  strain.  With  the  large  group  of 

i 

' 

90°  plies,  the  additional  constraint  of  the  0°  plies  in 
the  [02/9Q2]  laminate  had  little  effect.  The  constrain- 
ing  effect  is  evident  in  the  [0/90/0]  laminate.  The  FPF 

3 

strain  is  almost  the  same  as  the  ultimate  strain  of  a 
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Table  5.  Stress  and  Strain  at  First  Ply  Failure  ( FPF ) 


TYPE 


LAY-UP 


FPF 

STRENGHTH 


s 


FPF 

STRAIN 


C90s]a 

ii.i 

± 

1  . 1 

.77 

±  .05 

[0/90„]a 

14.5 

± 

3.8 

•  30 

±  .08 

[0/90j]3 

19.0 

± 

1  . 1 

.31 

±  .01 

[ 0/90  2 ] g 

44.1 

± 

8.0 

.55 

±  .09 

[0/90] 

72.9 

± 

7.7 

.63 

±  .09 

[02/902]g 

64.4 

± 

6.5 

.56 

±  .06 

[0/90/0] 

S 

101.4 

± 

5.6 

.74 

±  .03 

[0/90/0/90] 

77.9 

± 

2 . 2 

.69 

±  .02 

<!• 


unidirectional  90°  laminate. 

As  described  previously,  a  matrix  crack  is  a  crack 
that  runs  parallel  to  the  fibers  in  any  ply.  A  transverse 
crack  is  a  matrix  crack  in  a  ply  whose  fibers  are  not 
oriented  in  a  principal  loading  direction.  This  was  the 
predominant  mode  of  damage,  as  should  be  expected  for 
cross-ply  laminates.  Of  significant  importance  was  the 
observation  of  two  different  types  of  transverse  cracks. 
These  will  be  referred  to  as  "straight"  cracks  and  "angle" 
cracks.  A  description  follows. 

Figure  18  shows  a  typical  transverse  crack  in  a  group 
of  90°  plies.  As  a  rule,  this  type  of  crack  was  contin¬ 
uous  across  the  width  of  a  specimen  and  grew  in  a  self¬ 
similar  manner  through  the  thickness  of  the  90°  plies.  As 
seen  in  an  edge  replica,  the  crack  generally  passed  around 
fibers  in  its  path  and  extended  into  the  resin  rich  region 
at  the  0°/90°  ply  interfaces.  The  crack  extended  from 
interface  to  interface,  nearly  perpendicular  to  the  plane 
of  the  laminae.  It  rarely  deviated  from  this  course  more 
than  3  or  4  fiber  diameters.  When  the  front  surface  of 
the  specimen  was  viewed  in  an  X-ray  radiograph,  the 
transverse  crack  appeared  as  a  sharp,  well  defined  line 
running  across  the  specimen  width  from  edge  to  edge.  A 
transverse  crack  of  this  type  wa:.  labeled  as  a  "straight 
crack."  This  type  of  crack  was  found  in  all  cross-ply 
laminates  and  was  the  exclusive  type  found  in  laminates 
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that  had  groups  of  90°  plies  of  less  that  two. 

The  second  kind  of  transverse  crack  was  called  an 
"angle"  crack.  Several  features  distinguished  this  type 
of  crack  from  a  straight  crack.  The  most  obvious  of  these 
is  its  shape  as  seen  on  an  edge  replica.  As  the  name 
implies,  the  crack  propagates  at  an  angle  to  the  straight 
crack  from  a  0V90  0  interface.  One  form  of  the  angle 
crack  is  a  par  tial  crack  which  does  not  span  the  full 
thickness  of  a  group  of  90°  plies.  A  "curved"  crack  is  an 
extension  of  the  angle  crack  in  which  cracks  from  the  two 
interfaces  meet  at  some  point  in  the  interior  of  the  90° 
group.  An  example  of  a  partial  crack  (a)  and  a  curved 
crack  (b)  are  shown  in  Figure  19. 

When  viewed  closely  in  an  X-ray  radiograph,  an  angle 
crack  can  be  easily  distinguished  from  a  straight  crack. 
When  the  crack  is  tilted  in  relation  to  the  X-ray  beam  it 
presents  a  greater  projected  area  but  is  less  effective  in 
absorbing  the  radiation.  As  a  result,  the  angle  crack 
appears  as  a  fuzzy  band  traversing  the  specimen  while  a 
straight  crack  appears  as  a  sharp  fine  line  (Figure  20). 

Figure  21  establishes  a  coordinate  system  to  describe 
the  angle  crack  phenomenon.  As  seen  from  the  edge,  the 
specimen  is  laying  on  its  front  surface  and  a  straight 
crack  is  parallel  to  the  z-axis  and  is  assigned  an  angle 
of  0°.  Rotation  clockwise  from  the  z-axis  is  positive  and 
counterclockwise  is  negative.  All  crack  angles  are 
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between  -90°  and  +90°.  The  curved  crack  could  grow  from  4 
to  30  fiber  diameters  in  the  x-direction,  depending  on  the 
initial  angle  and  the  thichness  of  the  90°  layer.  A 
single  ply  thickness  has  between  18  and  22  fiber  diameters 
through  the  thickness. 

Another  characteristic  of  the  angle  crack  is  its 
proximity  to  a  straight  crack.  There  was  not  a  single 
incidence  of  a  curved  or  partial  crack  forming  as  an 
isolated  crack.  They  always  formed  relatively  close  to  a 
straight  crack  or  another  angle  crack.  From  observation, 
it  was  found  that  the  following  rule  applied.  The  ratio 
of  the  thickness  of  the  90°  layer  to  the  distance  beween 
two  cracks  was  generally  less  than  1.0  for  a  straight 
and  curved  crack.  However,  the  ratio  was  usually  greater 
than  1.5  between  two  straight  cracks  or  two  groups  of 
cracks . 

As  can  be  seen  in  Figure  19,  these  cracks  also  grow 
toward  the  straight  crack.  In  the  case  of  a  curved  crack, 
the  convex  side  always  faced  toward  the  straight  crack. 

It  was  not  uncommon  for  several  curved  cracks  to  nest 


together;  all  associated  with  one  straight  crack  (Figure 
22).  This  observation  prompts  the  definition  of  a  "crack 
family."  A  crack  family  is  a  group  consisting  of  two  or 
more  cracks,  where  one  and  only  one  can  be  considered  a 
straight  crack.  Figure  23  shows  an  example  of  several  of 
these  groups. 
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In  the  [0/90„]  and  [0/90,]  laminates,  some  delam- 
s  s 

ination  at  the  0°/90°  interface  did  occur.  The  delamina¬ 
tions  propagated  from  the  point  where  an  angle  crack  met 
the  0°  ply  (Figure  24).  A  delamination  was  never  observed 
to  originate  from  a  straight  crack. 

The  process  by  which  the  transverse  crack  density 
(number  of  cracks  per  unit  length)  increases  has  been  ex¬ 
plained  as  a  shear-lag  effect.  That  is,  after  an  initial 
transverse  crack  forms,  the  stress  field  in  the  surround¬ 
ing  area  is  relieved.  Moving  away  from  the  crack  in  the 
load  direction,  the  load  carried  by  0°  plies  is  redistrib¬ 
uted  into  the  90°  layers  through  shear.  The  stress  field 
then  approaches  that  of  an  undamaged  laminate  at  some 
finite  distance  from  the  crack.  At  this  point  another 
crack  may  form.  This  explanation  does  not  account  for  the 
formation  of  angle  cracks  near  straight  cracks.  Obviously 
another  damage  mechanism  is  driving  the  growth  of  these 
cracks . 

By  comparing  two  laminates  with  similar  geometry  but 

different  damage  modes,  some  possible  explanations  can  be 

attempted.  Laminate  theory  predicts  that  a  [0/90]  and 

s 

[0,/90,]  will  have  identical  mechanical  behavior.  One  is 
*  s 

a  scaled  version  of  the  other.  In  reality,  the  larger 
laminate  developed  a  complex  pattern  of  curved,  partial, 
and  straight  cracks,  while  the  thinner  laminate  developed 
a  dense,  evenly  spaced  pattern  of  straight  cracks.  Some 
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Figure  24.  Microphotograph  of  delamination  propagating 
from  an  angle  crack  at  the  0°/90o  interfaca 
of  a  [0/90,  ]_  laminate 
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alteration  of  the  stress  field  must  be  involved  in  the 
scaling  process.  By  examining  the  factors  which  are  not 
scaled,  an  explanation  may  be  found.  Several  features  1o 
not  change  from  one  laminate  to  another.  Fiber  diameter 
and  ply  thickness  remain  the  same.  The  resin  rich  region 
at  a  0°/90°  is  also  nearly  the  same.  Recent  work  [23] 
indicates  that  the  resin  rich  area  is  responsible  for  most 
of  the  load  transfer  from  the  0°  plies  to  the  90°  plies. 
When  the  laminate  is  scaled  up,  the  constant  thickness 
resin  rich  region  must  transmit  a  correspondingly  larger 
load.  This  coupled  with  a  stress  concentration  near  a 
straight  transverse  crack  could  possibly  cause  initiation 
of  an  angled  crack. 

Figure  25  demonstrates  the  relationship  between  the 
number  of  straight  and  curved  cracks  with  increasing 
stress  for  a  [02/902]  laminate.  As  expected,  the  total 

w 

number  of  cracks  increased  with  each  additional  increase 
in  stress.  Initially,  the  only  cracks  present  are 
straight.  Curved  cracks  begin  forming  at  higher  levels  of 
3tress.  Note  that  as  the  ultimate  stress  is  approached, 
the  straight  cracks  practically  cease  to  form  while  the 
curved  cracks  are  actively  increasing.  The  saturation 
of  straight  cracks  may  correspond  to  a  char ac t er  i  s t  i  c 
damage  state.  In  some  specimens  it  appears  that  the 
straight  cracks  have  a  consistent  and  regular  spacing. 

The  belief  that  the  typical  straight  transverse  crack 


Stress  (ksi) 

Figure  25.  Crack  development  by  type  for  [02/902]  laminate 
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propagates  in  a  rapid,  unstable  fashion  is  widely  accept¬ 


ed.  This  is  supported  by  the  fact  that  a  partially  formed 


crack  is  rarely  found.  In  laminates  with  large  groups  of  •  j 


90°  plies,  the  formation  of  a  transverse  crack  in  tne 


extensometer  gage  length  causes  a  clear  discontinuity  in 


the  stress-strain  curve.  This  phenomenon  would  support 


the  idea  of  unstable  crack  growth.  The  clear,  distinct 


pop  that  is  heard  when  the  crack  forms  would  also  indicate 


an  almost  instantaneous  crack  propagation. 


This  is  not  the  case  with  angle  cracks.  At  stress 


levels  where  angle  cracks  form,  the  familiar  popping  sound 


becomes  less  distinct  and  begins  to  sound  more  like  a  tear 


than  a  pop.  This  phenomenon,  along  with  the  observation 


of  partial  cracks,  indicates  that  the  fracture  process  is 


considerably  slower. 


Although  it  is  not  clear  at  what  point  a  straight 


crack  initiates,  this  is  not  the  case  for  an  angle  crack. 


The  evidence  presented  indicates  that  the  angle  crack 


initiates  at  the  0°/90°  interface.  Figures  26  through  29  I 

i 

are  a  sequence  of  edge  replicas  at  increasing  load  levels.  I 


The  curved  crack  is  formed  by  the  merging  of  two  partial 


cracks  growing  from  opposite  ply  interfaces.  Figure  30  I 

< 

i 

shows  the  result  of  two  partial  cracks  "overshooting"  one  < 

< 

another.  A  partial  crack  from  one  interface  has  inter-  < 


cepted  a  partial  crack  growing  from  the  opposite  interface  i 

i 

( 

at  a  point  behind  the  original  crack  tip.  i 
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Figure  28.  Edge  replica  of  [0/903]  laminate  showing 
development  of  a  curve!  crack,  8 1  %  F. 
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Several  other  forms  of  minor  damage  were  observed. 

In  some  instances,  a  transverse  crack  would  fracture  a 
fiber  lengthwise  rather  that  propagate  in  the  matrix 
around  the  fiber.  Figure  31  shows  an  example  of  this. 
Sectioning  along  the  specimen  length  might  reveal  the 
extent  of  fiber  splitting.  Fiber  fracture  could  sometimes 
be  seen  at  the  edge  of  the  specimen  in  the  0°  plies. 

These  breaks  were  typically  near  the  0°/90°  interface  and 
adjacent  to  termination  point  of  straight  transverse 
cracks.  The  fiber  breaks  were  not  generally  associated 
with  angle  cracks.  Figure  32  shows  a  detail  of  a  typical 
fiber  break  zone.  This  phenomenon  was  not  monitored 
closely  and  the  effect  on  stiffness  measurements  is  not 
known.  Oven  deply,  another  destructive  evaluation  tech¬ 
nique,  may  be  used  to  find  the  extent  of  fiber  fracture 
across  the  width  of  the  specimens.  Fiber  breaks  and  fiber 
splits  were  found  in  all  the  laminates  tested. 

In  the  [0/902]  and  [02/902]  laminates,  a  different 
type  of  crack  pattern  was  observed  on  the  edges  of  the 
specimens.  This  phenomenon  was  termed  "crack  branching." 
An  example  is  shown  in  Figure  33.  This  occurred  at  the 
point  where  a  transverse  crack  met  the  0°/90°  interface 
and  was  in  the  laminate  plane.  The  branches  grew  from 
both  angle  cracks  and  straight  cracks.  Branches  were  seen 
at  the  crack  tip  of  every  transverse  crack.  These 
branches  appeared  to  be  delaminations;  although,  a  close 
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inspection  of  the  radiograph  revealed  no  evidence  of 
delaminations.  Unlike  a  delamination,  which  develops  as  a 
separate  form  of  damage,  these  branches  formed  as  an 
integral  part  of  the  transverse  cracks.  The  crack  open¬ 
ing,  as  seen  at  the  edge,  is  less  than  half  that  of  the 
transverse  crack.  Because  nothing  can  be  seen  in  the 
radiograph,  the  distance  that  these  branches  extend  into 
the  width  of  the  specimen  is  unknown.  It  i3  possible  that 
this  form  of  damage  is  restricted  to  the  edge.  Longitudi¬ 
nal  sectioning  of  the  specimen  may  further  explain  these 
cracks . 

Type  A  C0/90^]a 

This  laminate  type  had  the  largest  grouping  of  90° 
plies.  As  a  result,  a  very  complex  damage  pattern  devel¬ 
oped.  Although  the  damage  mode  of  interest  was  matrix 
cracking,  the  chief  form  of  damage  which  developed  was 
delami na t i on .  The  X-ray  radiographs  in  Figure  3^  show  the 
incremental  growth  of  the  damage.  Straight  cracks,  angle 
cracks,  longitudinal  splits,  and  delami na t i ons  can  all  be 
seen  in  these  radiographs. 

Initial  transverse  cracking  was  sporadic.  No  pattern 
in  the  spacing  developed.  In  most  cases,  the  crack  pat¬ 
tern  which  developed  in  the  gage  length  beneath  the  exten- 
someter  was  not  representat  i  ve  of  the  damage  elswhere  in 
the  specimen. 


Figure  34.  X-ray  radiographs  for  a  pr ogr es s i vel y  loaded 

[0/90H]  laminate  (dark  rectangles  are 
aluminum  extensometer  mounting  tabs) 


This  laminate  was  susceptible  to  angle  cracks. 

Straight  cracks  stopped  forming  at  a  relatively  low  stress 
level  and  a  rapid  increase  in  angle  cracks  followed.  The 
cracks  were  counted  at  each  loadstep  and  the  number  of 
cracks  formed  were  plotted  against  their  respective  stress 
levels.  Figure  35  shows  that  the  angle  cracks  far  out¬ 
number  the  straight  cracks,  near  failure.  Large  groups  of 
angle  cracks  developed  throughout  the  specimen. 

Delaminations  grew  readily  in  this  laminate.  A 
delamination  originated  at  the  free  edge  and  propagated 
along  the  length  and  width  in  the  plane  of  the  laminate. 

A  delamination  formed  at  the  edge  soon  after  an  angle 
crack  formed.  The  delaminations  then  propagated  as  narrow 
bands  across  the  width  of  the  specimen.  These  bands  can 
easily  be  seen  in  Figure  34.  A  typical  delamination 
originated  at  the  edge  of  the  specimen  and  at  the  tip  of 
angle  crack.  The  cause  of  the  delaminations  is  not  clear. 
The  possible  causes  could  be  interlaminar  stress  at  the 
edge  or  the  localized  stress  field  at  the  matrix  crack 
tip.  The  direction  of  growth  along  the  length  of  the 
specimen  was  always  away  from  the  angle  crack  group.  An 
example  is  shown  in  Figure  36. 

The  last  radiograph  in  Figure  34  indicates  the  wide¬ 
spread  delaminations  in  this  laminate.  A  radiograph  does 
not  show  the  depth  at  which  the  delaminations  are  located. 
In  fact,  the  delaminations  generally  form  in  an  alternat- 


>y  type  for  [0/90„]  laminate 


ing  pattern  along  the  length  of  the  specimen.  What  ap¬ 


pears  to  be  side  by  side  delaminations  in  the  radiograph 
are  actually  at  opposite  C°/90°  interfaces.  This  pattern 
was  verified  by  examining  the  edge  replica. 

Longitudinal  splits  also  formed  throughout  this  lami¬ 
nate.  These  splits  were  found  only  in  areas  of  delamina¬ 
tions.  No  split  was  seen  crossing  a  transverse  crack  in 
an  adjacent  ply  except  when  a  transverse  crack  formed 
after  the  split.  A  split  appeared  in  the  0°  ply  over  a 
delamination  area  and  did  not  extend  beyond  this  boundary. 
No  split  was  seen  without  an  accompanying  delamination 
zone.  Sufficient  evidence  was  not  found  to  conclude  that 
the  delaminations  initiated  at  the  intersection  of  a 
transverse  crack  and  a  longitudinal  split.  There  is, 
however,  evidence  that  some  longitudinal  splits  formed 
after  the  growth  of  a  delamination  zone.  Figure  31*  shows 
the  formation  of  a  split  over  an  already  existing  delami¬ 
nation.  In  this  case,  the  delaminations  may  drive  the 
growth  of  the  longitudinal  splits. 

Just  prior  to  ultimate  failure,  this  laminate  incur¬ 
red  as  much  as  a  15%  reduction  in  stiffness.  This  reduc¬ 
tion  was  a  result  of  all  the  damage  modes  described  above. 
Figure  37  graphically  represents  the  stiffness  reduction 
as  a  function  of  stress  level.  The  sharp  decrease  in  the 
stiffness  occurred  at  the  stress  at  which  de 1  am i na t i ons 
appeared . 


Since  the  measurement  of  crack  surface  area  was  per¬ 
formed  only  for  matrix  cracks,  a  correlation  between  the 
total  surface  area  (ircluding  delaminations)  and  the 
stiffness  reduction  was  not  made.  Figure  38  does  show 
this  correlation  prior  to  the  formation  of  delaminations 
in  the  gage  length.  The  crack  lengths  and  angles  at  which 
the  transverse  cracks  grew  are  listed  in  Table  6. 

Type  B  [0/90, 

Like  the  previously  described  laminate,  the  [0/903] 

s 

also  developed  some  edge  delaminations,  although  not  on 
such  a  large  scale.  Typically  three  or  four  edge  delami¬ 
nations  appeared  within  the  full  length  of  the  specimen. 
The  series  of  radiographs  in  Figure  39  shows  the  progres¬ 
sion  of  damage  with  increased  load.  As  was  the  case  with 
the  [0/90,,]  laminate,  the  damage  in  the  gage  length  was 
not  necessarily  representative  of  damage  growth  elswhere 
in  the  laminate. 

This  laminate  developed  some  longitudinal  splits  in 
areas  other  than  over  delamination  zones.  These  splits 
were  found  in  groups  of  three  to  five,  along  individual 
transverse  cracks.  These  splits  were  limited  in  size  to 
approximately  .025  inches  long.  Unlike  the  longitudinal 
splits  associated  with  edge  delaminations,  these  splits 
could  be  found  bisecting  single  straight  transverse 
cracks,  rather  than  just  growing  from  one  side  of  a  cur  /ed 


transverse  crack. 


A  few  cases  of  interlaminar  delaminations  were  found 
at  the  intersection  of  longitudinal  splits  and  transverse 
cracks.  An  interlaminar  delamination  is  one  which  devel¬ 
ops  in  the  interior  of  a  laminate  with  no  connection  with 
the  free  edge.  These  were  extremely  small  compared  to  the 
size  of  the  edge  delaminations.  Angle  cracks  constituted 
more  than  half  of  the  total  number  of  cracks  formed  at 
failure.  Figure  40  shows  that  the  straight  cracks  formed 
at  a  steady  rate,  while  curved  cracks  began  forming  very 
rapidly  at  a  higher  stress  level.  Table  7  numerically 
documents  the  crack  pattern  for  comparison  to  stiffness 
reduction. 

Type  C  [0/90, 1 

s 

No  delaminations  or  longitudinal  splits  were  gener¬ 
ated  in  this  laminate.  The  transverse  cracking  pattern 
was  very  similar  to  the  [0/903]  laminate,  although  the 
crack  spacing  was  slightly  more  dense  and  angle  cracks 
were  not  as  predominant.  The  radiographs  in  Figure  41 
3how  the  progression  of  damage  with  each  load  step.  There 
seemed  to  be  a  more  even  distribution  of  transverse  crack¬ 
ing  along  the  length  of  the  specimen.  Figure  42  shows 
that  at  95%  of  the  ultimate  strength,  approximately  half 
of  the  transverse  cracks  are  straight  cracks.  In  this 
laminate  there  appears  to  be  a  more  clearly  defined  point 
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Table  7.  Crack  Angle  and  Length  Distribution  in 

Specimen  B-5,  [0/903]  ,  94J  F 

s  tu 


Length  of  cracks  (L)  in  inches* 


Angle 

1st 

2nd 

3rd 

4th 

5th 

6th 

COS0 

cosG 

0 

ply 

ply 

ply 

ply 

ply 

ply 

x  L 

x  L. 

-45 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

-40 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

-35 

.001  4 

.0 

.0 

.0010 

.0061 

.0061 

.0120 

.0083 

-30 

.0030 

.0 

.0 

.0 

.0029 

.0055 

.0099 

.0057 

-25 

.0 

.0 

.0 

.0028 

.0 

.0 

.0025 

.0012 

-20 

.0 

.0 

.0 

.0034 

.0 

.0 

.0032 

.0012 

-15 

.0 

.0 

.0028 

.0028 

.0028 

.0028 

.0108 

.0029 

-10 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

-  5 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

0 

.0884 

.0915 

.0829 

.0825 

.0886 

.0819 

.5158 

.0 

5 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

10 

.0 

.0 

.0055 

.0 

.0028  . 

.0 

.0082 

.001  4 

15 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

20 

.0 

.0 

.0026 

.0028 

.0028 

.0028 

.0103 

.0038 

25 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

30 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

35 

.0030 

.0 

.0 

.0 

.0 

.0 

.0025 

.0017 

40 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

45 

.0030 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

Total 

.0948 

.0915 

.0938 

.0953 

.1060 

.991 

.5752 

.0262 

Length 

Area 

Generated 

by  Transverse  Cracking 

Area 

-  2  x  L  x 

W 

(W  -  1.005 

in. ) 

Total 

area 

of  all 

transverse 

cracks 

1.1662 

in.2 

Projected  area  perpendicular  to  laminate  plane  .  .  .  1.1562  in.2 

Projected  area  parrallel  with  laminate  plane . 0527  in.2 

‘Each  entry  is  the  total  length  of  all  cracks 
in  the  gage  length  at  the  specified  angle. 
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at  which  straight  cracks  saturate  and  curved  cracks  begin. 
This  laminate  showed  a  transverse  crack  density  of  48 
cr\ck  per  inch  at  92$  of  the  ultimate  stress.  Typically 
very  few  of  these  cracks  would  be  partial  cracks.  Table  8 
was  used  to  determine  the  surface  area  generated  by  trans¬ 
verse  cracks. 

If  the  curved  and  partial  cracks  were  ignored,  the 
straight  cracks  appeared  to  show  a  consistent,  even  spac¬ 
ing.  This  would  correspond  to  the  saturation  of  straight 
cracks  seen  in  Figure  42.  This  might  be  considered  a 
characteristic  damage  state  with  approximately  29  cracks 
per  inch. 

Type  D  [0/90], 

This  laminate  was  the  most  difficult  to  test.  Be¬ 
cause  it  had  only  four  plies,  it  was  susceptible  to  vari¬ 
ous  gripping  and  alignment  problems.  The  total  stiffness 
loss  at  ultimate  failure  was  only  three  percent. 

This  laminate  developed  a  very  consistent  transverse 
crack  spacing.  The  radiographs  in  Figure  43  show  the 
evenly  distributed  crack  spacing  during  each  load  step. 

At  93$  of  the  ultimate  stress  the  crack  density  was  about 
45  cracks  per  inch.  The  spacing  averaged  .022  inches 
between  transverse  cracks.  No  curved  cracks  were  seen 
until  the  last  load  3tep,  at  which  time  one  or  two  were 
found.  These  cracks  were  located  less  than  .01  inches 
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Table  8.  Crack  Angle  and  Length  Distribution  in 

Specimen  C-8,  [0/902]  ,  91  %  F 

s  tu 


Length  of  cracks  (L,  in  inches1 


Angle 

3 

1  St 

ply 

SIGHT 

2nd 

ply 

EDGE 

3rd 

ply 

4tn 

ply 

1  St 

ply 

LEFT 

2nd 

ply 

EDGE 

3rd 

ply 

4th 

ply 

COS0 

x  L 

-45 

.0 

.0 

.0 

.3 

.0 

.0 

.0 

.0 

.0 

-40 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0059 

.0022 

-35 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

-30 

.0053 

.0053 

.0 

.0 

.0 

.0 

.0055 

.0057 

.0094 

-25 

.0051 

.0051 

.0051 

.0055 

.0047 

.0108 

.0051 

.0051 

.021 1 

-20 

.0053 

.0053 

.0053 

.0102 

.0053 

.0104 

.0110 

.0104 

.0297 

-15 

.0049 

.0049 

.0102 

.0152 

.0157 

.0 

.0053 

.0104 

.0322 

-10 

.0108 

.0 

.0096 

.0098 

.0 

.0100 

.0209 

.0010 

.0306 

-  5 

.0 

.0 

.0 

.0098 

.0047 

.0148 

.0049 

.0098 

.0219 

0 

.1116 

.1384 

.1577 

.1209 

.1429 

.1504 

.1384 

.1343 

.5473 

5 

.0281 

.0248 

.0100 

.0189 

.0 

.0 

.0047 

.0047 

.0454 

10 

.0102 

.0 

.0047 

.0150 

.0051 

.0051 

.0 

.0 

.0197 

15 

.0 

.0049 

.0102 

.0047 

.0051 

.0 

.0047 

.0047 

.0120 

20 

.0055 

.0055 

.0049 

.0 

.0 

.0 

.0 

.0 

.0075 

25 

.0104 

.021 1 

.0 

.0 

.0102 

.0051 

.01 10 

.0047 

.0283 

30 

.0161 

.0075 

.0 

.0047 

.0114 

.0077 

.0 

.0059 

.0231 

35 

.0057 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0023 

40 

.0059 

.0 

.0 

.0 

.0059 

.0 

.0 

.0 

.0045 

45 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

Total 

.2262 

.2228 

.2179 

.21  48 

.2112 

.21  44 

.2116 

.21 18 

.8372 

Length 


Area  Generated  by  Transverse  Cracking 
Area  -  2  x  L  x  W  (H-  1.004  in.) 


Total  area  of  all  transverse  cracks .  1  .7376  in.’ 

Projected  area  perpendicular  to  laminate  plane  .  .  .  1.6811  in.’ 

Projected  area  parrallel  with  laminate  plane . 1791  in.’ 


sms 
x  L 


.0 

.001  9 

.0 

.0055 

.0098 

.0108 

.0086 

.0054 

.0019 

.0 

.0040 
.0035 
.0032 
.0027 
.01  32 
.0133 
.0016 
.0038 
.0 


.0892 


Each  entry  Is  the  total  length  of  all  cracks 
In  the  gage  length  at  the  specified  angle. 


from  a  straight  crack. 

Type  E  [0  2/902 ] 

9 

The  moat  notable  phenomenon  during  the  teating  of 
thia  laminate  was  the  formation  of  longitudinal  aplita. 
Theae  were  of  a  different  form  than  thoae  found  in  the 
[0/90,,]  (Type  A)  and  [0/90  3  ]_  (Type  B)  laminatea.  The 
a p 1 i 1 3  were  very  prominent,  typically  running  the  full 
length  of  the  7  inch  gage  length  (Figure  44).  Unfortu¬ 
nately,  all  of  the  aplit3  had  at  leaat  one  end  within  the 
grip,  ao  the  effect  of  the  grip  on  the  formation  of  the 
aplita  ia  uncertain.  Thia  laminate  was  the  only  one  to 
exhibit  thia  damage  form.  The  reason  thia  occurs  in  the 
[0a/90a]  laminate  and  not  in  the  [0/902l  ia  not  clear. 
Thia  laminate  was  the  only  one  tested  with  more  than  one 
0°  ply  on  the  outside  surface.  It  ia  possible  that  the 
residual  thermal  stress  state  in  this  laminate  may  be 
responsible  for  the  long  splits.  No  del  ami na t i ons  formed 
at  the  intersection  of  the  aplita  and  transverse  cracks  as 
expected.  One  would  expect  del  am i na t i ons  to  form  if  the 
laminate  were  subjected  to  cyclic  loading. 

The  splits  began  forming  at  approximately  70$  of  the 
ultimate  load.  With  the  growth  of  the  longitudinal 
splits,  a  distinct  tearing  sound  could  be  heard  (as  com¬ 
pared  to  the  "pop"  of  a  transverse  crack).  This  implies 
that  the  growth  of  the  crack  is  relatively  slow  and 
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stable. 

At  approximately  95$  of  the  failure  strength,  almost 
half  of  the  cracks  formed  could  be  considered  either 
curved  or  partial  cracks.  Figure  25  is  a  graphic  repre¬ 
sentation  of  the  formation  of  all  the  cracks.  Notice  that 
the  formation  of  straight  cracks  has  nearly  ceased  at 
about  70$  of  the  ultimate  strengh,  while  the  number  of 
curved  cracks  is  rapidly  increased  from  75$  to  failure. 
Partial  and  curved  cracks  were  counted  together.  The 
number  of  longitudinal  splits  formed  is  also  plotted. 

Table  9  was  used  to  determine  the  crack  surface  area  from 
the  measured  length  and  angle  of  cracks  in  an  edge  replica 
of  a  one  inch  gage  section. 
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Type  F  [0/90/0] 

3 

This  laminate  was  tested  to  determine  the  behavior  of 
a  single  isolated  90°  ply.  An  edge  replica  with  typical 
transverse  cracks  is  shown  in  Figure  45.  At  about  95$  of 
the  ultimate  load,  this  laminate  developed  a  transverse 
crack  density  of  65  cracks  per  inch.  Figure  46  indicates 
that  the  cracks  did  not  reach  a  saturation  level  before 
ultimate  failure.  Some  crack  branching  was  seen  in  the 
edge  replicas,  but  since  the  crack  opening  was  so  slight, 
these  were  very  difficult  to  see. 

Prior  to  failure  the  [0/90/0]  laminate  had  a 

s 

stiffness  loss  of  about  3 $ •  This  is  almost  identical  to 


j- 


Table  9.  Crack  Angle  and  Length  Distribution  in 

Specimen  E-9,  [02/902]  ,  92?  F 

s  tu 


-engtn  of  cracks  (L)  in  inches1 


Angle 

3 

1  St 

ply 

RIGHT 

2nd 

ply 

EDGE 

3rd 

ply 

v*Tl 

ply 

*  3 1 

ply 

-EFT 

2nd 

ply 

EDGE 

3rd 

ply 

4th 

ply 

COSO 

X  L 

sin  3 

X  w 

-45 

.01  1  2 

.0023 

.  0 

.  0 

.0 

.0 

.0 

.0 

.0048 

.0048 

-40 

.0 

.0071 

.0 

.0 

.0128 

.0063 

.0 

.0 

.0100 

.0084 

-35 

.01  71 

.0 

.0053 

.0053 

.0 

.0 

.0049 

.0049 

.0154 

.0108 

-30 

.0222 

.0153 

.0053 

.0104 

.0329 

.0122 

.0 

.0057 

.0450 

.0260 

-25 

.0098 

.0 

.0049 

.0 

.0049 

.0 

.0049 

.0157 

.0183 

.0086 

-20 

.0057 

.0 

.0055 

.0200 

.0 

.0047 

.0045 

.0094 

.0235 

.0086 

-15 

.0 

.0051 

.01  41 

•  0053 

.0106 

.0152 

.0 

.0053 

.0269 

.0072 

-10 

.0100 

.0102 

.0198 

.0049 

.0150 

.0148 

.0138 

.0148 

.0509 

.0090 

-  5 

.0051 

,0'-96 

.0289 

.0236 

.0049 

.01  42 

.0 

.0 

.0430 

.0037 

0 

.1118 

.1441 

.1396 

.1  404 

.1270 

.1443 

.1783 

.1  400 

.5628 

.0000 

5 

.0289 

.0238 

.0047 

.0 

.0094 

.0230 

.0090 

.01  42 

.0564 

.0049 

10 

.0194 

.01  49 

.0051 

.0 

.01  48 

.0098 

.0148 

.0047 

.0412 

.0073 

15 

.0 

.01  45 

.0051 

.0 

.0205 

.0 

.0045 

.0096 

.0263 

.0070 

20 

.01  45 

.0055 

.0165 

.0110 

.0102 

.0 

.0049 

.0049 

.0318 

.0116 

25 

.0206 

.0102 

.0047 

.0102 

.0 

.0057 

.0 

.01  1  4 

.0285 

.0133 

30 

.0 

.0 

.0 

.0059 

.0055 

.0112 

.0106 

.0185 

.0224 

.01 30 

35 

.0 

.0 

.0 

.01 16 

.0 

.0 

.0 

.0 

.0047 

.0033 

40 

.0 

.0 

.0 

.0106 

.0169 

.0039 

.0 

.0 

.0120 

.0100 

45 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0000 

.0000 

Total 

Length 

.2763 

.2626 

.2595 

.2592 

.2854 

.2653 

.2502 

.2591 

1.0239 

.1575 

Area  Generated  by  Transverse  Cracking 
Area  -  2  x  L  x  W  (W  •  1.002  In.) 


Total  area  of  all  transverse  cracks . 2.1218  In.1 

Projected  area  perpendicular  to  laminate  plane  .  .  .  2.0519  In.1 

Projected  area  parrallel  with  laminate  plane . 3156  ln.J 


Each  entry  Is  the  total  length  of  all  cracks 
In  the  gage  length  at  the  specified  angel. 


the  stiffness  loss  in  the  [0/9Q]3  even  though  it  has  a 
lower  percentage  of  0°  plies.  This  may  have  happened 
because  the  [0/90/0]  laminate  generated  a  third  more 

3 

crack  surface  area  than  the  [0/90]  laminate,  thus  making 

3 

up  for  the  original  difference  in  stiffne33. 

Type  G  [0/90/0/90] 

3 

This  laminate  was  the  most  representative  of  an  aero 
space  structure  of  all  the  laminates.  Typical  composite 
applications  usually  have  alternating  orientations  in 
adjacent  lamina,  except  when  unusual  structural  proper¬ 
ties  are  needed. 

Figure  47  shows  a  typical  edge  replica  of  this  lami¬ 
nate.  At  95$  of  the  ultimate  load,  each  single  90°  ply 
had  between  95  and  100  transverse  cracks  per  inch.  The 
double  90°  layer  had  about  60  cracks  per  inch.  These 
densities  are  higher  than  both  the  [0/90/0]3  laminate 
(single  90°)  and  the  [0/90]  laminate  (double  layer). 

This  indicates  that  a  saturation  level  had  not  been 
reached  in  these  two  laminates.  This  probably  occurred 
because  the  [0/90/0/90]  laminate  was  able  to  carry  a 
disproportionately  larger  load  than  the  two  other  lami- 
inates.  The  ratio  of  ultimate  load  to  number  of  0°  plies 
for  the  [0/90/0/90]  laminate  was  about  1750  pounds  per 

3 

ply.  This  ratio  was  about  1450  pounds  per  ply  for  the 
other  two  laminates. 


The  ability  of  the  [0/90/0/90]  laminate  to  have  such 

a  high  ultimate  strength  must  be  due  the  alternating 

stacking  sequence.  Consider  a  [02/902]  cross-ply 

s 

laminate  whose  only  damage  mode  prior  to  failure  is  trans¬ 
verse  cracking.  Rupture  of  a  0°  ply  (and  total  failure) 
probably  occurs  where  a  transverse  crack  in  a  90°  contacts 
the  0°  ply.  This  rupture  would  be  caused  by  a  stress 
riser  at  the  crack  tip.  A  qualitative  representation  of 
the  through- the- thi ckness  variation  of  the  normal  stress 
is  shown  in  Figure  48a.  The  area  under  the  curve 
represents  the  load  carried  by  the  laminate.  Now  consider 
a  [0/90/0/90]g  subjected  to  the  same  load.  The  mechanical 
properties  predicted  by  laminate  theory  are  identical. 
Figure  48b  shows  transverse  cracks  in  each  of  the  three 
90°  plies  lying  in  the  same  plane.  In  this  laminate  the 
load  is  redistributed  at  6  interfaces  rather  than  two. 

This  means  that  the  stress  concentration  at  any  interface 
will  be  smaller  than  in  the  first  case,  thus  requiring 
additional  load  to  develop  a  critical  stress. 

Laminate  Comparisons 

The  stiffness  losses  of  three  laminates  are  compared 
in  Figure  49.  The  three  laminates  have  identical  theoret¬ 
ical  mechanical  properties  (ratio  of  0°  plies  to  90°  plies 

is  1).  The  [0/90/0/90]  and  [0,/902]  laminates  have  the 

s  s 

same  number  of  plies,  but  have  different  stacking 


i  nates 


sequences.  The  [02/902]  and  [0/90]  laminates  have  the 
same  stacking  sequence,  but  a  different  number  of  plies. 
The  stiffness  loss  is  plotted  as  a  function  of  3tress 
level.  The  difference  between  the  [02/902]  and 
[0/90/0/90]s  curves  can  be  explained  by  the  constraint 
conditions  on  the  90°  plies.  Although  the  [0/90/0/90]3  has 
a  greater  crack  surface  area  at  all  stresses  (Figure  50), 
a  larger  portion  of  the  load  is  shared  by  the  90°  plies. 
The  increased  0°/90°  interface  surface  area  allows  a 
greater  proportion  of  the  load  to  be  redistributed  by 
shear  into  the  transverse  layers.  This  effectively  makes 
the  [0/90/0/90]  laminate  stiffer  at  all  stress  levels. 

3 

Figure  51  shows  the  correlation  of  crack  surface  area  to 
stiffness  loss.  The  large  difference  in  behavior  can  be 
attributed  to  the  90°  constraint  effect  which  also  affects 
the  crack  opening  displacement  in  the  transverse  cracks. 

The  constraint  conditions  in  the  [0/90]  and 

s 

[0/90/0/90]3  laminates  are  very  nearly  the  same.  In  fact, 
a  [0/90]  type  laminate  makes  up  the  core  of  the  thicker 

3 

laminate.  The  ratio  of  crack  surface  area  to  number  of 
90°  plies  at  each  stress  level  was  compared  for  the  two 
laminates,  and  it  was  found  that  the  ratio  was  only 
slightly  higher  in  the  [0/90/0/90]s  laminate.  This  can  be 
attributed  to  the  high  crack  density  in  the  isolated  90° 
plies.  All  of  these  phenomenon  support  the  observation 
that  the  two  laminates  develop  nearly  the  same  stiffness 


[0/90/0/901  and  [02/902]„  laminates 
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Figure  51.  Relationship  between  crack  surface  area  and 

stiffness  reduction  for  the  [02/902] 
and  [0/90/0/901  laminates  3 


loss  with  respect  to  stress  level. 

Laminates  A  through  D  were  designed  to  observe  the 
effects  of  90°  ply  thickness  on  laminate  response.  All 
the  laminates  contained  two  0°  plies.  The  load  at  failure 
for  this  group  was  usually  between  2600  and  3000  pounds. 
The  damage  development  in  each  laminate  has  been  described 
in  detail  in  the  previous  section. 

Several  investigators  [9,1^,20,22]  have  studied  the 
effect  of  90°  ply  thickness.  In  every  case,  it  was  found 
that  for  a  given  strain,  as  the  thickness  of  the  90°  layer 
increased,  the  crack  density  decreased.  In  the  AS4/3502 
material  system,  this  trend  was  obscured  by  the  formation 
of  curved  and  partial  cracks. 

Table  10  indicates  the  number  of  straight  cracks  ob¬ 
served  just  prior  to  ultimate  failure  for  each  laminate. 
The  [0/90/0]  laminate  can  be  considered  to  be  two 

3 

adjacent  [0/90/0]t  laminates.  This  assumption  allows  the 
comparison  with  the  other  laminates.  The  relationship 
between  straight  crack  density  and  90°  ply  thickness  is 
obvious.  When  only  straight  cracks  are  considered,  the 
familiar  relationship  is  established.  In  each  case  the 
straight  transverse  cracks  appear  to  approach  a  saturation 
crack  density.  No  instance  of  curved  crack  saturation  is 


observed . 


Table  10.  Comparison  of  Total  Crack  Density  and 

Straight  Crack  Density 


TOTAL  NO. 


NO.  OF  STRAIGHT 


TYPE 

LAY-UP 

OF  CRACKS 

CRACKS 

A 

[0/90„]s 

32 

1  0 

B 

[0/90  3  ]g 

38 

18 

C 

[0/90  2  ] 3 

45 

30 

D 

[0/90  ] 

w 

45 

43 

F 

[0/90/0] 

3 

63 

63 
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CONCLUSIONS 


The  objective  of  this  investigation  has  been  to  iden¬ 
tify  the  mechanisms  of  initiation  and  growth  of  matrix 
cracks  in  a  graphite/epoxy  laminate  and  to  identify  the 
effect  of  matrix  cracking  on  the  material  response.  This 
information  has  been  used  in  the  development  of  a  consti¬ 
tutive  model  based  on  internal  state  variable  theory.  The 
model  currently  utilizes  an  internal  3tate  variable  based 
on  the  surface  area  generated  by  matrix  cracking.  The 
response  of  a  laminate  with  an  isolated  90°  ply  has  been 
used  as  input  for  the  model;  and,  the  theoretical  predic¬ 
tions  of  the  model  have  been  compared  to  the  experimental 
results  from  several  cross-ply  laminates. 

In  this  experimental  program,  two  distinct  types  of 
transverse  cracks  were  observed.  Typical  transverse 
cracks  formed  at  regular  intervals  along  the  length  of  the 
specimen  and  generally  spanned  the  90°  layer  perpend i cular 
to  the  laminate  plane.  The  accepted  explanation  for  this 
widely  observed  phenomenon  is  the  concept  of  shear  lag. 
This  type  of  transverse  crack  was  termed  a  "straight 
crack."  The  second  type  of  crack  had  many  features  that 
indicated  that  it  possessed  a  totally  different  growth 
mechanism.  The  most  obvious  of  these  features  was  the 
incidence  angle  to  the  plane  of  the  laminate  at  which  the 
crack  initially  propagated.  This  crack  was  termed  an 
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"angle  crack."  Several  important  observations  of  this 
phenomenon  are  enumerated  as  following: 


1.  The  angle  crack  always  formed  in  conjunction  with 
a  straight  crack.  An  angle  crack  did  not  form  as 
an  isolated,  independent  crack. 

2.  The  crack  was  always  oriented  toward  the  straight 
crack. 

3.  The  angle  crack  formed  within  a  characteristic 
distance  from  a  straight  crack  or  another  angle 
crack.  This  distance  was  found  to  be  less  than  or 
equal  to  the  thickness  of  the  90°  ply  group. 

4.  The  proportion  of  angled  cracks  to  straight  cracks 
increased  with  an  increase  in  the  90°  ply  thick¬ 
ness. 

5.  Two  angle  crack  at  opposite  interfaces  would  com¬ 
monly  intersect  at  the  midplane  to  form  a  "curved 
crack . " 

6.  The  acoustic  emission  of  an  angle  crack  was  audib¬ 
ly  different  than  that  of  a  straight  crack. 

7.  Upon  unloading,  a  number  of  angled  cracks  could  be 
seen  in  intermediate  stages  of  growth.  This  was 
not  the  case  with  straight  cracks. 


Two  damage  phenomena  were  found  to  be  associated  with 
the  cracks.  Edge  delaminations  formed  as  bands  across  the 


width  of  the  specimens.  In  each  case,  the  delamination 
propagated  along  the  angle  crack  interface.  Longitudinal 
splits  formed  over  the  delamination  areas  with  om  end  of 
the  split  terminating  at  the  angle  crack  interface. 

Several  observations  were  made  concerning  laminate 
response.  A  laminate  with  alternating  0°  and  90°  plies 
developed  less  stiffness  loss  than  a  similar  laminate  with 
the  90°  plies  grouped  together.  Possible  causes  are  the 
difference  in  90°  constraint  and  the  difference  in  crack 
opening  displacement.  Another  observation  was  that  the 
strain  at  which  first  ply  failure  occurred  wa3  inversely 
related  to  the  number  of  90°  plies  grouped  together. 
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ABSTRACT 

A  continuum  mechanics  approach  is  utilized  herein  to 
develop  a  model  for  predicting  the  thermomechanical  constitution 
of  elastic  composites  subjected  to  both  monotonic  and  cyclic 
fatigue  load),  ng.  In  this  model  the  damage  is  characterized  by  a 
set  o F  second  order  tensor  valued  internal  state  variables 
representing  locally  averaged  measures  of  specific  damage  states 
such  as  matrix  cracks,  fiber-matrix  debonding,  interlaminar 
cracking,  or  any  other  damage  state.  Locally  averaged  history 
dependent  constitutive  equations  are  posed  utilizing  constraints 
imposed  from  thermodynamics  with  internal  state  variables. 

In  Part  I  the  thermodynamics  with  internal  state 
variables  is  constructed  and  it  is  shown  that  suitable 
defini tions  of  the  locally  averaged  field  variables  will  lead  to 
useful  thermodynamic  constraints  on  a  local  scale  containing 
statistically  homogeneous  damage.  Based  on  this  result  the 
Helmholtz  free  energy  is  then  expanded  in  a'  Taylor  series 
in  terms  of  strain,  temperature,  and  the  internal  state  variables 
to  obtain  the  stress-strain  relation  for  composites  with  damage. 
In  Part  II  the  three  dimensional  tensor  equations  developed  in 
Part  I  are  simplified  using  material  symmetry  constraints  and  are 
written  in  engineering  notation.  The  resulting  constitutive 
model  is  then  cast  into  laminate  equations  and  an  example  problem 
is  solved  and  compared  to  experimental  results. 

It  is  concluded  that  although  the  model  requires  further 
development  and  extensive  experimental  verification  it  may  be  a 
useful  tool  in  characterizing  the  thermomechanical  constitutive 
behavior  of  continuous  fiber  composites  with  damage. 


INTRODUCTION 


Ultimate  failure  of  composite  structural  comp'nents  is 
preceded  by  a  sequence  of  m i c r o s t r uc t ur a  1  and  ma c r o s t r u c t ur a  1 
events  such  as  microvoid  growth,  matrix  cracking,  fiber-matrix 
debonding,  interlaminar  cracking,  edge  de 1  am i nat  i  on  ,  and  fiber 
fracture,  which  are  all  loosely  termed  damage.  Considerable 
experimental  research  has  been  performed  in  the  last  decade 
detailing  the  growth  of  damage  under  both  monotonic  and  cyclic 
loading  conditions  [1-7].  The  significance  of  this  damage  lies 
in  the  fact  that  numerous  global  material  properties  such  as 
stiffness,  damping  and  residual  strength  may  be  substantially 
altered  during  the  life  of  the  component,  as  shown  in  Fig.  1  [8]. 
It  has  been  found  that  the  first  phase  of  fatigue  is  typified  by 
development  of  a  characteristic  damage  state  (CDS)  [9]  which  is 
composed  primarily  of  matrix  cracking  in  off-axis  plies.  During 
the  second  phase  of  damage  development  the  CDS  contributes  to 
fiber-matrix  debonding,  d e 1  am i n a t  i  on ,  and  fiber  microbuckl ing. 
These  phenomena  in  turn  contribute  to  a  tertiary  damage  phase  in 
which  edge  delamination  and  fiber  fracture  lead  to  ultimate 
failure  of  the  specimen  [6]. 

Analytical  modeling  of  the  damage  state  appears  to  be 
only  recently  studied.  The  earliest  attempts  fall  under  the 
general  heading  of  laminate  analysis,  in  which  various  empirical 
schemes  have  been  developed  to  discount  ply  properties  in  the 
presence  of  damage  [10-12].  Axial  stiffness  reduction  and  stress 
distribution  in  the  CDS  have  also  been  predicted  using  a  one¬ 
dimensional  shear  lag  concept  [5].  Several  researchers  have 
obtained  analytic  solutions  for  effective  moduli  of  elastic 
bodies  with  distributed  cracks  [1 3~1 7]-  In  the  case  where  cracks 
are  either  aligned  [15]  or  randomly  oriented  [13,14,16]  the 
total  crack  surface  area  is  found  to  cause  a  first  order  effect 
on  the  stiffness.  These  theoretical  results  apply  only  to 
elastic  bodies  with  cracks  of  homogeneous  and  predetermined 
dimensions. 

Fracture  based  concepts  have  recently  been  utilized  to 
model  damage  development  [18-21 j.  Although  the  first  of  these 
studies  [18]  contains  a  general'  theory  which  may  be  applied  to 
fibrous  composites,  it  has  so  far  only  been  utilized  for  quasi¬ 
isotropic  random  particulate  composites  such  as  solid  rocket 
propellant  [9],  and  as  such  has  not  been  applied  to  continuous 
fiber  composites.  Although  the  theory  in  the  latter  two  [20,21] 
has  not  been  utilized  to  predict  reduction  of  off-axis  stiffness 
components,  it  ha3  been  utilized  to  develop  fatigue  crack  growth 
laws.  Kachanov's  modulus  reduction  technique  [22]  has  also  been 
applied  to  fibrous  composites  [23]  and  although  promising  results 
were  obtained,  the  model  was  utilized  in  uniaxial  form  only. 

A  complex  interactive  experiment  and  analysis  model 
(called  a  mechanistic  model)  has  been  proposed  [8]  for  prediction 
of  life  of  damaged  composites.  The  mechanistic  model  appears  to 
require  numerous  experimental  results  for  each  geometric  layup  in 
order  to  determine  which  damage  mode  results  in  failure. 

Finally,  extended  forms  of  Miner's  rule  [24]  have  been 
proposed  for  life  prediction  [18,25].  However,  they  are  based  on 
simplified  microphysical  models  at  this  time. 


RESIDUAL  STRENGTH 
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TIME 


Damage  Accumulation  in  a  Continuous  Fiber 
Composite  Subjected  to  Monotonic  Load  or 
Strain  Controlled  Cyclic  Fatigue  (from  ref.  8) 


The  concept  of  damage  as  an  internal  state  variable  [26] 
has  been  previously  utilized  in  continuum 
m e c h an i c s / t her  mod yn am i c s  based  theories  for  crystalline  and/or 
brittle  materials  [27-3^],  as  well  as  for  nonlinear  viscoelastic 
materials  [18].  A  study  has  been  made  of  the  effect  of  vector¬ 
valued  damage  parameters  on  various  compliance  terms  [35],  and 
this  methodology  is  currently  undergoin;  further  development  [36, 
37]. 

The  research  reviewed  above  indicates  that  although 
important  progress  has  been  made  in  characterizing  damage  in 
composites,  substantial  and  continued  research  is  warranted 
before  several  issues  can  be  resolved. 

In  this  paper  an  attempt  will  be  made  to  assemble  many  of 
the  concepts  embodied  in  the  research  efforts  mentioned  briefly 
above  and  to  utilize  these  concepts  to  develop  a  thermomechanical 
constitutive  model  for  damage  in  composites  which  is  rigorously 
based  in  continuum  mechani os/ thermodynami cs  and  is  generic  with 
regard  to  material  type,  load  spectrum,  and  specimen  geometry. 

The  model  will  utilize  the  concept  of  a  local  volume  element 
with  statistically  homogeneous  damage  to  construct  constitutive 
equations  relating  stress,  strain,  and  damage.  Unlike  current 
analytical  methods,  the  local  volume  element  will  be  modelled 
experimentally.  The  model  will  therefore  not  be  restricted  to 
linear  elastic  media  with  homogeneous  elastic  properties. 
Furthermore,  the  model  will  be  applicable  to  cracks  which  are 
oriented  and  of  heterogeneous  and  irregular  size  and  shape.  The 
effect  of  the  cracks  will  be  reflected  through  locally  averaged 
quantities  describing  the  kinematics  of  the  cracks. 

The  output  of  the  model  will  be  a  set  of  constitutive 
equations  which  apply  on  a  scale  that  is  small  compared  to  the 
boundary  value  problem  of  interest.  Therefore,  it  will  be 
applicable  to  the  analysis  of  bodies  with  stress  gradients  and 
heterogeneous  damage  states  as  long  as  the  scale  of  the 
microcracks,  voids,  etc.  is  relatively  small. 


CHARACTERIZATION  OF  DAMAGE  AS  A 
SET  OF  INTERNAL  STATE  VARIABLES 

Consider  an  initially  unloaded  and  undamaged  composite 
structural  component,  denoted  B,  as  shown  in  Fig.  2a,  where 
undamaged  is  defined  here  to  mean  that  the  body  may  be 
considered  to  be  continuous  (without  cracks)  on  a  scale  several 
orders  of  magnitude  smaller  than  the  smallest  external  dimension 
of  the  component.  Although  cracks  may  exist  in  the  initial 
state,  their  total  surface  area  is  assumed  to  be  small  compared 
to  the  external  surface  area  of  the  component.  Under  this 
assumption  the  body  is  assumed  to  be  simply  connected  and  we  call 
the  initial  bounding  surface  the  external  boundary  S.  Although 
the  component  is  undamaged,  there  may  exist  local  heterogeneity 
caused  by  processing  and  second  phase  materials  including  fibers, 
matrix  tougheners  and  voids.  In  addition,  the  body  may  be 
subjected  to  some  residual  stress  state  due  to  processing,  cool 
down,  etc. 

Now  suppose  that  the  component  is  subjected  to  some 
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traction  and/or  deformation  history,  as  shown  in  Fig.  2b.  The 
specimen  will  undergo  a  thermodynamic  process  which  will  in 
general  be  in  some  measure  irreversible.  This  irreversibility  is 
introduced  by  the  occurrence  of  such  phenomena  as  material 
inelasticity  (even  in  the  absence  of  damage),  fracture  (both 
micro-  and  macroscale),  friction  (due  to  rubbing  and/or  slapping 
of  fractured  surfaces),  temperature  flux,  and  chemical  change. 
While  all  of  these  phenomena  can  and  do  commonly  occur  in 
composites,  in  the  present  research  it  will  be  assumed  that  all 
irreversible  phenomena  of  significance  occur  in  small  zones  near 
crack  surfaces.  Outside  these  zones,  the  behavior  will  be 
considered  to  be  elastic  and  therefore  reversible  under  constant 
temperature  conditions.  All  fracture  events  will  be  termed 
damage.  Due  to  these  fracture  events,  the  body  will  necessarily 
become  multiply  connected,  and  all  newly  created  surfaces  not 
intersecting  the  external  boundary  will  be  termed  internal 
boundaries.  Because  of  the  above  assumptions  the  model  may  be 
limited  to  polymeric  and  ceramic  matrix  composites  at 
temperatures  well  below  the  glass  transition  temperature  T  or 
melting  temperature,  where  viscoelasticity  in  matrix  materials  is 
small.  Metal  matrix  composites  may  have  to  be  excluded  due  to 
complex  post-yield  behavior  of  the  matrix. 

While  fracture  involves  changes  in  the  boundary  conditions 
governing  a  complex  field  problem,  it  is  hypothesized  that  one 
may  neglect  boundary  condition  changes  caused  by  creation  and 
alteration  of  both  internal  and  external  surfaces  created  during 
fracture  as  long  as  the  resulting  damage  in  the  specimen  is 
statistically  homogeneous  on  a  local  scale  which  is  small 
compared  to  the  scale  of  the  body  of  interest.  However,  the  total 
newly  created  surface  area  (which  includes  internal  surfaces)  may 
be  large  compared  to  the  original  external  surface  area.  Under 
the  condition  of  small  scale  statistical  homogeneity  all 
continuum  based  conservation  laws  are  assumed  to  be  valid  on  a 
global  scale  in  the  sense  that  all  changes  in  the  continuum 
problem  resulting  from  internal  damage  are  reflected  only  through 
alterations  in  constitutive  behavior.  Typical  m i c r o s t r u c t ur a  1 
events  which  may  qualify  as  damage  are  matrix  cracking  in  lamina, 
fiber/matrix  debonding,  localized  interlaminar  delamination  and 
fiber  fracture.  Large  scale  changes  in  the  external  surface  such 
as  edge  de 1  ami na t i ons ,  however,  are  treated  as  boundary  effects 
which  must  be  reflected  in  conservation  laws  via  changes  in  the 
external  boundary  conditions  rather  than  in  constitutive 
equations  [35,38]. 


THERMODYNAMICS  OF  MEDIA 
WITH  DAMAGE 

We  now  proceed  to  construct  a  concise  model  of  the 
composite  with  damage.  To  do  this,  consider  once  again  the 
structural  component,  denoted  B  in  Fig.  2a.  The  body  B  is 
assumed  to  be  of  the  scale  of  some  appropriate  boundary  value 
problem  of  interest.  Now  consider  some  local  element  labelled 
and  with  external  surface  faces  S1  arbitrarily  chosen  normal  to  a 
set  of  Cartesian  coordinate  axes  (x1t  x2>  x^)»  as  shown  in  Fig. 


2c.  The  element  V  L  extracted  from  B  and  the  newly  created 
surfaces,  denoted  S2  and  with  volume  VQ,  are  subjected  to 
appropriate  boundary  conditions  so  that  the  element  response  is 
identical  to  that  when  it  is  in  B.  Furthermore,  the  volume  of  the 
element  is  defined  to  be  V^,  which  includes  the  volume  of  any 
initial  voids.  The  scale  of  is  chosen  so  that  its  dimensions 
are  small  compared  to  the  dimensions  of  B,  but  at  the  same  time, 
the  dimensions  of  are  large  enough  to  guarantee  statistical 
homogeneity  of  the  material  heterogeneities  and  defects  in  V , 
even  though  the  total  surface  area  of  defects  may  be  of  the  same 
order  of  magnitude  as  [39].  Suppose  furthermore  that  in  the 
absence  of  defects  or  at  constant  damage  state  the  material 
behavior  is  linearly  thermoelastic.  Now  consider  the  local  volume 
element  V^.  For  the  case  where  tractions  or  displacements  are 
applied  uniformly  to  the  external  boundary  of  V ^ ,  the  average 
stresses  and  strains  in  VL  will  be  determinable  from  the  external 
boundary  tractions  or  displacements. 

Although  the  damage  process  actually  involves  the  conversion 
of  strain  energy  to  surface  energy,  the  fact  that  the  damage  is 
reflected  in  the  local  constitutive  equations  rather  than 
boundary  conditions  suggests  that  it  be  treated  as  a  set  of 
energy  dissipative  local  state  variables  which  are  not 
discernible  on  the  external  boundary  of  the  local  element. 


Review  of  Thermodynamic  Constraints  on  Linear  Thermoelastic  Media 

The  following  notation  is  adopted.  Quantities  without 
capitalized  subscripts  denote  pointwise  quantities.  Those  with 
subscripts  L  denote  quantities  which  are  averaged  over  the  local 
element  V^.  Finally,  the  subscript  E  denotes  linear 
thermoelastic  properties. 

Under  the  conditions  described  in  the  previous  section 
the  pointwise  Helmholtz  free  energy  per  unit  volume  h  of  the 
undamaged  linear  elastic  medium  may  be  expressed  as  a  second 
order  expansion  in  terras  of  strain  and  temperature  T  as 
follows  [40]:  "J 


h  =  u  -  Ts  «  h  ( ei j , T) - 

*  *  *  2  CijklEijEkl  *  D4T  *  *  2  F4T2  •  <') 

where  u  and  3  are  the  internal  energy  and  entropy  per  unit 
volume,  respectively,  and  A,  B^,  c  j.  j  1  »  D,  Eij  and  F  are 
material  parameters  which  are  independent  of  strain  and 
temperature  and  AT  et-Tr,  where  TR  is  the  reference  temperature 
at  which  the  strains  are  zero  at  zero  external  loads.  In 
addition,  we  assume  here  that  all  motions  are  associated  with 
small  deformations.  Furthermore,  inertial  effects  and 
electromagnetic  coupling  are  assumed  to  be  negligible. 

Pointwise  conservation  laws  appropriate  to  the  body  are 
as  follows: 


1)  conservation  of  linear  momentum 


Thermodynamic  Constraints  with  Local  Damage 

It  is  our  intention  to  construct  locally  averaged  field 
equations  which  are  similar  in  form  to  the  pointwise  field 
equations  discussed  above.  In  performing  this  averaging  process 
the  pointwise  Helashol  tz  free  energy  described  in  equation  (1) 
will  undergo  a  natural  modification  to  include  the  energy 
conversion  due  to  crack  formation. 

Now  consider  the  local  element  shown  in  Fig.  2c  with 
traction  boundary  conditions  on  the  external  surface  S-j.  In 
addition,  the  interior  of  V  L  is  assumed  to  be  composed  entirely 
of  linear  elastic  material  and  cracks  (which  may  include  thin 
surface  layers  of  damage).  Integrating  pointwise  equations  (1) 
through  (6)  over  the  local  volume  will  result  in 

hEL  “  AL  +  BLijeLij+2CLijklGLij  eLkl  +  DLATL+ELije  LijATL  +  2FLATL2  ’(11) 

where  A  ^ ,  BLij,  cLijkl'  °L’  ELij’  and  FL  are  l°cally  averaged 
material  constants.  J Also, 


°LJi.J  “  0 

°Lij  ‘  aL  j  i 

UL  "  aL  i  j  eL  i  j  +  qLj  ,  j  "  rL 

and 


(12) 

(13) 


sL  -  j  L  ♦  (f  Lj  )  ,  J  >_  0 
L 


(15) 


where  u£,  called  the  effective  local  internal  energy,  is  given  by 


u£L  represents  the  internal  energy  of  the  equivalent  uncracked 
body,  given  by 


if  -  if  T?“i 


(17) 


where  are  called  equivalent  tractions,  representing  tractions 
in  the  uncracked  body  acting  along  fictitious  crack  faces,  as 
described  in  detail  in  the  appendix,  and  u£  is  the  mechanical 
power  output  due  to  cracking,  given  by 


(18) 


where  T?  are  fictitious  tractions  applied  to  the  crack  faces 
which  represent  the  difference  between  the  actual  crack  face 
tractions  and  T^.  Furthermore,  the  locally  averaged  stress  is 
given  by 


hf° 


(19) 


and  the  locally  averaged  strain  is  given  by 


VLJ  2(uinj  +  ujn: 


(20) 


Equations  (11)  through  (15)  are  identical  in  form  to  equations 
(1)  through  (5)  ,  respectively.  Further  details  on  this 
similarity  are  given  in  the  appendix. 

On  the  basis  of  this  similarity  we  now  define  the  locally 
averaged  Helmholtz  free  energy  [18,  38]: 


TLsL 


lEL  ~  TLsL  +  UL 


hEL  +  UL 


(21  ) 


where  it  can  be  seen  from  definition  (17)  that  hEL  is  the  locally 
averaged  elastic  Helmholtz  free  energy  for  which  residual  damage 
is  zero. 

The  similarity  between  the  pointwise  and  local  field 
equations  leads  to  the  conclusion  that 
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gLj  =  TL , j 


where 


(25) 
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Note  the  similarity  between  equations  (7)  through  (10)  and  (22) 
through  (25),  respectively. 

Equations  (23)  will  nerve  as  the  basis  for 
thermomechanical  stress-strain  relations  in  damaged  composites. 
All  damage  will  be  reflected  through  the  local  energy  due  to 
cracking  u£.  This  term  will  be  modelled  with  internal  state 
variables  characterizing  the  various  damage  modes. 


Description  of  the  Internal  State 

In  order  to  describe  the  internal  state,  we  first 
consider  the  kinematics  of  a  typical  point  0  with  neighboring 
points  A  and  B,  as  shown  in  Fig.  3.  Before  deformation  lines  0A 
and  OB  are  orthogonal,  as  shown  in  (a).  After  deformation  we 
imagine  that  lines  joining  O', A',  and  B*  are  as  shown  in  (b),  d 
just  at  the  instant  that  deformation  is  completed,  a  crack  forms 
normal  to  the  plane  of  AOB  through  point  O',  as  shown  in  (c). 
Furthermore,  point  O’  becomes  two  material  points  O’  and  0"  on 
opposite  crack  faces  and  points  A'  and  B'  deform  further  to 


points  A"  and  B".  It  is  assumed  that  all  displacements, 
including  displacement  jumps  across  crack  faces,  are 
infinitesimal,  so  that  strain  gages  attached  at  points  0,  A,  and 
B  record  only  the  deformation  A"0'B".  However,  the  actual  strain 
is  associated  with  A"0"B".  Therefore,  it  is  essential  to 
construct  an  internal  state  variable  which  will  relate  these  two 
strain  descriptions.  We  therefore  construct  the  vectors  uc 
connecting  O'  and  0”  and  nG  describing  the  normal  to  the  crack 
face  at  O',  as  shown  in  (c).  It  should  be  noted  that  uG  can  be 
used  to  construct  a  pseudo-strain  representing  the  difference  in 
rotation  and  extension  of  lines  A"0'B"  and  A"0"B". 

Now  recall  that  the  mechanical  power  output  during 
cracking  is  given  by  equation  ( 1 8 ) .  We  assume  that  at  any  point 
in  time  t1  tractions  T^  can  be  applied  along  the  crack  faces 
which  will  result  in  an  energy  equivalent  to  that  produced  by  the 
damage  process: 


u£( ti ) 


V  l  fT  iui' 


(27) 


S  2  (  t  -j  ) 


The  quantities  T,  do  not  necessarily  coincide  with  the  terms  in 
the  integrand  or  (18)  since  the  process  is  in  some  measure 
irreversible.  However,  we  define  them  such  that  the  total 
energies  in  equations  (18)  and  (27)  are  equivalent.  For 
convenience  we  will  call  them  crack  closure  tractions,  although 
they  do  not  necessarily  result  in  complete  crack  closure. 

Guided  by  the  fact  that  u°  and  nG  describe  the  kinematics 
of  the  cracking  process  at  point  0,  we  now  define  the  following 
second  order  tensor  valued  internal  state  variable: 


ai j  s  uinj  =>  [ctij] 


uGnG  uGn|  uGn^ 
u2 n G  u^n^  u^n^ 
u^nG  u|n|  u^n^ 


(28) 


The  above  description  has  been  previously  proposed  by  M. 
Kachanov[4l  ].  Substituting  the  above  into  (27)  and  utilizing 
Cauchy's  formula  gives 
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«  <dS 


°i  j  ai  j 


where  it  should  be  pointed  out  that  integration  is  performed  with 
respect  to  undeformed  coordinates. 

Note  that  the  components  of  uc  can  be  recovered  from  (28) 
by  using  simple  row  multiplication  on  aij: 
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u^  -  u?njU°rij  (no  sum  on  i)  .  (30) 

Similarly,  nc  can  be  recovered  by  using  column  multiplication  on 
°ij  : 

n2.  == -^c.2  (no  sum  on  j  )  .  (31  ) 

3  (u  )  i  j  i  j 

Therefore,  although  it  would  not  be  necessary  to  actually  perform 
the  operations  described  in  equations  (30)  and  (31).  the  normal 
and  shear  modes  of  crack  displacement  can  be  recovered  from 

Note  furthermore  that  is  generally  an  asymmetric  tensor, 

and  that  a  symmetric  alternative  to  equations  (28)  could  not  be 
utilized  to  recover  normal  and  shear  modes  as  described  in  (30) 
and  (31).  As  an  example,  consider  the  following  decomposition  of 
(28)  into  symmetric  and  anti-symmetric  components 
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In  order  for  the  anti-symmetric  tensor  w2ij  t0  be  zero»  uc  and 
must  be  parallel  vectors,  implying  pure  made  I  fracture.  In  this 
case  wiij  could  be  decomposed  into  a  vector  (in  local 
coordinates),  thus  resulting  in  v e c t or - v a  1 ue d  internal  state 
variables.  For  the  case  where  the  cracks  in  the  local  volume  VL 
are  randomly  oriented  and  of  statistically  homogeneous  shape  and 
3ize,  the  surface  integral  in  equation  (29)  may  be  carried  out 
over  all  cracks.  However,  if  various  groups  of  cracks  in  V,  are 
distinguished  by  markedly  different  crack  normals  n^3  or 
geometries,  then  it  will  be  necessary  to  distinguish  between  the 
damage  modes  in  order  to  retain  the  kinematic  features  of  the 
damage  process.  Therefore,  define  the  locally  averaged  internal 
state  var iab 1 e' a  £ ^ j  for  the  nth  damage  mode  as  follows: 


u°nj dS 


(  35) 
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where 


E  s2" 


and  N  is  the  number  of  damage  modes.  For  a  continuous  fiber 
laminated  composite,  the  modes  might'  be  represented  by  matrix 
cracks,  interply  delamination,  fiber  fracture,  and  fiber-matrix 
debond  (N-4).  For  a  quasi-isotropic  chopped-f i ber  metal  matrix 
composite,  a  single  isotropic  damage  tensor  might  suffice  for 
randomly  oriented  matrix  cracking  (N-1). 

Therefore,  if  we  define  a£ij  to  be  the  average  crack  closure 
stress  for  the  nth  damage  mode  such  that 


acn  an 

1,1  ^  T  I 


Lkl  Lkl  - 
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ij  ij 


it  follows  from  equations  (29),  (35),  (36),  and  (37)  that 


aL?jaZij 


(38) 


where  we  have  assumed  that  repeated  indeces  n  imply  summation 
over  the  range  N.  It  is  clear  from  the  above  discussion  that  the 
value  of  N  must  be  sufficiently  large  to  recover  the  essential 
physics  of  the  damage  process.  In  a  mathematical  sense,  this 
implies  that,  whereas  the  mapping  from  ctji  to  oc1^,  i  j  is  unique, 
the  inverse  should  also  be  true  in  an  approximate  sense. 
However,  there  is  no  clearcut  definition  for  the  range  N  which 
will  lead  to  an  accurate  description  of  the  internal  damage 
state.  Due  to  the  inclusion  of  crack  opening  displacements  in 
the  description  of  the  internal  state  variables  a^ii’  these 
quantities  should  not  be  interpreted  strictly  as  ""damage" 
parameters.  Since  cracks  may  actually  close  on  unloading 
(without  healing),  ct£ij  should  be  interpreted  as  damage 
parameters  only  during  loaJdup.  Note  also  that  both  u^^  and  n,  in 

equations  (35)  will  be  affected  by  crack  interaction  in  the  local 

vol ume . 

As  an  example,  consider  the  case  of  mode  I  opening  of  an 
elliptic  crack.  For  this  case,  equation  (35)  will  result  in 
dependence  of  on  the  volume  of  the  inclusion.  Although 

analytic  models  for  linear  elastic  bodies  with  cracks  result  in 
response  which  is  dependent  on  the  surface  area  of  cracks  only 

[14-16],  it  should  be  pointed  out  that  they  also  require  the 

average  crack  diameter.  This  quantity  is  replaced  herein  by  the 
crack  opening  displacement,  which  is  proportional  to  the  crack 
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diameter  in  a  linear  elastic  body.  Therefore,  specifying  the 
crack  opening  displacement  is  equivalent  to  specifying  the  crack 
diameter . 

Now  consider  equation  (38)  in  further  detail.  The 
kinetic  quantities  °^  •  :  may  be  interpreted  as  generalized 
stresses  which  are  energy  conjugates  to  the  kinematic  strain-like 
internal  state  variables  ct^  ji  •  We  infer  from  this  that  there 
exists  a  constitutive  relation  between  these  variables  of  the 
form 


0cn  0cn  /  e  T  a  ^  \ 

Lij  Lij'’  Lkl *  1 L ’  Lkl ' 


which  is  history  dependent  via  the  explicit  dependence  on  the 
internal  state  variables. 

Therefore,  substituting  (39)  into  (38)  will  give 


uL^t1^  -Ju£(t)dt  «  u£  (  eLkl  ^  1  ^  *  TL^ 1 1  ^  *  aLkl  ^ 1 1  ^ 


.  (  40) 


It  is  now  proposed  that  u£  be  expanded  in  a  Taylor  series 
which  is  second  order  in  each  of  the  arguments  in  equation  (40) 
as  foil ows : 
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(  41  ) 


where  all  terms  are  at  least  linear  in 
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that  u  £  depends  explicitly  on  damage,  and  ATleTl~Tr.  Thus, 
substituting  (11)  and  (41)  into  equations  (23)  and  neglecting 
higher  order  terms  yields: 


Lij 


BL i j  +  ELij  ^L  +  CLijkleLkl  +  IijklaLkl 


(  42) 


Equations  (42)  may  be  written  in  the  following  alternate  form  for 
isothermal  conditions 


Lij 


a +  n  *  £ 

Lij  uLijkl  Lkl 


,  (43) 


where 


°Lij  H  BLij 


is  the  residual  stress  tensor;  and 


(44) 


CL i j  kl ekl 


—  n  £  ,  T  ^  n  ^ 

CL  i  j  kl  Lkl  ILijklaLkl 


(  45) 


defines  the  effective  modulus  tensor  for  any  damage  state. 
Note  that  although  equations  (43)  are  similar  to  Kachanov's  model 
[22],  the  stiffness  reduction  is  a  first  order  effect  of  damage. 
Note  also  that  the  inclusion  of  higher  order  terms  will  result  in 
damage  dependent  residual  and  thermal  stresses,  as  well  as 
nonlinear  stiffness  loss  as  a  function  of  damage. 

Equations  (42)  are  the  completed  description  of  the 
stress-strain  relationship.  Note  that  these  equations  reduce  to 
the  standard  linear  thermo  elastic  equations  in  the  absence  of 
damage  (  ±  j  -  0  ) . 

Damage  Growth  Laws 

The  model  is  completed  with  the  construction  of  the 
damage  growth  laws,  which  may  be  described  in  the  following 
differential  equation  form: 


•  n 
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,  (46) 


or  equivalently,  when  ^ j  are  single  valued  functions  of  time, 
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(t1  ) 


fiij  ^  £Lkl  ^ ^  *  Tl/  ^  *  “Lkl  ^  ^  )dt 


(  47) 


Although  the  above  equations  are  called  "growth"  „  aws  they 
have  the  more  general  capability  to  model  su  :h  phenomena  as 
heal ing . 

The  precise  nature  of  equations  (47)  is  determinable  only 
through  a  concise  experimental  program  coupled  with  an 
understanding  of  the  micromechanics  of  the  medium.  Indeed,  these 
growth  laws  constitute  the  single  most  complex  link  in  the  model 
development . 

In  this  section  an  example  of  a  first  generation  growth 
law  will  be  constructed  for  predicting  damage  up  to  the  CDS  in 
continuous  fiber  composites.  Experimental  evidence  suggests  that 
matrix  cracks  dominate  development  of  the  CDS  [4-6].  Guided  by 
this  observation,  a  single  damage  tensor  is  considered  in  this 
section:  Lij  rePresentin8  matrix  cracking. 

In  order  to  completely  define  equations  (47),  it  is 
necessary  to  construct  indicators  of  both  the  magnitude  and 
direction  of  the  damage  tensor.  In  this  first  generation  model 
it  is  assumed  that  the  direction  of  the  damage  tensor  is  known  a 
priori  and  does  not  vary  as  the  damage  state  changes. 
Specifically,  in  a  typical  laminate,  it  is  assumed  that,  for  this 
simple  example,  in  accordance  with  equation  (35),  the  locally 
averaged  resultants  of  uc  and  nc  are  normal  to  the  fiber 
direction  in  each  ply,  as  shown  in  Fig.  4.  Thus,  for  example,  in 
a  0°  ply  122  /  and  a1^  other  components  are  zero,  whereas 
in  a  90°  ply,  .  /  0,  and  all  other  components  are  zero  (in 
global  coordinates).  In  Part  II  a  somewhat  more  general  case  of 
the  damage  state  for  matrix  cracking  will  be  considered. 

Under  the  above  assumptions,  the  magnitude  of  the  damage 
tensor  is  the  sole  repository  for  history  dependence  in  each  ply. 
Experimental  evidence  indicates  that  for  matrix  cracking  in 
randomly  oriented  particulate  composites  [42]  and  matrix  cracks 
in  fibrous  composites  [20,21]  the  growth  of  damage  surface  area 
is  related  to  the  energy  release  rate  G  by 
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,  (48) 


where  S2  represents  crack  area,  N  represents  the  number  of  cycles 
in  a  fatigue  test,  and  n  is  some  material  parameter.  Guided  by 
these  results,  a  similar  law  is  constructed  here.  Equation  (48) 
may  be  rewritten  in  the  following  form: 
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so  that  it  follows  that 
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Assuming  that  the  energy  release  rate  is  essentially  mode 
I  and  therefore  depends  on  the  maximum  normal  strain,  the  damage 
growth  law  for  matrix  cracking  is  thus  hypothesized  to  be  of  the 


form 

® [22  - 

aL22 

'  dtn 

*  f  enmin  <  £n 

,  and 

“L22  =  k2en 

^  Gnmin  —  £n 

t 

(51  ) 

> 

> 

V 


T 


is 

9 


V 


J? 


where  en  is  the  local  normal  strain  component  which  is  normal  to 
the  fibers.  Furthermore,  enmin  is  the  value  of  en  at  which 
matrix  cracking  initiates.  k1f  k2,  and  n  are  experimentally 
determined  material  parameters  which  may  crepend  on  the  initial 
damage  state  or  on  history  dependent  damage  other  than  matrix 
cracks.  The  use  of  en  presupposes  that  the  fracture  mode  is 
predominantly  mode  I  in  nature,  which  may  not  be  the  case  in  some 
complex  layups.  In  these  cases,  mode  II  and  mode  III  terms  may 
be  required.  Note  that  all  components  of  “^ij  are  zero  except 
L22  *  which  is  nonzero  in  the  local  ply  coordinate  system 

wherein  the  fibers  are  aligned  parallel  to  the  local  x-|  axis. 

Experimental  evidence  [43]  indicates  that  in  crossply 
laminates  with  multiple  adjacent  crossplies  in  sequence,  it  is 
not  uncommon  to  observe  matrix  cracks  which  are  curved  rather 
than  normal  to  the  plane  of  the  ply.  For  these  cases  it  is 
necessary  to  carry  components  of  aL,ij  in  both  the  x2  and  x^ 
coordinate  directions.  Although  it  is  hypothesized  that  these 
components  may  perhaps  be  determinable  from  the  orientation  of 
the  maximum  normal  strain  £n ,  this  issue  is  under  further 
investigation  by  the  authors. 

Equations  (51)  complete  the  description  of  the  damage 
model  for  the  case  of  matrix  cracking.  Integration  of  these 
equations  in  time  will  lead  to  current  values  of  the  damage 
tensor  which  is  input  to  constitutive  equations  (42).  Fig.  5 
shows  a  typical  growth  history  for  a  specimen  subjected  to 
monotonical 1 y  increasing  deformation  u(L).  It  should  be  pointed 
out,  however,  that  these  equations  may  be  extremely  nonlinear  and 
as  such  must  in  some  cases  be  integrated  numerically  with  stiff 
integration  schemes  [44]. 


S’ 

.<  CONCLUSION 

v  Stress-strain  relations  have  been  developed  herein  which 

^  account  for  various  forms  of  damage  in  continuous  fiber 

composites.  Furthermore,  a  damage  growth  law  has  been  proposed 
for  matrix  cracking  in  fibrous  composites.  The  model  developed 
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The  i  o  t  u.t  1 


L  are)  . 
t  h  f.  a 


1  e in ;>  l  i  oated 


r  e  q  u  1  r  erne  n  t  p^r  n  u  ;ne  r’  o  u  a  ex  per  1  rner.L-j  1  1.  y  d  u  l  e  r  4  ru-  i  quanti  t  *  os, 
we  1  1  is  the  necessity  to  determine  1  o<:«)  l  y  bused  eoser  vabl  e  state 
v  a  r  4  a  b  l  e  s  L  y  a  na  1  y  t  >  e  in  ethods.  The  c  o  n  s  t  r  u  ot  ion  o  f  t  a  e  s  e 
parameters  constitutes  on  entire  separate  research  effort  wh  4  ch 
is  considered  mi  Part  II. 

Finally,  it  should  be  pointed  oat  that  although  an 
internal  state  variable  growth  law  has  been  proposed  herein  only 
''or  matrix  cracks,  the  model  is  in  principal  applicable  to  more 
complex  damage  states  in  laminated  composites,  and  research  is 
underway  to  consider  other  damage  modes. 
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APPENDIX 


Consider  a  local  volume  element  with  some  damage  state, 
where  the  crack  faces  are  defined  as  traction  free  surfaces,  as 
shown  in  Fig.  A1  (al).  For  convenience  we  show  only  one  crack, 
although  in'  actuality  the  damage  must  be  statistically 
homogeneous  in  V^.  Now  replace  the  actual  cracks  with  fictitious 
cracks  which  a;  e  described  by  the  bounding  surface  between 
elastic  and  inelastic  response  near  cracks,  as  shown  in  (a2).  We 
define  this  surface  to  be  S2-  In  order  to  insure  that  the  total 
mechanical  states  in  the'  two  systems  are  identical,  the 
fictitious  case  must  include  tractions  labelled  on  $2* 

Now  suppose  that  VL  is  subjected  to  boundary  tractions  on 
S  1  _i  _n  _t  _h  e_  undamaged  state  as  shown  in  (bl).  We  define  an 
equivalent  elastic  problem  in  which  the  surface  S2  described  in 
(a2)  is  cut  from  V L  and  elastic  tractions  are  applied  on  S2  so 
that  the  total  mechanical  states  of  systems  (bl)  and  (b2)  are 
equivalent. 

The  actual  system  of  interest  is  described  in  (al). 
However,  for  pragmatic  reasons  we  wish  to  replace  the  actual 
system  with  a  fictitious  system  with  equivalent  mechanical  state. 
To  do  this,  we  first  replace  (al)  with  (a2),  which  by  definition 
has  equivalent  mechanical  state.  Next,  we  define  a  system 
equivalent  to  (a2),  such  that 


T=  =  if  -  T?  ->  Tj  -  T?  ♦  T?  ,  (Al) 

as  shown  in  (c).  Integrating  the  balance  of  energy  (4)  over 
the  local  volume  and  dividing  through  by  the  local  volume  results 
in 


ij  ^jS^dv  * 

Jvl 


( A2  ) 


Now  consider  the  second  term  in  equation  (A2).  Recall 
that  since  is  a  symmetric  tensor 


ai  j  e  i  j  "  2°ij  (ui  ,  j+uj  ,i  }  “  aijui,j 


(A3) 


Thus,  assuming  that  the  stresses  are  negligible  in  Vc,  the  volume 
enclosed  by  S2,  using  the  divergence  theorem  and  substituting 
Cauchy’s  formula  gives 


inelastic  zone 


Kf° 


•Jijeljav 


\J° ij"t ,  jav  ■  \JTi  “ids  *  1L/Ji 


dS 


VL'Vc 
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Jr t ut dS  +  J  JrfuidS  +  J  ^tJuj 


dS 


( Ail) 


where  n,  are  the  components  of  the  unit  outer  normal  vector  to 
the  surface  S  -  S-j  +  $2*  Now  define 


u 


c 

L 


v/T 


T?uidS 


( A5) 


which  is  the  effective  specific  mechanical  power  output  of  the 
continuum  due  to  the  crack  surface  tractions.  This  term  contains 
both  the  mechanical  power  due  to  crack  extension  as  well  as  the 
mechanical  power  due  to  apparent  stiffness  loss  caused  by 
existing  cracks.  For  the  special  case  of  a  reversible  process 
this  is  the  time  rate  of  change  of  surface  energy  release  per 
unit  local  volume  due  to  cracking  in  V L.  Furthermore,  define 
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Lij 


+  Uj  n^ )  dS 


(A6) 
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and 


°L  i  j 


aimnmxj+  ajmnmxi>dS 


( A7) 


Therefore,  for  the  case  of  either  spacial 1 y  uniform  surface 
tractions  or  displacements  which  are  linear  in  coordinates  on  S1 
one  readily  obtains 


°Lij  ^ij 


“injdS 


dS 


(AS) 


Although  it  will  be  assume d  in  the  remainder  of  this  paper  that 
the  above  conditions  are  satisfied,  they  need  only  be 
approximately  true  if  is  statistically  homogeneous.  Thus, 
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equation  (A4)  becomes 


Kf° 


°ljeijav 


°L  i  j  e  L  i  j  +  V 


1l/!S 


.  (  A  9  ) 


Define  also 


qLi,j  = 


vL  Jqinj 


,  (A10) 


.  (All) 


Now  define 


UEL  H 


,  (A12) 


which  can  be  seen  from  Fig.  A1  (b)  to  be  the  equivalent  internal 
energy  rate  that  would  be  produced  in  the  body  without  cracks. 
Note  that  Ugg  is  not  path  dependent  since  it  represents  elastic 
response.  Substituting  equations  (A4),  (A5)»(A9),  (A10),  (All), 
and  ( A 1 2)  into  equation  (A2)  yields  the  following  locally 
averaged  balance  of  energy: 


UEL  +  UL  aL  i  j  eL  i  j  +  qLj,j  “  rL 


,  (  A 1  3  ) 


We  now  define  the  effective  internal  energy  u^  (which  may  be  path 
dependent)  such  that 


UL  =  UEL  +  UL 


Substitution  of  (A14)  into  ( A 1 3 )  results  in 


,  (  A 1  4 ) 


UL  ‘  °L  i  j  e  L  i  J  +  qLJ,j  ’  rL 


,  (  A 1  5  ) 
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which  can  be  seen  to  be  equivalent  in  form  to  energy  balance  law 
(4). 

In  order  to  construct  a  similar  statement  for  entropy 
production  inequality  (5),  first  multiply  through  by  T  and  then 
integrate  over  the  local  volume  V L  and  divide  by  this  quantity  to 
obtain 


}  j's  TdV  -  }  Jr  dV  +  J  Jt(  qj/T)fjdV  _>  0 


.  (  A 1  6  ) 


Now  define 


TL  =  ijTdV 


,  (  A 1  7  ) 


Thus,  for  the  case  when  T  is  a  linear  function  of  coordinates  in 
V k ,  definitions  (A10)  and  (A21)  may  be  substituted  into  (A20)  and 
this  result  into  ( A 1 9 )  to  obtain 


TL+  (“t 


V) 


j  >  sc  >  0 


,  (A22) 


where 


3C  E 


/T)dV  - 


(1/T lVl) 


/q JdV  :  /8J 


(A23) 


sc  can  be  shown  to  be  strictly  nonnegative  with  the  assumption 
that  T  is  everywhere  nonnegative,  along  with  equation  (9). 

We  now  assume  that  the  local  volume  is  small  enough 
compared  to  B  that  the  standard  procedure  may  be  utilized  to 
obtain  the  linear  conservation  of  momentum  equations  [39] 


°Lji.j 


(  A2  4) 


similar  to  pointwise  equations  (2),  and  the  conservation  of 
angular  momentum  may  also  be  obtained 


CTLij 


,  (A25) 


similar  to  equations  (3)*  Thus,  it  is  assumed  that  no  body 
moments  are  introduced  via  material  inhomogeneity  or  other 
sources.  This  assumption  must  be  relaxed  when  the  model  is 
utilized  for  interply  delamination,  since  in  this  case  the  local 
volume  element  goes  through  the  entire  laminate  thickness. 

Equations  (A24),  (A25),  (A15),  (A22),  ( A1 4) ,  (A18),  (A5), 
(A7),  and  (A6)  are  rewritten  as  equations  (12)  through  (20), 
respectively,  in  the  main  text. 


rv 


"VW 


unci assi f ied 

Sf  Cl’Ritv  CLASSIFICATION  OF  This  PAGE 


la  HEPOHT  SUOHITV  L  L  ASS)  F  i  CAT  JON 

Unclassified 


2a  SECURITY  CLASS  IFICATlGN  AUThOK.IV 


2b.  DE  CLASS  I  F  ICATION,  DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBERiS) 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


3.  DiSTRI  BUT  ION/A  VA:  LAB  I  LIT  Y  OF  REPORT 

Uni imited 


5.  MON  IT  OR  I  NO  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Aerospace  Engineering  Dept. 


6c.  AOORES5  (City,  State  and  ZIP  Code ) 

Texas  A&M  University 
College  Station,  TX  77843 


S«  NAME  OF  FUNOING/SPONSORING 
ORGANIZATION 

AFOSR 


Be  ADDRESS  (City,  State  and  ZIP  Code ) 


b.  OFFICE  SYMBOL  7#.  NAME  OF  MONITORING  ORGANIZATION 
(If  applicable ) 

Air  Force  Office  of  Scientific  Research 


7b.  ADORESS  (City,  State  and  ZIP  Code) 

Bol ling  AFB 

Washington,  D.C.  20332 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable ) 

Grant  No.  AF0SR-84-0067 


10  SOURCE  OF  FUNDING  NOS. 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

same 

E  LE  ME  NT  NO. 

NO. 

NO. 

NO. 

11  TITLE  ilnciude  security  Classification i  i 

A  Thermomechanical  Constitutive  Theory  for _ 

12.  PERSONAL  AUTHORIS) 

D.H.  Allen,  C.E.  Harris,  and  S.E. Groves 


13a  type  of  REPORT  13b.  TIME  COVERED 

Interim  I  FROM  TO 


14.  DATE  OF  REPORT  (Yr  ,  Mo..  Day )  15.  PAGE  COUNT 

October  1985  38 


COSAT l  CODES 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number ) 

composites  internal  state  variables 

damage  continuum  mechanics 

laminate  analysis 


19  ABSTRACT  i Continue  on  reverse  if  necessary  and  identify  by  block  number / 

A  continuum  mechanics  approach  is  utilized  herein  to  develop  a  model  for  predicting 
the  thermomechanical  constitution  of  initially  elastic  composites  subjected  to  both  mono¬ 
tonic  and  cyclic  fatigue  loading.  In  this  model  the  damage  is  characterized  by  a  set  of 
second  order  tensor  valued  internal  state  variables  representing  locally  averaged  measures 
of  specific  damage  states  such  as  matrix  cracks,  fiber-matrix  debonding,  interlaminar 
cracking,  or  any  other  damage  state.  Locally  averaged  history  dependent  constitutive 
equations  are  constructed  utilizing  constraints  imposed  from  thermodynamics  with  internal 
state  variables.  In  Part  I  the  thermodynamics  with  internal  state  variables  was  constructe 
and  it  was  shown  that  suitable  definitions  of  the  locally  averaged  field  variables  led  to 
useful  thermodynamic  constraints  on  a  local  scale  containing  statistically  homogeneous 
damage.  Based  on  this  result  the  Helmholtz  free  energy  was  then  expanded  in  a  Taylor  serie 
in  terms  of  strain,  temperature,  and  the  internal  state  variables  to  obtain  the  stress-stra 
relation  for  composites  with  damage.  In  Part  II,  the  three  dimensional  tensor  equations  n 
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PART  II:  Application  to  Matrix  Cracking 
i n  Lam i nated  Compos i tes 

by 

D . H .  Allen 
C.E.  Harris 
S . E .  Groves 


ABSTRACT 

A  continuum  mechanics  approach  is  utilized  herein  to 
develop  a  model  for  predicting  the  thermomechanical  constitution 
of  initially  elastic  composites  subjected  to  both  monotonic  and 
cyclic  fatigue  loading.  In  this  model  the  damage  is 
characterized  by  a  set  of  second  order  tensor  valued  internal 
state  variables  representing  local ly  averaged  measures  of 
specific  damage  states  such  as  matrix  cracks,  fiber-matrix 
debonding,  interlaminar  cracking,  or  any  other  damage  state. 
Locally  averaged  history  dependent  constitutive  equations  are 
constructed  ut i 1 i z i ng  constraints  imposed  from  thermodynamics 
with  internal  state  variables.  In  Part  I  the  thermodynamics  with 
internal  state  variables  was  constructed  and  ft  was  shown  that 
suitable  definitions  of  the  locally  averaged  field  variables  led 
to  useful  thermodynamic  constraints  on  a  local  scale  containing 
statistical ly  homogeneous  damage.  Based  on  .this  result  the 
Helmholtz  free  energy  was  then  expanded  in  a  Taylor  series  in 
terms  of  strain,  temperature,  and  the  internal  state  variables  to 
obtain  the  stress-strain  relation  for  composites  with  damage. 
In  Part  II,  the  three  dimensional  tensor  equations  from  Part  I 
[1]  are  simplified  using  symmetry  constraints.  After  introducing 
engineering  notation  and  expressing  the  constitutive  equations  in 
the  standard  laminate  coordinate  system,  a  special  i  zed 
constitutive  model  is  developed  for  the  case  of  matrix  cracks 
only.  The  potential  of  the  model  to  predict  degradation  of 
effective  stiffness  components  is  demonstrated  by  solving  the 
problem  of  transverse  matrix  cracks  in  the  90°  layer  of  several 
crossply  laminates. 

To  solve  the  example  problems,  the  undamaged  moduli  are 
determined  from  experimental  data.  The  internal  state  variable 
for  matrix  cracking  is  then  related  to  the  strain  energy  release 
rate  due  to  cracking  by  utilizing  li near  elastic  fracture 
mechanics.  These  values  are  then  utilized  as  input  to  a  modified 
laminate  ana  lysis  scheme  to  predict  effective  stiffnesses  in  a 
variety  of  crossply  laminates.  The  values  of  effective  (damage 
degraded)  stiffnesses  predicted  by  the  constitutive  model  are  in 
agreement  with  experimental  results.  The  agreement  obtained  in 
these  example  problems,  while  limited  to  transverse  matrix  cracks 
only,  demonstrates  the  potential  of  the  constitutive  model  to 
predict  degraded  stiffnesses. 


1 


INTRODUCTION 


In  Part  I  it  was  hypothesized  that  damage  can  be  modeled  in 
continuous  fiber  composites  by  a  set  of  second  order  tensor 
valued  internal  state  variables  which  represent  locally  averaged 
measures  of  matrix  cracking,  interlaminar  de 1  am i nat ion,  and  other 
damage  mechanisms  on  a  scale  which  is  assumed  to  be  smal I 
compared  to  the  scale  of  the  boundary  value  problem  of  interest. 
Continuum  mechanics  with  internal  state  variables  [I]  was  then 
utilized  to  construct  stress-strain  relations  in  which  all 
components  of  the  degraded  modulus  tensor  can  be  determined  for  a 
given  damage  state. 

The  purpose  of  this  paper  (Part  II)  is  to  demonstrate  how 
the  model  may  be  utilized  to  predict  the  reduced  stiffness  of 
laminates  which  are  subjected  to  known  damage  states.  This 
procedure  is  illustrated  via  specific  examples  in  which  there  is 
a  single  damage  mode  consisting  of  matrix  cracking.  It  is  shown 
that  the  properties  of  a  single  lamina  with  known  damage  can  be 
utilized  as  given  properties  to  obtain  favorable  comparisons  to 
experimental  results  for  specific  laminates. 

The  model  application  is  accomplished  by  first  imposing 
symmetry  constraints,  reducing  to  generalized  plane  stress,  and 
finally  performing  the  laminate  integration. 


SIMPLIFICATION  OF  THE  MODEL 


We  now  consider  the  stress-strain  relation  described  in 
equations  (42)  through  (45)  of  Part  I  (see  Appendix  A).  For  the 
examples  to  be  considered  herein,  it  is  assumed  that  all  residual 
stress  components  are  zero  ( aLi  j  =  and  that  there  are  no 
temperature  changes  (AT^=0). 


Reduct i on  to  S i ng 1 e- I ndex  Not at i on 

By  incorporating  the  symmetry  of  the  stress  and  strain 
tensors,  the  quadratic  dependence  of  the  Helmholtz  free  energy  on 
strain,  and  the  Voigt  single  index  notation  [2],  the  constitutive 
equations  reduce  to  (see  Appendix  A) 


Although  we  have  dropped  the  subscript  L,  all  quantities  in 
equations  (1)  represent  locally  averaged  measures.  The 
subscripts  i  and  j  range  from  1  to  6,  the  subscript  k  ranges  from 
1  to  9,  and  the  superscript  n  ranges  from  1  to  N,  the  number  of 
damage  modes. 

At  this  point,  further  reductions  can  be  made  to  the  number 
of  unknown  constants  in  equations  (1)  only  by  specifying  the 
material  symmetry  and  specific  damage  modes  of  interest. 


'Vv 


Mater i a  1  Symmetry  Constra  i  nts 


Material  symmetries  may  now  be  utilized  to  further  simplify 
the  constitutive  equations.  The  material  in  question  is  assumed 
to  be  initially  transversely  isotropic  in  the  undamaged  state  on 
the  local  scale,  where  the  plane  of  isotropy  is  the  x?  -  x,  plane 


shown  in  Fig.  2 
is  given  by  [3] 
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where  C44=2(S22“S23)»  and  S22  and  S23  are  the  associated 
compl fance  terms. 

It  is  assumed  that  the  crack  induces  orthotropy  in  three 
planes;  the  plane  of  the  crack,  the  plane  in  which  the  crack 
opening  displacement  uc  and  crack  normal  nc  lie,  and  a  third 
plane  which  is  orthogonal  to  the  first  two.  Therefore,  the 
damage  tensor  ljk  is  an  orthotropic  tensor  containing  15  unknown 
constants  in  the  coordinates  described  by  the  crack  geometry  (see 
Appendix  B),  given  by 
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Thus,  the  complete  constitutive  equations  (1) 
(assuming  the  damage  growth  law  is  known)  require  the 
determination  of  5  independent  material  constants  for  the 
undamaged  modulus  tensor  C,  jj,  and  15  independent  constants  for 
the  damage  tensor,  I  }  It  should  be  noted,  however,  that  the 
planes  of  these  symmetries  will  not  coincide  unless  the  crack 


cl  f  spl  ac^ment  u  Is  orfented  parallel  to  the  Xj,  x2»  or  axis  in 
ply  cc  di nates. 

App 1 i cat i on  to  Matr i x  Crack i no  i n  Cont i nuous  F f ber  Laminates 

As  discussed  in  the  int  -oduction,  the  capability  of  the 
constitutive  model  will  be  demonstrated  by  considering  the  case 
of  matrix  cracking  in  continuous  fiber  laminated  composites.  An 
example  of  this  damage  state  is  shown  schematically  in  Fig.  1. 
In  order  to  apply  the  proposed  constitutive  model  to  this  system 
we  first  examine  the  response  of  a  single  ply  subjected  to 
transverse  matrix  crack i ng,  as  shown  in  Fig.  2.  Assuming  that 
the  crack  geometry  is  symmetric  about  normals  to  each  of  the  ply 
coordinates,  the  internal  state  variable  associated  with  matrix 
cracking  is  represented  in  ply  coordinates  by 


[a*]  =  [  0  al2  0  0  0  0  0  ]  ,  (4) 


where  the  single  subscripted  notation  is  defined  by  equations 
(7a).  This  implies  that  the  crack  normal  nc  in  a  single  ply  is 
parallel  to  the  local  x2  coordinate.  Furthermore,  the  crack¬ 
opening  displacement,  ijc,  may  contain  three  components. 

Note  that  a  second  order  tensor  representation  of  the 
internal  state  variable  may  be  insufficient  if  the  crack 
displacement  uc  or  normal  nc  rotates  during  the  load  history.  In 
this  case  a  higher  order  tensor  may  be  required  [4].  However, 
since  the  crack  is  matrix  dominated  and  constrained  by  fibers, 
time  dependent  rotation  is  assumed  to  be  negligible  and  the 
second  order  tensor ial  representat ion  is  considered  adequate  in 
the  current  model.  Recent  experimental  evidence  [5]  indicates 
that  cracks  do  indeed  change  planes  sometimes  in  multi-ply 
laminates  with  several  adjacent  crosspl ies  at  the  same 
orientation.  However,  the  crack  plane  is  essentially  straight  in 
each  ply,  the  level  at  which  the  local  volume  is  constructed  for 
matrix  cracks. 

For  the  single  damage  mode  of  matrix  cracking  described  in 
Fig.  1,  equations  (1)  reduce  to 

°i=Ci Jc j+I f2a2+I i5a5+I I8a8  '  (5) 

For  relatively  thin  laminates  it  is  useful  to  apply  the 
conditions  of  generalized  plane  stress  where  the  out-of-plane 
shear  stresses  a4  and  a ^  are  neglected.  Applying  these 
conditions  to  equations  (5),  imposing  the  symmetry  constraints 
described  in  equations  (2)  and  (3),  and  using  matrix  notation 
results  in 


C1 1  C 1 2  C 1 2  0 


C  1 2  C22  C23  0 


C i 2  C23  C33  0 
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+  ClhCa1} 


,  (6) 


where 


[  I  ]  = 


I{2  0 


1  32  0 


I  eZQ  0 


Note  that  the  fifth  column  of  the  coefficient  matrix  in  equations 
(5)  is  zero  due  to  the  fact  that  a  5  does  not  contribute  to  the 
in-plane  stresses  in  the  generalized  plane  stress  reduction. 
F  urt  her  more,  note  that  I  !  ?  »  1  22  and  J  32  are  the  coe^^  *  c  '  ent  s  of 
the  effect  of  loss  of  stiffness  on  the  normal  stresses  o . ,  o2  and 
Oo.  respectively.  Finally,  note  that  I is  the  coefficient 


respectively,  r  mca  i  i  y  »  ritice  in  at  i  gQ  '  »  trie  t  cer  r  i  t  i  er  1 1 

determining  the  influence  of  stiffness  loss  on  the  in-plane  shear 


stress  a, 


It  is  apparent  from  the  above  equations  that  for 


generalized  plane  stress  conditions  there  are  ten  unknown 
material  constants  to  be  determined  for  the  case  of  matrix 
crack i ng . 


Determ i nat i on  of  the  I  Matr i x 


Theoretically  it  is  possible  to  determine  the  I  matrix 
directly  from  experimental  data.  This  may  be  accomplished  by 
subjecting  test  coupons  to  prescribed  deformation  histories, 
removing  the  deformations,  and  nondestruct i ve 1 y  evaluating  the 
damage  state.  The  residual  stresses  will  determine  the  I  matrix. 
However,  in  graphite/epoxy  laminates  this  procedure  breaks  down 
due  to  the  fact  that  although  the  crack  surfaces  may  be 
determined  nondestruct i ve 1 y  using  x-rays  and  edge  replicas,  the 
crack  opening  displacements  cannot  be  accurately  determined 
experimentally.  Therefore,  an  alternative  approach  is  used 
herein  to  evaluate  the  I  matrix. 

As  described  in  Appendix  C,  for  the  case  considered  in  this 
paper,  at  least  to  a  first  approximation,  it  can  be  shown  that 
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Therefore,  the  number  of  unknown  material  constants  is 
reduced  to  a  total  of  six  for  the  case  considered  herein. 

Lam i nate  Eouat i ons 


In  order  to  utilize  single  lamina  equations  to  characterize 
the  response  of  multilayered  laminates,  it  is  necessary  to 
globally  average  the  local  ply  constitutive  equations.  This  is 
accomplished  herein  by  imposing  the  Kirchhoff  hypothesis  for  thin 
plates  [5].  However,  higher  order  plate  or  shell  theories  could 
be  utilized  also.  Generalized  plane  strain  conditions  are 
imposed  rather  than  plane  strain  because  this  is  consistent  with 
the  stress  state  in  equations  (6)  (A  detailed  description  of  the 
global  averaging  is  given  in  Appendix  D).  The  resulting 
equations  are  as  follows: 


[N}=[A] {e°}+{0} 


,  (9) 


or 


{e°}=[A  1 ] ( {N}~ (D) ) 


,  (10) 


where 


Aij  5  *  <Cij>k  *k  •  (11) 

k=  1 


and 


CD]  s  l  ([1  ’  (  12) 
k=  1 

[I1]^.  and  (a1}k  are  in  laminate  coordinates  as  defined  in 
Append ix  E ,  and  (e°)  contains  the  laminate  midplane  strains. 
Furthermore,  k  specifies  the  ply  and  t  k  is  the  ply  thickness. 
For  convenience,  we  have  assumed  that  no  moments  are  produced  by 
the  damage  (in  the  absence  of  curvature),  which  is  assumed  to 
hold  for  symmetric  laminates. 

In  order  to  determine  the  effective  stiffness  for  any  damage 
state,  we  evaluate  the  rate  of  change  of  (N)  with  respect  to  the 
midsurface  strains  {e°}  during  unloading,  that  is. 
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where  Sfm 


fs  defined  to  be  the  effective  stiffness, 
work  on "graph i t e/epoxy  laminates  has  shown  that 
nearly  a  constant  during  unloading,  implying  that, 
first  approximation  for  crossply  laminates. 


Exper i menta 1 
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THE  INTERNAL  STATE  VARIABLE  FOR  MATRIX  CRACK  DAMAGE 


Equation  (41)  in  Part  I  [1]  gave  the  second  order  Taylor 
series  expansion  of  the  local  energy  per  unit  volume  due  to 


in  terms 


of  strain, 
n 


Li 


to 


temperature,  AT,  and  the 
or  demonstration  purposes, 
the  predictions  of  this  paper  are  being  confined  to  symmetric 
cross-ply  laminates  loaded  in  uniaxial  tension.  For  this  case. 


cracking, 

internal  state  variable  (ISV),  aj^.j. 


1 


aj>  is  the  only  component  of  the  ISV  of  interest  and  is  defined 


(in  the  ply  coordinate  system)  as 


=  L 

vr  J- 


,dS 


(15) 


is  the  crack-opening  displacement,  V,  is  the  local 


where  U2 

element  volume  and  S2  is  the  surface  area  of  matrix  cracks. 
Furthermore,  we  will  consider  only  the  case  of  load-up  in  the 
fixed  grip  mode  where  matrix  crack  extension  occurs  at  constant 
strain.  Therefore,  if  the  higher  order  terms  in  the  Taylor 
series  expansion  are  neglected,  the  local  energy  due  to  cracking 
reduces  to 


u£= ( A+ I  22s 


2)<x  2 


(16) 


where  A  i s 
the  fixed 


a  constant, 
grip  mode. 


release  rate,  G  ,  for  matrix 
related  to  G^  as  follows: 


Since  equation 
a2  can  be  related 


(16)  applies  to  load- up  in 
to  the  strain  energy¬ 


cracking  by  noting  that  u,  is 
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GmdS 


(  17) 


where  it  is  assumed  that  the  initial  matrix  crack  surface 


zero  and  S^(t)  is  the  surface  area  at  time  t.  If 


e  area  15 
we  make  the 

energy  stored  due  to  residual  damage 
negligible,  then  the  constant  A  in  equation 
equating  (16)  and  (17)  yields 


assumption  that  the 
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■For  stable  crack  growth.  In  order  to  properly  account  for  crack 
interaction  an  expression  for  Gm  will  be  determined 
exper i menta 1 1 y . 

Using  the  standard  definition  of  strai..  energy  release  rate 
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J_  3_U 
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(  19) 


and  noting  that  matrix  crack  surface  area  is  defined  as 


S  =  2Ba 


(20) 


the  strain  energy  release  rate  can  be  expressed  in  terms  of 
matrix  crack  surface  area  as  follows 


(21  ) 


Furthermore,  the  strain  energy  in  a  90°  layer  is  defined  by 


U  =  E22  e f 


(22) 


for  a  uni  form  applied  strain,  e  2  *  in  an  elastic  material. 
Assuming  the  applied  strain  is  constant,  substituting  (22)  into 
(21)  resu Its  in 


■  -vi  Ihi 

as 


(23) 


It  is  noted  that  Equation  (23)  is  similar  to  the  expression  for 
strain  energy  release  rate  written  in  terms  of  test  specimen 
comp  I ?  ance . 

The  rule  of  mixtures  yields  the  following  expression  for  the 
loading  direction  modulus  of  a  cross-ply  laminate: 


E  =  pEl 1  *  qE22 
P+q 


(24) 


10 


where  p  is  the  number  of  0  pi ies  and  q  is  the  number  of  90 
plies.  Assuming  that  matrix  cracks  are  confined  to  the  90° 
pi ies. 


aEx  =  *  3E; 


Substituting  (25)  into  (23)  gives 


G  —  _w  c2  P+q  3  E. 
um  "  vLe  2  — - 


If  the  right  hand  side  of  equation  (26)  is  determined 
experimental  ly  from  a  laminate  that  has  a  90°  layer  that  is  a 
single  ply  thick,  the  resulting  strain  energy  release  rate  can  be 
utilized  for  other  layups  by  observing  that  the  strain  energy 
re  1  ease  rate  for  a  p  1  y  i  n  a  1  ayer  that  is  n  plies  thick  is  g  i  ven 
by  (See  Appendix  F.) 


(Gm>I 


-n  VL 


,  (27) 


1  ply  layer 


where  VL  is  the  volume  of  a  single  ply.  Equation  (27)  can  now  be 
substituted  into  equation  (18)  to  obtain  the  expression  for  the 
I SV  of  a  single  90°  ply  for  matrix  crack  extension  during  load- 
up.  The  resulting  equation  is 


■S2(t ) 


„1  _  e,  n ( p+q ) 

2  '  7i - “ 

1 2?  d 


IQx  dS 
3  S 


.  (28) 


Tne  integral  term  in  equation  (28)  can  be  evaluated  as  follows: 


-S2(t ) 


•x  dS 


dExl  Exl  Exo 


,  (29) 


where 


is  the  undamaged  elastic  modulus  and  Exl  is  the 


degraded  modulus  corresponding  to  damage  state  S  2 ( t j  ) . 
Substituting  (29)  into  (28)  and  rearranging,  the  I  SV  for  load-up 


is  expressed  as 


.  (  30  ) 


a2 


=  £ 


1  oad-up 


at 


S2(tj  ) 


-  J 

s2(t,) 


Although  it  is  possible  for  matrix  crack  surface  area  to  increase 
during  unloading,  in  the  current  development  this  effect  is 
assumed  to  be  negligible.  Therefore,  on  unloading  ai>  depends 
only  on  the  crack-closure  displacement,  u2  in  equation  (15),  and 
would  go  to  zero  on  complete  crack  closure.  Assuming  that  the 
crack-c I osure  displacement  is  linear  with  strain  and  the  matrix 
crack  surface  area  is  constant,  equation  (15)  can  be  rewritten  as 


a 


1 

2 


=  ce 


un 1 oad i ng 


(31  ) 


where  c  is  a  constant  of  proport i ona 1 i ty. 

The  constants  in  equations  (30)  and  (31)  must  be  determined 
from  experimental  data.  Considering  a  tensile  test  with  a  load 
and  unloading  cycle,  at  the  instant  of  load  reversal  the 
expressions  for  the  I SV  for  load-up  and  unloading  must  be  equal. 
Therefore,  setting  equation  (30)  equal  to  (31)  and  rearranging, 
gives  the  following  relationship 


c 


=  n ( p+g ) 

q 


S2(t  j  ) 


(32) 


It  should  be  noted  that  all  matrix  cracking  information  is 
contained  in  the  term  Ex|/Exo-  Since  equation  (32)  applies  only 
to  a  single  90°  ply,  Ex 1  was  determined  from  the  experimental 
results  of  the  [0/90/0Js  laminate.  The  following  expression  was 
obtained  from  a  least  squares  curve  fit  to  the  experimental 
values  of  Ex|/Exo  versus  S^: 


Exl/Exo  =  0-^969-0.06  1  607*S2(t  j  )  +  0.046230-S2(t  j  )2  .  (33) 


Reca  1  1  i  ng  equat  i  on  (13),  it  is  seen  that  the  effect  i  ve 
stiffness  of  a  laminate  can  be  obtained  by  specifying  3a2/3e2. 
On  unloading  this  is  given  by  equation  (31)  to  be 


3 g  2 / ^ £  2 


(34) 


Therefore,  equations  (32)  and  (33)  are  used  with  laminate 
equations  (13)  to  predict  the  effective  stiffness  of  any  laminate 


MODEL  COMPARISON  TO  EXPERIMENT 


The  model  has  been  utilized  to  predict  the  damage  dependent 
reduced  stiffness  of  several  crossply  laminates.  This  has  been 
accomplished  by  utilizing  the  laminate  stiffness  equations  (13), 
in  conjunction  with  the  damage  evaluation  procedure  described  in 
the  previous  section.  The  reduced  s;iffnesses  predicted  by  the 
model  have  been  compared  to  experimental  results  obtained  from 
graphite/epoxy  coupons  composed  of  Hercules  AS4/3502. 

The  coupons  were  obtained  from  laminates  fabricated  from 
prepreg  tape  using  a  hot  press  technique  in  the  Mechanics  and 
Materials  Center  at  Texas  A&M  University.  Quasi-static  tensi le 
tests  were  conducted  on  an  Instron  1128  screw  driven  uniaxial 
testing  machine.  Matrix  crack  damage  states  were  evaluated  by  x- 
ray  radiography  and  edge  replication.  Further  details  of  these 
procedures  are  contained  in  reference  [5].  Initial  undamaged 
1  am i na  properties  are  given  in  Table  1.  These  properties  were 
obtained  exper i menta 1 1 y  from  [0]g,  [90]g  and  ±  C  4  5 ] ^  s  laminates. 
Energy  release  rates  were  obtained  from  [0*90,0]  control  coupons. 
Experimental  results  for  stiffness  loss  as  a  function  of  surface 
area  of  matrix  cracks  are  given  in  Fig.  3. 

The  I  am i na  properties  and  strain  energy  release  rate 
described  above  were  used  as  input  to  the  model  to  produce  the 
predicted  results  for  the  [0,90]  ,  [0,902]s*  [0,90,0]s,  and 
[ 0 , 9 0 ^ ] s  laminates  shown  in  Figs.  4  through  7.  As  shown  in  the 
figures,  the  comparisons  to  the  exper i menta 1 1 y  determined  values 
of  reduced  stiffness  are  quite  good  in  all  cases. 


SUMMARY  AND  CONCLUSIONS 

A  model  for  predicting  the  stiffness  loss  in  laminated 
composites  as  a  function  of  microstructural  damage  has  been 
proposed  in  this  two  part  paper.  In  part  I  the  general 
theoretical  framework  was  constructed  for  elastic  composites  with 
damage.  In  part  I  I  the  model  has  been  special  i  zed  for  the  case 
of  matrix  cracks  in  crossply  laminates.  In  this  process  the 
following  key  developments  have  been  reported: 

1)  material  symmetry  constraints  have  been  imposed  on  the 
damage  constant  tensor  I-j^; 

2)  the  damage  tensor  oj j  has  been  reduced  for  the  case  of 
plane  stress; 

3)  an  approximate  procedure  has  been  proposed  for  obta i n i ng 
the  damage  constant  tensor; 

4)  damage  dependent  laminate  equations  have  been 
constructed;  and 

5)  the  internal  state  for  any  crossply  layup  has  been  found 
to  be  derivable  from  energy  release  rates  experimentally  obtained 
from  a  single  layup. 

The  model  has  been  demonstrated  to  be  accurate  in  predicting 
the  damage  dependent  reduced  stiffness  of  several  graph i te/epoxy 
crossply  laminates  with  matrix  cracks. 


Current  and  future  development  of  the  model  will  deal  with 
the  following  issues: 

1  )  app li cation  of  the  mode  1  to  laminates  with  both  matrix 
cracks  and  interply  de 1  ami  net  ions ; 

2)  application  of  the  model  to  layups  more  complex  than 
crossply  laminates;  and 

3)  development  of  internal  state  variable  growth  laws  for 
matrix  cracking  and  interply  del  am? nation. 
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Tab 1 e  1  . 


Material  Properties  for  Hercules  AS4/3502 


Lami na  Propert i es 


Longitudinal  Modulus 

E11 

21.0  x  1 06  ±  2.07. 

PS  i 

Transverse  Modulus 

E22 

1 .39  x  106  ±  2.17. 

ps  i 

Shear  Modulus 

G1  2 

0.694  x  106  psi 

Poisson's  Ratio 

v  12 

0.310  ±  3.77. 

Longitudinal  Strength 

Ftu  1 

326,000  ±  3.57.  ps  ? 

Transverse  Strength 

Ftu2 

1  1 , 085  ±  9.87.  ps  i 

Long.  Failure  Strain 

etul 

0.0144  ±  4.67.  in/ in 

Tran.  Failure  Strain 

etu2 

0.00773  ±  6.77.  in/ in 

! 
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Vv 


Fig.  4.  Model  vs.  Experiment 
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APPENOIX  A:  APPLICATION  OF  SYMMETRY  CONSTRAINTS 


The  damage-dependent  constitutive  model  (equations  (42) 
through  (45)  of  Part  I  [1])  is  defined  as  follows: 


=a  U+C 


i  jkl  (ckl  -£k  1  )  +  lfi  jklakl 


(  la) 


where  a j  ,  Is  the  local  stress  tensor,  ek)  Is  the  local  strain 
tensor,  ®ij  *3  the  residual  stress  in  the  absence  of  strain  and 
temperature  change,  Cjjk1  is  the  undamaged  modulus  tensor,  is 

the  thermal  strain  tensor,  aj^j  is  the  internal  state  variable 
tensor,  and  I^j^i  is  the  damage  modulus  tensor.  Furthermore,  we 
have  dropped  Jthe  subscript  L  (denoting  locally  averaged 
quantities)  for  convenience. 

For  demonstrat ion  purposes,  the  residual  stress  tensor  and 
the  temperature  change  are  assumed  to  be  negligible,  resulting  in 


the  local  stress  tensor,  ek)  is  the  local  strain 


is  the  residual  stress  in  the  absence  of  strain 
lange,  Cjjkj  is  the  undamaged  modulus  tensor, 


is  the  i nterna 1 


e  of  strain  and 
s  tensor,  is 
state  variable 
Furthermore,  we 


<,ij=Cijklekl  +  llijklakl 


(2a) 


Note  that  Ifjk]  is  a  fourth  order  tensor  with  81  coefficients  for 
each  value  of  n-  It  is  assumed  here  that  the  constitutive 
equations  given  by  (2a)  are  statistically  homogeneous. 
Therefore,  the  conditions  of  stress  and  strain  symmetry  as  well 
as  the  existence  of  an  elastic  potential  can  be  applied  to 
equations  (2a)  to  obtain 


C i jk 1 =C j i k 1 • 


C i jk 1 =C  f  j 1 k ’ 


C i jk l =Ck 1 i j 


*  (3a) 


■■'In 


( 4a) 


It  is  most  convenient  at  this  point  to  refndex  the  constitutive 
tensors  using  the  Voigt  notation  [2]  where 


Q 

III 

Q 

0 45  °  23=°  32 

°2~° 22 

°5S°  '  3=a 3  1 

,  (5a) 

Q 

u> 

hi 

Q 

U) 

OJ 

°650 12=°21 
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22 


Furthermore,  for  all  values  of  n 


0 

III 

a 

a4Sa23 

0 

in 

o 

LU 

a2  5a  22 

a  5Sa  32 

OqSo  J2 

0) 

a3  5a  33 

a6Sa  13 

a9=a21 

Using  the  contracted  notation,  equations  (2a)  can  be  written  as 

ai=CiJe j+I?kak  ’  (8a) 


where  i  and  j  range  from  1  to  6,  k  ranges  from  1  to  9,  and  n 
ranges  from  1  to  N,  where  N  i  s  the  number  of  damage  modes. 


^ *  IT-  -C* 


APPENDIX  B:  SYMMETRY  CONSTRAINTS  ON  THE  DAMAGE  MODULUS  TENSOR 


Consider  the  following  component  of  internal  energy  due  to 
crack i ng : 

u?  5  1 i jkle i jak 1  *  ( lb) 

Since  the  strain  tensor  is  symmetric 


I 


1  -I  1 
i  J  k  1  “  1  j  i  k  1 


.  (2b) 


Therefore,  there  are  54  independent  constants  in  the  damage 
modulus  tensor  I | j k 1  *  Expanding  out  equation  (lb)  thus  gives 


.  ,c-  t  1  «•  „  1 

ul_1  1  1  1  le  1  la  1  1 
+I 1 132e 1 la32+I 

+  1  },2i«i !a2 l  +  I 

+I 2223e22a23+I 

+  I  22 1 2e  22a I 2+ 1 

+,3333e33a33+I 
+ 1 333 1 e  33a  3 1 +  1 
+  I 2322e 23a  22+ ^ 
+I2313e23a 13+I 
+  1  1 3 1  1 e 1 3a  1  1  +  1 
+  1  l332e i 3a  32+ 1 
+  1  } 32 i e 1 3a  2 1  +  1 
+1 i223e 12a23+I 
+  I  1 2  1  2C  12a 1 2+ 1 


+I 1 122e 1 la22+I 1 !33e 1 la33 
1 1 1 3e 1 la 13+I 1 1 3 1 c 1 la31+1 
221 1 e  22° 1 1+I2222e22a22+1 
2232£22a32+I2213e22a 13+I 
2221e22°21+I331 Ie33al 1+I 
3323e33a23+I 3332e33a32+I 
3312e33a 12+I3321e33a2I+I 
2333e23a33+I2323e23a23+I 
2331e23a31+12312e23a 12+I 
1 322e 1 3a  22+ 1  1 333e 13a33+I 
1 3 1 3e 1 3a 1 3  + 1  1  33  1  e  1 3a  3 1 +  1 
121 le 12a 1 1+I 1 222e 12a22+I 
1 232e 12a32+I  1  2  1  3e  1 2a 1 3+ 1 

1  r  «  1 

1 22 1 e 1 2a  2 1 


.  ,  1  P  „  1 
1 123e 1 la23 

1  P  „  1 

1  1  1  2e  1  1 a  1  2 

1  p  „  1 

2233e  22a33 

1  p  Q  1 

223 1 £  22a3 1 

1  e  a  1 

3322e33a22 

1  E  Q  1 

33 1 3e  33a  1  3 

1  p  „  1 

23 1  1 6  23a  1  1 

1  p  „  1 

2332e23a32 

1  p  „  1 

232 1 e  23a2  1 

1  p  „  1 

1 323e 1 3a23 

1  p  a  1 

13 12e  13a  12 

1  p  „  1 

1 233e 1  2a  33 

1  p  „  1 

123 le 1 2a  3  1 


.  (3b) 


We  now  wish  to  impose  orthotropic  symmetry.  In  order  to  do 
this,  first  rotate  180°  about  the  x^  axis  [2],  The  direction 
cosines  for  this  transformat i on  are 


[aik']  = 


-1  0  0 

0-1  0 
0  0  l 


.  (4b) 


Therefore,  since  ejj  is  a  second  order  tensor, 

ek' 1 '=ei jaik'ajl ' 
it  follows  that 


,  (5b) 


Furthermore , 


1 1 

c  12 

_e  13 

12 

e  22 

_e  23 

13 

-e  23 

e  33 

,  (6b) 


all 

a  1  2 

a  1  3 

a2  1 

q22 

-a23 

1 

a3  1 

1 

a  32 

„  1 
a  33 

(7b) 


Since  u?  must  be  independent  of  coordinate  system 


■  ,c_  r  1  ,  ^  1 

ul-‘p'q'r's'ep'q,Qr's' 


Substituting  (6b)  and  (7b)  into  (8b)  and  comparing  this  result  to 
(3b)  will  result  in 


1  1 1 23” 1  1  1  32"  1  1  1  13" 1  1  131  =  12223"l2232“u 

I  1  _  i  1  _  I  1  _  |  1  _  i  1  _  i  1  _o 

*2213" *223  1" *3323" *3332" *33 13" *333 l"u 

I  1  _  T  1  _  I  1  _  I  1  _  |  1  _  I  1  _0 

*231 1" *2322" *2333" *23 12"* 232 1"* 131 l"u 

I  1  _  i  1  _jl  _  r I  _  r  1  _ |  1  _  n 

1 1322" *  1333" *  131 2" 1 1321" 1 1223" 1 1232"u 
1  1213=I 1231=0 

Rotating  180°  about  the  axis  gives 


.  (9b) 
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Therefore 


C  ek/  ]  /  1 


Furthermore , 
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al  1 
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1 
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1 

„  1 

a3  1 

-a  32 
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.  (lib 


(12b) 


Substituting  (11b)  and  (12b)  into  (8b)  and  comparing  to  (3b) 
results  in 


I  1  _jl  _ 

1 1 1 12~J 1 121_12212_ 

I  1  _  j  1  _  i  1  _ 

23 1 3“ 1 233 1 ~ 1 1323~ 


i  1  _ 

*222  1~ 

i  1  _ 

1 1332" 


i  1  _  ]  1  _  g 

13312“J3321_U 

I  1  _  t  1  _  t  1  _n 

1  121 l-1 1 222~ 1 1 233-u 


.  (13b) 


Rotating  180°  about  the  Xj  axis  yields  no  additional  constraints. 
Therefore,  imposition  of  orthotropic  symmetry  on  ljjk,  reduces 
the  number  of  constants  to  15.  These  are 


Therefore,  the  orthotropfc  damage  modulus  matrix  is  given  by 


Cl1]  = 


1  1  1  1  1  2  1  1  3  0 
1 2 1  *22  1 23  0 
*31  1 32  1 33  0 


0  0 


1 44  1 45  0 


I  1  !  1 

1 56  57 


I  1  i  1 

68  l69\ 


(15b) 


For  the  case  where 


[a1]  =  [0 


equation  (15b)  reduces  to 


,  (16b) 


I  ,  0 


loo  0 


0  0 


0  0 


Cl1]  = 


I  -3Q  0 


0  0 


0  0 


Ii*  0 


.  (17b) 


0  0 


0  0 


0  0 


0  0 


APPENDIX  C:  DETERMINATION  OF  THE  I  MATRIX 


At  a  material  point  In  V|_  the  stress-strain  relation  in  the 
absence  of  temperature  change  i s 


a  i  J  "  C1 jkl ekl 


.  (  lc) 


Integrating  over  the  local  volume  (excluding  cracks)  gives 


,  (2c) 


where  ^  =  is  the  average  stress  outside  the  damage  zones.  In 
this  section  this  value  of  the  stress  tensor  is  assumed  to  be 
identical  to  the  average  stress  a^jj,  which  includes  the  average 
of  the  stress  in  the  damage  zones.  Assuming  that  C.  jkJ  is 
spatial ly  homogeneous  in  VL,  the  above  may  be  written  J 


°Li  j 


jm  JekidV  -  fr-uni Jt<u k. 


l+ul,k>dV 


VL-VC 


VL"VC 


Using  the  divergence  theorem  on  the  last  term  gives 


(3c) 


L  f  j 


=  C 
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[J—  ti-(ukn 

VL  J2  k 


, +u,nk)dS+ 


1+u1nk)dS] 


.  (4c) 


Or ,  equ i va 1 ent l y , 


'Li  j 


=  Ci jkl (eLkl 


kl'2* 


1  k  ■ 


.  (5c) 


However,  the  Taylor  series  expansion  has  already  given  (for 
isothermal  conditions) 


°lli  j  =  CI  Jkl6kl  +  I  I  Jkl°kl 


.  (6c) 


Therefore,  equating  like  terms  in  equations  (5c)  and  (6c)  gives 


APPENOIX  0:  LAMINATE  EQUATIONS 


The  values  of  generalized  plane  strain  are  given  by 


/•* ' 

OX 

/<x  A 

,o 

cy 

'=  1 

ey 

+z 

r 

y  1 

e_ 

0 

z 

z 

1 

t  CXy J 

X  0 
< 

1 

Kxy  ' 

where  the  superscript  o  denotes  the  mfdsurface  strains  and  the 
k  matrix  denotes  the  midsurface  curvatures.  Under  the  condition 
of  generalized  plane  strain  there  is  no  warping  allowed  out-of¬ 
plane,  which  implies  that  *z=0. 

It  is  now  assumed  that  no  moments  or  curvatures  are  imposed 
and  that  all  laminates  studied  are  symmetric  through  the 
thickness  (including  damage).  Therefore,  in  order  to  determine 
the  resultant  forces,  it  is  necessary  only  to  integrate  the  given 
stress  state  over  the  laminate  thickness  to  obtain 


(2d) 


where  t  is  the  total  thickness  of  the  laminate. 

Substituting  equations  (8a)  and  (Id)  into  (2d)  for  the  case 
where  there  are  no  rotations  results  in 


{N} 


/ 


t/2 


-t/2 


([C](e°}  +  [I * ] {a  1 }dz 


,  (3d) 


where  {N}  denotes  the  force  resultants,  overbars  denote  that 
these  quantities  are  transformed  to  global  coordinates,  and  {e°} 
represents  the  mid-surface  strains.  Note  that  sfnce  transverse 
cracks  are  assumed  to  go  completely  through  the  thickness  of  the 
cracked  pi ies  the  stiffness  and  damage  are  assumed  to  be 
spacial ly  constant  through  the  thickness  of  a  single  ply. 
Therefore,  equation  (3d)  can  be  written  as 


(N)  =  I  (CC]k  (zk-zk_j ){e°}  +  [ I  1 ]k(zk-zk_, ){a‘}k) 
k=  l 


,  ( 4d) 


where  k  specifies  the  ply  and  zk~zk_ j  is  the  thickness  of  each 
ply.  One  can  define 


i j  5  E  (Ci j}k  (zk_zk-l } 
k=  1 


,  (5d) 


{D}  = 


n 

=  l 
k=  1 


<CI1]k(zk-zk_i)(a1}k) 


,  (6d) 


where 


represents  the  laminate  averaged  stiffness  matrix  and 


{D}  is  the  laminate  averaged  damage  term, 
averaged  constitutive  equations  become 


Thus,  the  laminate 


{N)=[AJ{e°}+{D} 


Experimental  testing  is  often  conducted  on  uniaxial  testing 
machines  in  which  the  applied  force  resultants  are  input  and  the 
strains  are  experimentally  determined  output.  Therefore,  at 
times,  it  is  more  convenient  to  express  the  strains  in  terms  of 
the  applied  force  resultants  as  follows: 


{e°}=[A]-1({N}-{D}) 


(8d ) 


Note  also  that  moments  wi  1  1  be  produced  even  in  the  absence  of 
strain  if  the  damage  state  is  not  symmetric  through  the 
thickness.  However,  for  the  case  considered  herein,  it  will  be 
assumed  that  all  damage  states  are  symmetric,  and  moments  are 
therefore  not  considered. 


APPENDIX  E:  TRANSFORMATION  EQUATIONS  FOR  THE  DAMAGE 
TENSOR  AND  THE  DAMAGE  MODULUS  TENSOR 


Consider  a  coordinate  rotation  0  in  the  laminate 
*2  plane)  measured  clockwise  From  the  ply  coordinate 
the  laminate  coordinate  system.  For  this  case  the 
cosines  are 


[aik']  = 


cos0  -sin6  0 
sfn0  cos0  0 
0  0  I 


Reca 1 1  that  since  a ! •  is  a  second  order  tensor 


i  J 


ak '  1  '  =  ai jaik'ajl ' 


Substituting  (le)  into  (2e)  thus  gives 

/ 


{a1}  = 


a  j  2COS0S i n0+a  22s ' 

+ 

0 

1  1  7 

-a j  2s i n0cose+a22cos  9 

+ 

0 

0 

+  0 

+ 

0 

0 

+  0 

+ 

0 

0 

+  0 

+a 

1 
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COS0 

0 

+  0 

+ 

0 

0 

+  0 

+a 

1 

32 

s  i  n0 

a } 2cos20 

+ct22s  i n0cos0  + 

0 

-a { 2S i n  0 

+a  22s ’ n0cos0  + 

0 

1  is 

p'q'r's'  1S 

given  by 

1 p ' q ' r ' s '  =  1 i jklaip'ajq'akr'als' 


p 1 ane  (x j - 
system  to 
d i r ect i on 


.  (  le) 


.  (2e) 


.  ( 3e ) 


.  ( 4e ) 


Substituting  the  nonzero  components  from  equation  (17b)  into  (4e) 


APPENDIX  F:  DIMENSIONAL  ANALYSIS  OF  STRAIN 
ENERGY  RELEASE  RATES 


Consider  the  cracked  90°  layer  shown  in  Fig.  FI.  The  total 
strain  energy  in  the  region  surrounding  the  crack  that  includes 
the  strain  energy  available  to  be  released  during  crack  extension 
is  g i ven  by 


U 


UA  +  Uc 


,  (If) 


where  UA  is  the  strain  energy  ahead  of  the  advancing  crack,  and 
Ug  is  the  strain  energy  behind  the  advancing  crack. 

In  terms  of  strain  energy  density,  U0,  equation  (If)  becomes 

U  =  (Uq)a  [ t (w-a K E ]  +  (U0)B[tatE]  ,  (2 f) 


where  i.  g  is  the  effective  length  of  the  material  from  which 
strain  energy  will  be  released  by  the  advancing  crack.  The 
effective  length  is  a  function  of  both  the  crack  spacing,  S,  and 
the  thickness  of  the  90°  layer.  Therefore,  l  g  can  be  expressed 
as 


*E  =  ct 


,  ( 3f ) 


where  c  is  a  nond i mens i ona 1  function  of  the  crack  spacing  and  the 
90°  layer  thickness.  In  addition,  the  strain  energy  density 
ahead  of  the  crack  and  behind  the  crack  can  be  expressed  as 
functions  of  the  strain  energy  density  in  the  90°  layer  in  the 
absence  of  cracks  multiplied  by  some  dimensionless  constant  that 
depends  on  the  existing  matrix  crack  spacing  and  the  thickness  of 
the  90°  layer.  Written  symbolically,  then 


and 


<Va 


Uo  fA 


(lVb 


Uo  f  B 


,  (4f ) 


,  ( 5f ) 


where  U0  is  the  strain  energy  density  in  the  90°  layer  in  the 
absence  of  cracks,  and  f A  and  fg  are  functions  of  the  crack 
spacing,  S,  and  time,  t.  Substituting  (3f),  (4f)  and  (5f)  into 
(2f)  yields 
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N*. V} 

.Vi 


■.-XV.t-KtVAV.V.L-.-.-, 


v-y 


u 


=  U0fA[ct2(w-a) ]+U0fB[ct2a] 

=  UDct2[fAW  +  a(fB-fA)]  .  (6f) 


Strain  energy  release  rate  at  each  crack  tip  is  defined  as 
fo 1 1 ows 


1  3  U 

2t  3a 


.  (7f) 


Substituting  (6f)  into  (7f)  thus  gives 

G  =  J-UQct  •  (8f) 


Notice  that  c,  fA  and  fg  are  all  functions  of  the  layer  thickness 
and  existing  matrix  crack  spacing.  Therefore,  define  a  new 
function  f  such  that 


f(S.t)  =  c(fA-fB)  .  (9f ) 

Since  the  function  f(S,t)  is  a  dimensionless  function,  it  must  be 
a  function  of  S/t.  Therefore, 


c(fA~fB)=f (S/t)  •  (10f> 

Substituting  (lOf)  into  (8f)  yields  the  following  expression  for 
the  available  strain  energy  release  rate  for  matrix  cracks 


G  =  -|t[UQf(S/t)] 


•  (Ilf) 


Therefore,  the  available  strain  energy  release  rate  due  to  matrix 
cracks  in  a  90°  layer  is  linearly  proportional  to  the  thickness 
of  the  90°  layer.  The  quantity  in  the  brackets  in  equation  (Ilf) 
is  related  to  the  properties  of  the  composite  material  system  and 
the  laminate  stacking  sequence.  This  quantity  can  be  determined 
from  experimental  data. 

If  it  is  assumed  that  cracking  occurs  when  the  available 
strain  energy  release  rate  is  equal  to  the  critical  strain  energy 
release  rate  which  is  constant  for  all  matrix  cracking,  then 

G  =  GcR  .  (12f) 


Furthermore,  for  a  linear  elastic  material  with  rigid  fibers,  the 


strain  energy  density  is  oiven  by 


_Le  e  ^ 
d??  e  22 


(  1  3f ) 


Substituting  (Ilf)  and  (13f)  into  (12f)  and  solving  for  strain 
results  in  the  following  expression  for  the  strain  in  the  90° 
layer  at  which  matrix  crack  extension  occurs: 


e22 


1/2 


tE22f (S/t ) 


{  14f ) 


Rearranging  gives 


4GcR 


tE22  c  22 


=  ffS/t) 


(  1  5f ) 


Finally,  the  thickness  of  the  90°  layer  is  given  by 


t  =  ntj  ,  ( 16f ) 

where  tj  is  the  thickness  of  one  ply  and  n  is  the  number  of 
consecutive  90°  plies.  Substituting  equation  (16f)  into  (15f) 
resu Its  in 


4G 


cR 


t  1e22  \?zh\ 


f (S/t) 


(  1  7f ) 


If  the  influence  function,  f(S/t),  is  constant  for  al 1  laminate 
stacking  sequences,  then  the  left  hand  side  of  equation  (15f) 
will  be  a  function  of  matrix  crack  surface  area  only. 

The  terms  in  parentheses  on  the  left  and  right  hand  sides  of 
equation  ( 1 7 f )  are  laminate  specific  while  all  other  terms  are 

1  am i nat es  then 


constants.  The  function  f(S/t)  is  constant  for  all 
a  plot  of  nt  22  versus  t/S  should  be  the  same  for  al  1  laminates. 
The  experimental  data  is  plotted  in  Fig.  F2  and  as  can  be  seen 
the  data  for  all  laminates  follow  the  same  trend  curve. 
Therefore,  it  can  be  seen  that  the  available  strain  energy 
release  rate  is  a  function  of  the  90°  layer  thickness  and  the 

matrix  crack  surface  area.  All  other  laminate  parameters  such  as 

>o 


the  number  of  consecutive  0'~'  plies  results 
effects  on  the  energy  release  rate. 


in  second  order 


rr.r.' 


rimental  strain  energy  release  rate  results  for  a  variety 
ross-ply  laminates.  (The  solid  lines  are  a  straight  line 
quadratic  least  squares  curve  fit) . 
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ABSTRACT 


The  development  of  damage  in  cross-ply  Hercules  ASV3502 

graphite/epoxy  laminates  has  been  investigated.  Specific  endeavors  were 

to  identify  the  mechanisms  for  initiation  and  growth  of  matrix  cracks 

and  to  determine  the  effect  of  matrix  cracking  on  the  stiffness  loss  in 

cross-ply  laminates.  Two  types  of  matrix  cracks  were  identified.  These 

include  both  "straight"  and  "angle"  type  cracks.  While  straight  cracks 

have  been  the  subject  of  extensive  research  efforts,  a  treatment  of  the 

angle  cracks  does  not  appear  in  the  literature.  The  growth  mechanisms 

of  the  angle  cracks  were  found  to  be  distinctly  different  from  those  of 

the  straight  cracks,  which  are  essentially  driven  by  shear  lag.  The 

experimental  study  of  matrix  crack  damage  revealed  that  the  angle  cracks 

formed  after  the  straight  cracks  and  followed  a  repeatable  pattern  of 

location  and  orientation  relative  to  the  straight  cracks.  Therefore,  it 

was  postulated  that  the  growth  mechanism  for  angle  crack3  is  driven  by 

the  3tress  state  resulting  from  the  formation  of  the  straight  cracks. 

This  phenomenon  was  analytically  investigated  by  a  finite  element  model 

of  straight  cracks  in  a  cross-ply  laminate.  The  finite  element  results 

substantiate  the  postulated  growth  mechanism.  Growth  mechanisms  for  both 

types  of  cracks  were  experimentally  related  to  the  number  of  consecutive 

90°  plies,  the  ratio  of  0°  to  90°  plies,  and  the  applied  stress  level. 

Finally,  experimental  results  indicate  that  as  much  as  a  10  percent 

degradation  in  axial  stiffness  due  to  matrix  cracking  in  cross-ply 

laminates  such  as  a  [0,90_]  . 

j  3 


INTRODUCTION 


Damage  accumulation  in  composite  laminates  has  become  an  extremely 
important  design  criterion  in  modern  aerospace  structures.  To  make  more 
efficient  use  of  composite  laminates  in  the  design  and  construction  of 
aircraft  and  large  space  structures,  the  scale  of  acceptable  damage  must 
be  increased.  Consequently,  the  composite  designer  must  be  aware  of  the 
effect  on  structural  behavior  of  damage  modes  3uch  as  matrix  cracking, 
interlaminar  delamination,  and  fiber  breakage.  Also,  he  must  be  able  to 
characterize  the  constitution  of  the  composite  laminate  throughout  its 
load  history.  In  spite  of  the  considerable  research  that  has  occurred, 
a  complete  understanding  and  formulation  of  the  mechanics  of  damage  in 
composite  materials  is  not  available. 

Most  of  the  previous  work  may  be  classified  as  experimental, 
emphasizing  the  various  types  of  damage  and  their  possible  causes  [1-8]. 
This  has  been  supported  by  extensive  analytical  investigations  which 
describe  the  stress  and  deformation  fields  in  the  damaged  regions. 
Methods  used  to  analyze  the  damage  regions  are  finite  elements  [9-1 6], 
finite  differences  [3,8,15],  shear  lag  [16],  modified  classical 
lamination  theory  [4,6,10],  and  general  analytical  techniques  ranging 
from  classical  fracture  mechanics  [17~24],  nonlinear  viscoelastic 
theories  [25],  to  general  boundary  value  problem  solutions  for  special 
crack  geometries  [18,26].  Only  recently  has  significant  attention  been 
given  to  the  development  of  cumulative  damage  models  capable  of 
predicting  the  resulting  damage  from  a  specified  load  history.  These 
models  proceed  from  three  basically  different  paths;  1)  fracture 
mechanics  [14,27_31],  2)  empirical  [32-39],  and  3)  internal  state 


variable  theories  which  are  based  on  thermodynamics  [40-49]. 

The  primary  objective  of  the  research  program,  for  which  the  results 
reported  herein  have  been  obtained,  is  the  development  of  a  damage 
dependent  constitutive  model  for  predicting  the  stiffness  loss  :n 
composite  laminates.  This  model  has  been  developed  from  a  continuum 
mechanics  approach  utilizing  thermodynamic  contraints  wherein  the  damage 
is  characterized  by  a  set  of  second  order  tensor  valued  internal  state 
variables  [48,49].  To  characterize  these  internal  state  variables,  the 
mechanics  of  damage  must  be  studied  experimentally  and  described 
analytically.  The  subject  of  this  paper  is  the  experimental  and 
analytical  correlation  of  matrix  crack  damage  in  continuous  fiber 
graphite/epoxy  laminates. 

After  a  brief  description  of  the  damage-dependent  constitutive 
model,  the  paper  will  focus  attention  on  the  following  research 
objectives: 

1 .  Experimental  study  of  the  mechanisms  of  initiation 
of  matrix  cracks  in  graphite/epoxy  laminates. 

2.  Correlation  of  experimental  observations  of  cracks 
with  finite  element  analysis  results. 

3.  Identification  of  the  effects  of  matrix  cracking  on 
the  axial  stiffness  loss  in  cross-ply  laminates. 


DAMAGE  MODEL 


A  continuum  mechanics  approach  has  been  utilized  by  the  author t  for 
predicting  the  stiffness  of  continuous  fiber  composites  subjected  to 
both  monotonic  and  cyclic  fatigue  loading  in  the  presence  of  damage 
[48,49].  In  this  model,  each  specific  damage  mode  is  characterized  by  a 
second  order  tensor  valued  internal  state  variable  representing  locally 
averaged  measures  of  damage.  Utilizing  thermodynamics,  the  constitutive 
equations  have  been  developed  as  follows: 


R  v\  n 

°LiJ  0  Lij  CLiJkl  eLkl  +  D  Lijkl  “  Lkl  , 


(1) 


where  oLi^  is  the  stress  tensor,  o  ^  represents  the  residual  stress 


tensor,  eLkl  is  the  strain  tensor,  CLijkl  is  the  undamaged  elastic 


modulus  tensor ,  D 


represents  the  damage  modulus  tensor,  and  a 


Lijkl  *  . . . . . . . . .  Lkl 

represents  the  internal  state  variables.  The  subscript  L  denotes 


locally  average  properties  with  all  other  subscripts  ranging  from  one  to 
three  and  the  superscript  n  ranges  from  one  to  the  number  of  damage 
modes . 


The  internal  state  variable  i3  represented  by  a  diadic  product 


between  the  crack  opening  displacement  vector,  u  ,  and  its 


corresponding  unit  normal  vector,  n  ^ ,  integrated  over  the  local  volume 


element  as  follows: 


Lkl 


{  [  < "? 
*  « 


dS 


(2) 


where  V.  is  the  local  volume  and  S  _  is  the  crack  surface  area  of  each 


damage  mode.  One  of  the  major  goals  of  the  research  is  to  establish  a 


£ 

H 
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direct  dependence  of  on  3ince  can  be  measured 

experimentally.  Complete  characterization  of  the  constitutive  equations 
will  require  development  of  the  growth  laws  for  the  internal  state 
variables  of  the  general  form 


ft 


u(V  V 


“l 


kl 


(3) 
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In  order  to  develop  accurate  growth  laws,  one  must  have  a  thorough 
understanding  of  the  microphysics  of  the  damage  process.  The  remainder 
of  this  paper  will  be  devoted  to  a  discussion  of  the  damage  mechanisms 
associated  with  the  first  fundamental  damage  mode  of  matrix  cracking. 


EXPERIMENTAL  PROGRAM 


The  purpose  of  the  experimental  program  was  twofold.  One  objective 
was  to  make  as  many  physical  observations  as  possible  concerning  the 
development  of  matrix  cracking  in  order  to  identify  the  mechanisms  of 
initiation  and  growth.  The  second  objective  was  to  identify  the  effect 
of  matrix  cracking  on  axial  stiffness  in  a  variety  of  angle-ply 
laminates.  To  isolate  this  specific  damage  state  cross-ply  laminates 
consisting  of  0°  and  90°  plies  were  studied  (3ee  Table  1).  This  array 
of  laminates  was  selected  in  order  to  examine  the  effect  of  the  number 
of  successive  90°  plies,  as  well  as  the  percentage  of  90°  plies  to  the 
0°  constraining  plies.  Another  reason  for  choosing  the  0°/90°  cross-ply 
laminate  was  that  it  minimizes  the  development  of  other  damage  modes, 
especially  edge  delamination.  The  material  system  was  Hercules  AS4/3502 
graphite/epoxy.  The  unidirectional  properties  are  given  in  Table  2. 

In  order  to  identify  the  various  aspects  of  matrix  cracking  one  must 
nondestructively  evaluate  the  damage  in  each  laminate.  This  was 
accomplished  by  using  edge  replication  [50]  and  X-ray  radiography.  These 
methods  have  become  quite  common  for  detecting  damage  in  composite 
laminates.  Since  graphite/epoxy  is  relatively  transparent  to  x-rays, 
zinc  iodide  [51  ]  was  introduced  into  the  damaged  laminate  to  produce 
enhanced  radiographs  of  the  damage. 

A  typical  edge  replica  showing  damage  in  a  cross-ply  [0,90,] 

5  3 

laminate  is  shown  in  Fig  1.  In  this  figure  one  can  see  two  distinctly 
different  types  of  matrix  cracks:  straight  cracks  and  angle  cracks.  The 
mechanisms  for  each  crack  will  be  discussed  in  the  next  section.  An 


SUCCESSIVE 
90°  PLIES 


PRECENTAGE 
9u  PLIES 


LAMINA  PROPERTIES 


LONGITUDINAL  MODULUS,  E1 

144.78  ±  2.0  %  GPa  (21.0  MSI) 

TRANSVERSE  MODULUS,  E22 

9.58  ±  2.1  %  GPa  (1.39  MSI) 

SHEAR  MODULUS,  G^1 

4.785  GPa  (0.694  MSI) 

POISSON'S  RATIO,  v12 

0.310  ±  3-7  % 

LONGITUDINAL  STRENGTH, F, 

lu 

2.247  ±  3.5  %  GPa  (326  KSI) 

TRANSVERSE  STRENGTH,  F„ 

tu 

0.0764  ±  9.8  %  GPa  (11.1  KSI) 

LONG.  FAILURE  STRAIN,  e, 

lu 

0.0144  ±  4.6  % 

TRAN.  FAILURE  STRAIN,  et 

tu 

0.00773  ±  6.7  % 

PREPREG  PROPERTIES 

FIBER  AREAL  WEIGHT2 

164  G/M2 

FIBER  VOLUME  FRACTION2 

0.58  ±  .03 

RESIN  CONTENT,  %  BY  WEIGHT2 

35  ±  3 

1  ANALYTICALLY  CALCULATED  FROM  A  [±45]  LAMINATE 

3 

2  HERCULES  INC.  SPECIFICATIONS 

Table  2.  Material  Properties  for  Hercules  AS4/3502. 


LOAD  DIRECTION 


Fig.  1.  EDGE  REPLICAS  OF  MATRIX  CRACKS; 

A)  PARTIAL  CRACKS,  B)  CURVED  CRACKS 
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x-ray  radiograph  of  this  same  damage  mode  is  shown  in  Fig  2.  It  is  very 
easy  to  distinguish  an  angle  crack  from  a  straight  crack.  A  straight 
crack  appears  as  a  sharp  narrow  band  while  an  angle  crack  appears  as  a 
wider,  fuzzy  band. 

The  mechanical  testing  was  conducted  using  an  Instron  1125  universal 
testing  system.  All  tests  were  conducted  under  monotonic  loading  at  a 
constant  cross-head  displacement  rate  of  0.05  in/min  following  ASTM 
standard  3039-76.  The  axial  strain  was  measured  using  a  one  inch  MTS 
extensometer( 632.1  IB-20).  The  data  were  recorded  with  a  digital  data 
acquisition  system  (at  35KHz)  using  a  Digital  1123  Plus  computer.  The 
loading  was  interrupted  periodically  to  measure  the  axial  stiffness  and 
to  document  the  damage  state  by  edge  replication  and  x-ray  radiography. 
A  typical  progression  of  damage  development  is  3hown  in  Fig.  3.  Other 
damage  modes,  including  interlaminar  delaminations  and  axial  splitting, 
are  visible  in  these  radiographs.  Interlaminar  delaminations  were 
observed  to  develop  at  the  Intersection  of  straight  and  angle  cracks 
with  the  0°  interface  (see  Fig.  4)  as  well  a3  at  the  intersection  of 
axial  splits  and  transverse  matrix  cracks.  (A  detailed  description  of 
the  experimental  methods  and  results  is  given  in  reference  [52].)  The 
remaining  portion  of  the  paper  will  be  concerned  only  with  matrix 
cracking.  Therefore,  all  mechanical  test  were  terminated  when  damage 
modes  other  than  matrix  cracking  initiatied. 
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STRAIGHT  CRACK 
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fig.  2,  detail  of  x-ray  radiograph  showing  difference  between 

STRAIGHT  AND  ANGLE  CRACKS 
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FIG.  3.  X-RAY  RADIOGRAPH  OF  DAMAGE  IN  A  [0/904lg  LAMINATE 
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EXPERIMENTAL  RESULTS 


IDENTIFICATION  OF  MATRIX  CRACK  MECHANISMS 

We  observed  in  Fig.  1,  two  distinct  types  of  matrix  cracks,  straight 

cracks  and  angle  cracks,  each  having  its  own  unique  growth  mechanism.  A 

schematic  diagram  of  a  well  defined  crack  system  for  a  [0,90  ] 

n  3 

laminate,  where  n  is  greater  than  1,  is  shown  is  Fig.  5.  A  straight 
crack  is  defined  to  be  a  matrix  crack  extending  completely  through  the 
off-axis  ply  oriented  normal  to  the  laminate  (x,y)  plane.  This  type  of 
crack  grows  very  suddenly  and  probably  emanates  from  a  preexisting 
microscopic  flaw.  This  behavior  is  generally  accepted  as  a  form  of 
brittle  fracture  and  can  be  described  by  fracture  mechanics  concepts.  It 
is  widely  accepted  that  the  density  (number  of  cracks  per  inch)  of 
straight  cracks  is  driven  by  shear  lag  [3]  between  the  off-axis  plies 
and  the  constraining  plies. 

An  angle  crack  is  a  matrix  crack  extending  partially  from  the  0°/90° 
interface  oriented  at  an  angle  to  the  laminate  plane.  Angle  cracks  are 
always  located  near  a  straight  crack  and  always  grow  toward  the  straight 
crack.  This  type  of  matrix  cracking,  though  probably  encountered  by 
others,  has  yet  to  be  reported  in  the  literature.  Sequential  edge 
replicas  indicate  that  the  growth  of  the  angle  crack  is  much  slower  than 
that  of  the  straight  cracks,  suggesting  a  ductile  mode  of  fracture. 
Hence,  this  would  represent  a  fracture  process  which  is  different  from 
that  for  straight  cracks.  The  angle  cracks  can  be  seen  to  extend  only 
partially  across  the  ply  thickness.  Angle  cracks  at  opposite  0°/90° 
interfaces  almost  always  coalesce  to  form  a  single  curved  crack,  as 
shown  in  Figs.  1  and  5.  Furthermore,  there  are  even  occurrences  of  two 


curved  cracks  located  on  one  side  of  a  straight  crack.  It  is 
hypothesized  that  the  orientation  and  location  of  the  curved  crack  is 
related  to  the  orientation  of  the  principal  stresses  in  the  vicinity  of 
a  straight  crack.  This  hypothesis  has  been  investigated  analytically  and 
will  be  discussed  in  a  later  section. 

The  experimented  results  showed  a  correlation  between  the  number  of 

curved  cracks  and  changes  in  the  local  stress  field  associated  with  the 

straight  cracks.  There  are  several  factors  that  influence  the  local 

stress  field  surrounding  the  straight  cracks.  One  factor  directly 

affecting  this  stress  field  is  the  number  of  consecutive  90°  plies  in 

the  laminate.  Figure  6  shows  the  relation  between  straight  cracks  and 

curved  cracks  for  am  increasing  number  of  90°  plies  in  a  [0,90  ] 

n  s 

laminate.  For  n-1  there  are  no  curved  cracks  because  most  of  the  stress 
is  carried  by  the  0°  ply.  As  n  increases  the  amount  of  stress  carried 
by  the  90°  plies  increases,  thus  resulting  in  a  corresponding  increase 
in  the  number  of  curved  cracks.  It  is  Interesting  to  note  that  for  n 
greater  than  2  there  aire  more  curved  cracks  than  straight  cracks, 
however  the  totad  number  of  cracks  for  all  three  laminates  is  48, 52, and 
50  respectively,  which  appears  to  be  fairly  constant. 


INCREASE  IN  CURVED  CRACKS  AS  NUMBER 
OF  CONSECUTIVE  90°  PLIES  INCREASES 
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Fig.  6.  Increase  in  Curved  Cracks  as  Number  of  Consecutive  90 
Plies  Increases. 
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EFFECT  OF  MATRIX  CRACKING  ON  AXIAL  STIFFNESS  LOSS 


The  axial  stiffness  loss  Jue  to  matrix  cracking  in  five  laminates 
was  monitored  during  step-wise  monotonic  loadings.  The  stiffness  was 
measured  only  on  the  unloading  portion  of  the  axial  stress-strain 
curves.  Thus,  stiffness  was  used  to  indicate  the  degree  of  damage 
development  during  any  monotonic  loading  step.  The  stiffness  results 
are  shown  in  Fig.  7.  All  curves  represent  an  average  of  two  or  more 
separate  specimen  results.  The  stiffness  values  are  normalized  to  the 
undamaged  stiffness  and  compared  to  the  total  crack  density,  which  is 
the  total  number  of  straight  cracks  and  curved  cracks  per  unit  length  of 
specimen. 

One  can  see  that  for  [0,90  ]  laminates,  where  n-1/2,  1,  2,  and  3, 

n  3 

stiffness  loss  increases  with  increasing  n.  This  is  explained  by  the 
fact  that  for  increasing  n  more  load  is  carried  by  the  90°  plies.  Hence, 
for  the  same  crack  density  (damage),  the  laminate  with  more  90°  plies 
will  lose  more  stiffness.  This  same  trend  can  be  explained  by  ply 
discount,  however,  the  ply  discount  method  over  estimates  the  total 
stiffness  loss.  It  is  noted  here  that  the  last  data  points  for  the 
[0,90,0]  and  [0,90]  were  taken  after  application  of  95  percent  of  the 
ultimate  load,  whereas  the  last  data  points  for  the  remaining  laminates 
correspond  to  the  test  termination  loads  when  other  damage  modes  began 
to  occur. 

For  laminates  with  identical  0°/90°  ratios  such  as  the  [0,90]  and 

s 

[0_,90_]  ,  the  laminate  with  the  higher  consecutive  number  of  90°  plies 

c  9 

shows  greater  stiffness  loss.  This  is  due  to  the  fact  that  the  laminate 
with  the  higher  consecutive  number  of  90' s  has  more  crack  surface  area 


Stiffness  Loss  in  Cross-Ply  Laminates  for  Matrix  Cracking  Only 


for  a  given  crack  density  resulting  in  a  larger  crack  opening 
displacment.  This  implies  that  the  internal  state  variable  an. .  in  the 

LK± 

constitutive  model  is  larger  for  increasing  values  of  consecutive  90° 
plies. 

A  final  trend  in  the  data  can  be  observed  by  combining  Figs.  6  and 
7.  The  stiffness  loss  increases  as  the  ratio  of  curved  cracks  to 
straight  cracks  increases.  This  result  suggests  that  the  curved  cracks 
cannot  be  considered  by  the  model  to  be  simply  additional  straight 
cracks. 


CORRELATION  OF  CURVED  CRACKS  WITH  FINITE  ELEMENT  RESULTS 


A  finite  element  model  utilizing  constant  strain  triangular  elements 
(CST)  incorporating  internal  boundary  generation  was  developed  for 
analyzing  the  stress  and  strain  fields  in  the  vicinity  of  a  straight 
matrix  crack.  The  ply  level  material  constitution  was  assumed  to  be 
transversely  isotropic  elastic.  The  objective  was  to  analytically 
correlate  the  observed  location  and  orientation  of  the  angle  crack  and 
to  investigate  the  influence  of  existing  straight  cracks  on  angle  crack 
formation. 

The  [0,90.,]  laminate  was  used  for  the  correlation  between  the 
experimental  and  analytical  results.  A  typical  crack  family  for  this 
laminate  is  shown  in  Fig.  8.  The  mean  experimental  distance  from  the 
straight  crack  to  the  nearest  curved  crack  at  the  ply  interface  (based 
on  25  separate  crack  measurements)  was  found  to  be  9.8  mils  (0.2489  mm), 
and  the  standard  deviation  was  found  to  be  2.08  mils  (0.05283  mm).  The 
minimum  distance  was  7.0  mils  (0.1778  mm).  The  mean  angle  between  the 
curved  cracks  and  the  normal  to  the  interface  was  found  to  be  30°.  The 
analytical  distance  measured  from  the  z-axis  was  determined  by  the 
location  of  the  maximum  principal  stress  relative  to  the  straight  crack. 
The  basic  finite  element  mesh  used  to  model  the  [0,90,]  laminate  with  a 
straight  matrix  crack  is  shown  is  Fig.  9.  The  squares  shown  in  the 
figure  actually  correspond  to  two  "CST"  elements.  The  broad  dark  line 
shown  between  the  0/90  interface  represents  the  resin  rich  region  where 
there  are  actually  four  rows  of  "CST"  elements.  More  elements  were  used 
in  the  interior  90°  ply  to  increase  the  accuracy  of  the  solution.  In 


all  plies,  the  aspect  ratio  per  element  was  kept  below  5.  The  crack 


density  was  varied  by  changing  the  length  of  the  overall  mesh.  The 
actual  mesh  used  corresponds  to  40  cracks  per  inch.  The  stress  field  was 
produced  by  applying  a  uniform  end  displacement  corresponding  to  a  net 
axial  strain  o^  one  percent. 

T.ie  orientation  and  magnitude  of  the  principal  stress  in  the  90°  ply 
is  significantly  influenced  by  the  properties  of  the  resin  rich  region 
since  the  majority  of  the  load  transferred  between  the  0/90  plies  occurs 
in  this  region.  The  actual  properties  of  the  resin  rich  region  were 
unknown  (partially  due  to  the  unavailability  of  neat  resin  properties), 
and  consequently  the  thickness  and  stiffness  properties  of  this  region 
had  to  be  assumed.  The  effects  of  these  assumed  properties  were 
parametrically  evaluated.  The  initial  thickness  of  the  resin  rich 
region  was  cnosen  to  be  1/10  the  ply  thickness,  where  one  ply  thickness 
was  0.005  inches  (0.127  mm).  The  constitutive  behavior  was  assumed  to 
be  isotropic  in  the  resin  rich  region,  with  E  -  0.5x10^  PSI  (3.4475  GPa) 
and  v  -  0 . 3 • 

The  analytical  results  for  the  above  properties  are  shown  in  Fig. 
10.  The  magnitude  of  the  principal  stress  for  the  first  three  rows  of 
elements  below  the  resin  rich  region  are  shown.  For  comparison,  Fig.  10 
also  shows  the  experimental  locations  of  the  curved  cracks.  Again,  it 
was  hypothesized  that  the  analytically  determined  location  of  the  curved 
crack  corresponds  to  the  peak  in  the  principal  stress.  Thus,  it  is 
found  for  this  case  that  the  peaks  do  not  fall  within  the  experimental 
range.  The  peaks  in  the  principal  stress  can  be  shifted  closer  to  the 
experimental  range  by  doubling  the  thickness  of  the  resin  rich  region  or 
by  decreasing  the  stiffness  of  the  resin  rich  region.  However,  the  peak 
in  row  1  does  not  fall  within  an  acceptable  distance  of  the  mean 
experimental  distance  using  reasonable  engineering  properties.  Matrix 


Fig.  10.  Variation  of  Principal  Stress  Near  the  0  /90  Interface 


yielding  in  the  resin  rich  area  is  another  factor  that  may  affect  the 
stress  distribution.  The  shear  stress  in  the  resin  rich  area  near  the 
straight  crack  is  much  larger  than  the  shear  strength  of  the  in-situ 
matrix.  ,<y  assuming  failure  in  those  elements  where  the  stress  exceeds 
the  shear  strength,  the  peak  in  the  principal  stress  can  be  shifted 
toward  the  experimental  results.  Finally,  the  inherent  numerical 
difficulties  associated  with  the  solution  of  this  problem  cannot  be 
ignored.  Thousands  of  "CST"  elements  are  required  to  adequately  model 
the  straight  crack  problem.  Furthermore,  the  "CST"  element  does  not 
model  large  gradients  in  shear  stress  very  accurately,  and  cannot 
account  for  the  large  singularities  in  the  3tress  field  near  the  crack 
tip. 

The  variation  in  the  principal  stress  angle  near  the  0/90  interface 
is  shown  in  Fig.  11.  The  experimental  results  are  also  shown  in  Fig. 

1 1  .  The  results  indicate  a  maximum  angle  of  about  40°  decreasing 
towards  0°.  The  principal  stress  angles  corresponding  to  the  peak 
principal  stresses  are  between  30°  and  40°.  This  is  in  close  agreement 
with  the  experimental  results.  The  complete  trajectory  (profile)  of 
angle  cracks  which  have  coalesced  to  form  "curved  cracks"  can  be 
estimated  by  examining  the  principal  stress  angle  throughout  the  90° 
plies  at  a  given  distance  from  the  straight  crack.  The  estimated 
profile  at  a  distance  of  0.004  inches  (0.1076  mm)  shown  in  Fig.  12,  has 
the  same  profile  as  the  curved  cracks  observed  experimentally.  While  it 
is  recognized  that  this  stress  state  will  be  altered  locally  once  the 
angle  crack  forms,  it  will  nonetheless  establish  the  general  direction 
of  crack  growth. 

For  completeness,  it  is  noted  that  similar  results  have  been 
obtained  for  the  analytical  location  and  orientation  of  the  angle  crack 
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ABSTRACT 

Experimental  evidence  has  shown  that  in  crossply  laminates  numerous 
matrix  cracks  often  become  curved  with  respect  to  the  plane  normal  to  the 
laminate  plane  when  three  or  more  90°  plies  are  stacked  in  succession. 
This  phenomenon  produces  damage  dependent  stiffness  losses  in  the  out-of¬ 
plane  direction.  The  authors  have  constructed  a  general  model  for 
predicting  stiffness  loss  as  a  function  op  damage  in  composite  laminates, 
and  the  current  paper  seeks  to  model  this  out-of-plane  stiffness  loss. 
The  model  predictions  are  compared  to  finite  element  results,  and  it  is 
shown  that  the  model  is  accurate  in  predicting  this  stiffness  loss  due  to 
curved  cracks. 

NOMENCLATURE 

Oj  -  ply  averaged  stress  in  Voigt  notation 
-  ply  averaged  strain  in  Voigt  notation 


nth  internal  state  variable  in  tensor  notation 


ctj  -  nth  internal  state  variable  in  extended  Voigt  notation 
V.  -  local  volume  element  over  which  properties  are  averaged 

L 


s!?  -  surface  area  of  cracks  in  nth  damage  mode  in  V. 

1 1  Li 

u.  -  crack  opening  displacement  vector 

n.  -  crack  '’ace  unit  normal  vector 
J 

C..  -  elastic  modulus  tensor  in  Voigt  notation 
1 J 

» 

3-m  -  damage  dependent  laminate  stiffness  component 

n  -  number  of  plies  in  laminate 

t,  -  thickness  o'-  nth  ply 
k 

o 

t -  -  midsurface  strain  components  of  laminate 

'j 


INT-tODUCTlCN 

Toe  authors  have  previously  developed  a  model  for  pred’ct’ag 
st'^fness  loss  as  a  function  of  damage  in  composite  media  [1].  The  model 
utilizes  a  sot  second  order  tensor  valued  internal  state  variables 
(ISVs)  to  describe  internal  damage.  These  ISVs  are  defined  by 


u.n  .dS 
i  J 


.  (1  ) 


Thermodynamic  concepts  have  been  utilized  to  construct  stress-strain 
relations  of  the  form 


o.  =  C.  .  e .  +  IV .  a.1 
1 J  J  tj  j 


where  the  above  equations  are  in  ply  coordinates  and  are  interpreted  to 
apply  at  the  ply  level.  Furthermore,  a11  is  contracted  to  single  index 
notation  using  tne  following  extended  Voigt  notation  [2]: 


C 

111 

n 

_  n 
=  a23 

n 

a? 

_  n 

=  a3l 

_  n 
=  a22 

0 

_  n 
-  a32 

a8 

E 

_  n 
=  a33 

< 

_  n 
=  “13 

*9 

■•J, 

.  (3) 


The  model  described  briefly  above  has  been  utilized  to  predict  axial 
stiffness  loss  as  a  function  of  matrix  crack  damage  in  several 
graphite/ epoxy  orossply  laminates  [2j.  Comparison  to  experiment  was  quite 
good.  The  model  has  also  been  utilized  to  perform  a  preliminary  study  of 
stiffness  loss  in  a  chopped  fiber  metal-matrix  composite  [3]. 

In  the  current  paper,  out-of-plane  stiffness  loss  is  modelled  for 
orossply  laminates  with  curved  cracks  of  the  type  shown  in  Fig.  1.  Since 
it  is  not  possible  to  obtain  accurate  experimental  measurements  of  this 
component  of  stiffness  loss,  the  model  is  compared  to  finite  element 
results  in  order  to  determine  model  accuracy. 

MODEL  FORMULATION 

In  reference  2  it  was  shown  that  equations  (2)  may  be  implemented  to 
a  lam'nate  analysis  scheme  to  produce  damage  dependent  stiffness 
components  in  global  laminate  coordinates.  The  resulting  stiffness 
equations  are 
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whore  the  overbars  indicate  that  the  quantities  are  in  the  laminate 
coordinates  shown  in  Fig.  2.  Therefore,  the  stiffness  predictions  require 
that  C.  .  and  a.  be  transformed  from  local  ply  coordinates  to  global 
laminatV  coordinates. 

The  last  term  in  equations  (4),  wherein  the  gradient  of  damage  Is 


taken  with 
conditions. 


st  term  in  equations  (4),  wherein  the  gradient  of  damage  is 
respect  to  midsurface  strain,  applies  only  to  unloading 


Experimental  evidence  indicates 
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graphite/epoxy  laminates  at  room  temperature  3ct./3e  are  to  a  first  order 
approximation  constant  during  unloading  [4].  "kurtPiermore ,  in  a  [0,90,0] 
laminate  all  components  of  3a./  3£  in  ply  coordinates  are  negligible 
except  This  term  md^y  bem  obtained  for  a  general  layup  by 
utilizing  experimentally  determined  energy  release  rates  from  a  single 
[0,90,0]  coupon  [2].  The  next  section  details  the  transformation  of  this 
quantity  from  the  coordinates  of  the  crack  to  laminate  coordinates. 


COORDINATE  TRANSFORMATIONS  FOR  CURVED  CRACKS 

Consider  now  a  local  volume  element  with  n  cracks  as  shown  in  Fig. 

c  0 

3.  Equations  (1)  may  be  written  in  the  following  form: 
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where  matrix  cracking  is  designated  by  the  superscript  1.  Now  define 
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so  that  equation  (9)  may  be  written 
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Now,  since  cx^  is  a  second  order  tensor, 
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where  a.  f  are  the  direction  cosines  relating  the  coordinates  of  the  kth 
crack  to?  the  laminate  coordinates.  Differentiating  (3)  with  respect  to 
the  midplane  strain  components  gives 


MATRIX  CRACKS 


LAMINATE  COORDINATE 
SYSTEM 


Matrix  Cracking  in  a  Laminated  Continuous  Fiber  Composite 
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transforming  to  the  coordinates  of  the  crack  gives 
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However,  since  it  is  assumed  that  3  a^/  3  £22  :i'3  the  on-Ly  non-zero 

component,  the  above  reduces  to 
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The  above  equations  may  be  utilized  to  obtain  the  last  term  in  reduced 
stiffness  equations  (ii),  where  it  is  assumed  that  30^,^,/ 3  is 

independent  of  crack  orientation  and  is  obtained  d  from  L0790.0] 
experimental  data  [2]. 
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ABSTRACT 


A  continuum  mechanics  approach  as  been  utilized  by  the  authors  to 
develop  a  model  for  predicting  the  stiffness  loss  of  continuous  fiber 
composites  subjected  to  both  monotonic  and  cyclic  fatigue  loading  in  the 
presence  of  matrix  cracking  and  interlaminar  delaminations.  In  this 
model,  each  specific  damage  made  is  characterized  by  a  second  order 
tensor  valued  internal  state  variable  representing  locally  averaged 
measures  of  both  matrix  cracks  and  interlaminar  delaminations.  The 
resulting  constitutive  equations  are  reduced  to  laminate  theory  by 
globally  averaging  the  stress  field  over  the  laminate  thickness.  In 
this  paper  the  internal  state  variables  are  defined  and  related  to  their 
respective  surface  area  of  damage.  For  simplicity  only  crossply 
laminates  [0,90,0]  ,  [0,90  ],  [0,90]g  and  t0„,902]g  are  studied.  It  is 
assumed  for  these  faminates  ^hat  only3  a  mode  i  fracture  process  need  be 
modelled  by  the  internal  state  variable. 

INTRODUCTION 

Two  common  forms  of  damage  in  continuous  fiber  composites  are 
matrix  cracks  and  interlaminar  delaminations.  Matrix  cracks  generally 
form  parallel  to  the  fibers  in  the  off-axis  plies  with  the  plane  of  the 
crack  extending  through  the  thickness  of  ply  [5,  51].  However,  the 
growth  through  each  ply  is  not  necessarily  self  similar  and  may,  in 
fact,  form  complex  patterns  such  as  curved  cracks  [50].  Interlaminar 
delaminations  usually  initiate  at  the  intersection  of  a  matrix  crack 
with  the  ply  interface.  These  delaminations  result  from  the 
interlaminar  stresses  which  develop  due  to  the  interaction  of  the  matrix 
crack  with  the  ply  interface.  A  small  but  significant  reduction  in  the 
stiffness  of  laminates  will  occur  due  to  these  damage  modes.  The 
stiffness  loss  is  an  important  parameter  because  it  has  been  related  to 
life  predictions  [5]. 

In  order  to  fully  utilize  the  potential  of  composites  one  must  be 
able  to  characterize  the  constitution  throughout  its  complete  load 
history.  With  regard  to  characterizing  the  constitution  of  composites 
with  damage,  there  has  been  a  considerable  amount  of  work.  This 
characterization  has  followed  three  basically  different  paths;  1  ) 
fracture  mechanics  [16,29-33];  2)  empirical  C 3 ^ 41  ) ;  and  3)  internal 
state  variable  theories  which  are  based  on  thermodynamics  {1,2,  42-49]. 

In  this  paper  the  internal  state  variable  theory  of  Allen,  et  al. 
[1,2],  is  utilized.  To  characterize  the  internal  state  variables,  the 
mechanics  of  damage  must  be  studied  experimentally  and  describe 
analytically.  The  subject  of  this  paper  is  to  characterize  the  damage 
model  and  compare  its  predictions  to  experimentally  determined  stiffness 
loss  in  crossply  laminates  with  both  matrix  cracks  and  interlaminar 
delaminations. 

DAMAGE  MODEL 


The  thermomechanical  constitutive  equations  with  distributed 
damage  as  developed  by  Allen,  et  al .  [1],  are  given  by 

(T  z  CrR  +  E  A 7 L  +  C 


where  n  ranges  from  one  to  two  for  each  damage  mode, 


are  the 
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residual  stresses,  e  ^  are  the  strains,  AT^  represents  the  temperature 
change  from  the  reference  state,  tensors L E,C ,  and  In  are  material 
constants,  and  ail  ,  are  the  locally  averaged  internal  state  variables. 
The  subscript  "L'^Ts  used  to  denote  locally  averaged  field  variables  and 
material  constants  and  for  notational  simplicity  will  be  omitted.  All 
other  subscripts  range  from  one  to  three.  In  this  research  all 
experimental  tests  will  be  conducted  at  constant  temperature. 
Furthermore,  it  will  be  assumed  that  the  chant  e  in  stress  due  to  each 
damage  mode  is  uncoupled,  thus  allowing  one  to  omit  the  cross  terms  on 
n.  It  will  also  be  assumed  that  the  residual  stresses  are  negligible  on 
unloading.  Under  the  above  assumptions  the  constitutive  equations 
simplify  to 
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where  for  notational  simplicity  the  superscript  "1"  on  I  has  been 
replaced  by  "m"  to  denote  matrix  cracking  and  the  superscript  "2"  on  I 
has  been  replaced  by  "D"  to  denote  interlaminar  delamination.  It  is 
noted  that  equations  (2)  represents  the  locally  averaged  constitutive 
equations  and  a  different  local  volume  must  be  used  for  each  damage 
mode.  Equations  (2)  are  reduced  to  single  index  notation  by 
incorporating  the  symmetry  of  the  stress  and  strain,  the  quadratic 
dependence  of  the  Helmholtz  free  energy  on  strain,  and  the  Voigt  single 
index  notation  [52].  The  resulting  equations  are 
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where  i  and  j  range  from  one  to  six  and  k  ranges  from  one  to  nine 
because  is  an  asymmetric  tensor 

The  constitutive  equations  can  be  simplified  further  by  imposing 
material  symmetry  constraints.  The  material  in  question  is  assumed  to 
be  initially  transversely  isotropic  in  the  undamaged  state.  It  is 
assumed  that  both  damage  modes  introduce  orthotropy  into  the  material 
system.  Therefore,  the  damage  tensor  I  J1  is  essentially  an  orthotropic 
tensor. 

In  order  to  complete  the  description  of  the  model  it  is  necessary 
to  characterize  the  internal  state  variables,  average  them  over  their 
respective  local  volumes,  and  then  introduce  the  constitutive  equations 
into  the  laminate  theory. 

Characterization  of  Internal  State  Variables 

The  locally  averaged  internal  state  variable  as  given  in  Reference 
[1 ]  is  defined  as  follows: 

<r  r  -L  f  «  drS  ,  (?, 


where  V.  is  the  local  volume,  S„  is  the  crack  surface  area  of  each 
u  c  2  G 

damage  mode,  u  is  the  crack  opening  displacement  vector,  and  n 

represents  the  unit  normal  vector  associated  with  the  crack  orientation. 


In  Reference  [2]  it  was  shown  that  for  predicting  the  axial 
stiffness  loss  in  crossply  laminates  due  to  matrix  cracking  that  the 
only  ISV  for  matrix  cracking  of  interest  is  Note  that  the  ply 

coordinate  system  used  here  is  identical  to  that  used  in  Reference  [53]* 
The  local  volume  for  this  damage  mode  is  a  single  ply.  Therefore,  all 
terms  associated  with  matrix  cracking  are  averaged  over  each  ply  before 
introduction  into  the  laminate  analysis  scheme. 

For  interlaminar  delaminations  there  are  onl-  th^ee  non  zero 
components  of  the  tensor  valued  internal  sta  ;e  variable,  aZ,  ot'T ,  and 
[54].  Since  only  the  axial  stiffness  loss  is  predicted  herein  only 
is  of  interest.  Furthermore,  since  interlaminar  delaminations  form  at 
the  ply  interface,  a  different  local  volume  must  be  used.  In  this  case 
the  entire  laminate  is  treated  as  the  local  volume. 


[2]. 


Note  that  it  can  be  shown  using  thermodynamics  that  I. .  ,  »  -C 
Therefore,  introducing  the  non  zero  internal  state  variables 


iM 


equation  (3)  will  reduce  to 


In  order  to  utilize  single  lamina  equations  to  characterize  the 
response  of  multilayered  laminates,  it  is  necessary  to  globally  average 
the  local  ply  constitutive  equations.  This  is  accomplished  herein  by 
imposing  the  Kirchhoff  hypothesis  for  plates  [53].  It  is  noted  that 
only  symmetric  layups  are  considered  and  that  only  the  resultant  forces 
are  desired.  The  resulting  equations  are  as  follows: 


where  [A]  represent^,  the  conventional  undamaged  extensional  stiffness, 
given  by  A.  -  ^  <Ci . ) k  t^,  {e}  are  the  midplane  strains,  {N™} 
represents  tne  resultant  K forces  due  to  matrix  cracking,  and  {N  } 
represents  the  resultant  forces  due  to  interlaminar  delamination.  It 
should  be  recalled  that  there  are  two  different  local  volumes  used  for 
averaging  the  properties  for  each  damage  mode.  For  matrix  cracking 


=  '  i  ( o,l  tK 


where  t  is  the  local  ply  thickness  nd  (a)k  represents  the  locally 
averagea  measure  of  the  ISV  per  ply.  Note,  equation  (7)  is  still  in 
local  ply  coordinates.  For  interlaminar  delamination  the  local  volume 
is  the  laminate,  therefore,  the  locally  averaged  properties  are  given  by 
the  respective  laminate  properties.  Thus,  the  resultant  forces  are 
given  by  _  n 


0 

A /■  = 


In  final  form  the  laminate  equations  are  thu3  given  by 


H=WW-4N>(^-W<-^ 


The  stiffness  in  the  axial  direction  is  defined  by 

s'  - 
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Substituting  equation  (10)  into  (11)  gives 
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For  crossply  laminates  with  matrix  cracking  in  only  the  off-axis  plies, 
equation  (12)  reduces  to  q  *\ 

where  p  denotes  the  number  of  0°  plies  and  q  denotes  the  number  of  90° 
plies.  Note,  this  equation  assume  the  same  damage  state  in  each  ply  and 
also  assumes  all  plies  the  same  thickness.  The  axial  modulus  is  thus 
given  by 
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Relation  of  Internal  State  Variable  to  Surface  Area 

The  internal  state  variables  can  be  related  to  the  surface  area  of 
damage  using  fracture  mechanics  and  thermodynamics.  Utilizing 
thermodynamic  constraints  it  was  shown  in  Allen,  et  al.  [1,2],  that  the 
free  energy  due  to  damage  reduces  to 

or  in  single  index  notation 
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From  fracture  mechanics  the  free  energy  due  to  damage  is  related  to  the 
strain  energy  release  rate  as  follows 
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where  is  the  surface  area  of  damage  for  each  mode.  It  is  assumed 
here  that  the  energy  release  rate  is  uncoupled  for  each  damage  mode. 
For  matrix  cracking,  then.  . 

rS*7/J 


r  ->  tv 
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and  for  interlaminar  delamination, 
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where  the  bar  over  u  denotes  a  laminate  free  energy.  Equating 
equations  (16),  (18),  and  (19)  and  noting  the  respective  local  volumes, 
the  internal  state  variables  are  related  to  the  surface  area  as  follows 


4c 
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The  energy  release  rate  for  each  damage  mode  can  be  defined  in  terms  of 
surface  area  by 


G  -ZX'. 

JS 
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where  the  conversion  S=2Ba  for  damage  surface  area  is  employed. 
Furthermore  the  strain  energy  in  the  body  is  given  by 

U  ~  ~  V  £<y  <fy  , 

Substituting  equation  (24)  into  equations  (22)  and  (23)  yields 


Therefore, 


Application  of  Internal  State  Variables 

The  constants  in  equations  (27)  and  (28)  must  be  determined 
experimentally.  In  Reference  [2]  it  was  shown  that  for  equation  (27) 
one  could  relate  the  strain  energy  release  rate  in  multiple  90°  plies  to 
the  energy  release  rate  in  one  90°  ply.  Thus,  one  needed  only  to 
conduct  experimental  test  on  one  laminate,  to  fully  characterize 
equation  (27)  for  any  multiple  stacking  sequence  of  90°  plies. 

In  work  to  date  if  has  been  found  that  for  equation  (28)  3  E  /3  SU 
for  any  crossply  laminate  is  directly  related  to  the  number  of 
constraining  0°  plies  and  only  indirectly  related  by  the  number  of 
consecutive  90°  plies.  Further  investigations  are  underway  to  analyze 
the  dependence  of  fracture  energy  for  delamination  on  the  laminate 
sequence. 


Model  Results 


At  this  time  no  experimental  data  have  been  obtained  for  the 
crossply  laminates  with  matrix  cracks  and  interlaminar  delaminations. 
However,  both  damage  modes  have  been  incorporated  into  a  (plane  strain^ 
finite  element  program  with  results  given  in  Fig.  1.  Valves  of  9e  /  9sU 
were  obtained  from  the  finite  element  results  and  are  given  below 


the  finite  element  results  shown  ir.  Fig.  1  are  produced  using  the  above 
values. 
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